VORTEX FORMATION AND DYNAMICS IN SUPERFLUID 3 He 
AND ANALOGIES IN QUANTUM FIELD THEORY 



V.B. Eltsov 

Low Temperature Laboratory, Helsinki University of Technology, Box 2200, FIN-02015 HUT, Finland 
and P.L. Kapitza Institute for Physical Problems, 119334 Moscow, Russia 

M. Krusius 

Low Temperature Laboratory, Helsinki University of Technology, Box 2200, FIN-02015 HUT, Finland 
G.E. Volovik 

Low Temperature Laboratory, Helsinki University of Technology, Box 2200, FIN-02015 HUT, Finland 
and L.D. Landau Institute for Theoretical Physics, 119334 Moscow, Russia 



The formation and dynamics of topological defects of different structure has been of central interest 
in the study of the 3 He superfiuids. Compared to superfiuid 4 He-II, the variability of the important 
parameters with temperature and pressure is wider and many features are more ideal. This has 
made experimentally new approaches possible. An example is the formation of quantized vortices 
and other topological defects in a rapid time-dependent phase transition from the normal state 
to 3 He-B. This vortex formation process is known as the Kibble-Zurek mechanism and is one the 
central topics of this review. It demonstrates the use 3 He-B as a model system for quantum fields, 
with detailed control from the laboratory. 

PACS: 67.57.Fg, 47.32, 05.70.Fh, 98.80.Cq. 

Keywords: superfiuid, quantized vortex, cosmic string, vortex formation, vortex dynamics, mutual 
friction, non-equilibrium phase transition, superfiuid turbulence, critical velocity. 



1 



CONTENTS 



1. Supcrfiuid 3 He and quantum field theory 3 

2. Defect formation in quench-cooled supcrfiuid transition 3 

2.1. Non-equilibrium phase transitions 3 

2.2. Cosmic large-scale structure 4 

2.3. Kibblc-Zurek mechanism 4 

2.4. Experimental verification of KZ mechanism 5 

2.5. Principle of supcrfiuid 3 He experiments 7 

2.5.1. Outline of experimental method 7 

2.5.2. Interpretation of 3 He experiments 8 

2.6. Measurement of vortex lines in 3 He-B 9 

2.6.1. Critical velocity of vortex formation 9 

2.6.2. Rotating states of the superfluid 10 

2.6.3. Experimental setup 10 

2.6.4. NMR measurement 11 

2.7. Vortex formation in neutron irradiation 12 

2.8. Volume or surface mechanism? 15 

2.9. Threshold velocity for vortex loop escape 17 

2.9.1. Properties of threshold velocity 17 

2.9.2. Influence of 3 Hc-A on the threshold 

velocity v cn 17 

2.10. Other defect structures formed in neutron 
irradiation 19 

2.10.1. Radiation-induced supercooled A — > B 
transition 19 

2.10.2. Vortex formation, AB interfaces, and KZ 
mechanism 20 

2.10.3. Spin-mass vortex 20 

2.11. Vortex formation in gamma radiation 23 

2.12. Bias dependence of loop extraction 24 

2.12.1. Experimental velocity dependence 24 

2.12.2. Analytic model of vortex loop escape 25 

2.13. Neutron-induced vortex formation at low 
temperatures 27 

2.13.1. Experimental techniques 28 

2.13.2. Measurement of vortex formation rate 29 

2.13.3. Supcrfiuid turbulence in neutron 

irradiation 30 

2.13.4. Calorimctry of vortex network 32 

2.14. Simulation of loop extraction 34 

2.14.1. Initial loop distribution 35 

2.14.2. Network evolution under scaling 
assumptions 36 

2.14.3. Direct simulation of network evolution 38 

2.15. Supcrfiuid transition as a moving phase front 40 

2.15.1. Neutron absorption and heating 40 

2.15.2. Thermal gradient and velocity of phase 

front 41 

2.16. Quench of infinite vortex tangle 43 

2.16.1. Vorticity on microscopic and macroscopic 
scales 43 

2.16.2. Scaling in equilibrium phase transitions 44 

2.16.3. Non-equilibrium phase transitions 45 

2.17. Implications of the quench-cooled experiments 45 

2.17.1. Topological-dcfect formation 45 

2.17.2. Phase transitions 46 

3. Vortex dynamics and quantum field theory analogues 47 

3.1. Three topological forces acting on a vortex and 

their analogues 47 

3.2. Iordanskii force 48 
3.2.1. Supcrfiuid vortex vs spinning cosmic 

string 48 



3.2.2. Gravitational Aharonov-Bohm effect 49 

3.2.3. Asymmetric cross section of scattering 

from a vortex 50 

3.2.4. Iordanskii force: quantized vortex and 
spinning string 50 

3.3. Spectral flow force and chiral anomaly 51 

3.3.1. Chiral anomaly 51 

3.3.2. Anomalous force acting on a continuous 
vortex and baryogenesis from textures 53 

3.3.3. Anomalous force acting on a singular 

vortex and baryogenesis with strings 54 

3.4. Analog of magnetogenesis: vortex textures 
generated in normal-supcrfluid counterflow 56 

3.5. Vortex mass: chiral fermions in strong magnetic 

field 58 

3.5.1. "Relativistic" mass of the vortex 58 

3.5.2. Contribution from bound states to the 

mass of a singular vortex 59 

3.5.3. Kopnin vortex mass in the continuous-core 
model: connection to chiral fermions in 
magnetic field 59 

3.5.4. Associated hydrodynamic mass 60 

3.5.5. Topology of the energy spectrum: gap 

nodes and their ramifications 61 
4. Final remarks 62 
Acknowledgements 62 
References 63 



2 



1. SUPERFLUID 3 He AND QUANTUM FIELD THEORY 

In recent times condensed matter and elementary par- 
ticle physics have been experiencing remarkable conver- 
gence in their developments, as many-body aspects have 
become increasingly more important in particle physics: 
The Early Universe is the ultimate application field for 
the theories on interacting particle systems which can 
be worked out for different energy regimes, or epochs 
of the expansion after the Big Bang. Compact astro- 
physical objects provide other "laboratories" with a nar- 
rower range of conditions in which to test theories. Ac- 
tual collider experiments are set up to study interacting 
particle systems, such as the quark-gluon plasma in heavy 
ion collisions or the pion condensate. 

Collective phenomena in interacting many-body sys- 
tems is what condensed matter physics is about, but el- 
ementary particle systems are adding to the picture ex- 
treme quantum behavior plus relativistic motion. So far 
systems, where all of these features would be of impor- 
tance, are not available for direct observation. During 
the nineties it turned out that valuable analogues can 
nevertheless be constructed from comparisons with non- 
relativistic many-body quantum systems of condensed 
matter, such as superconductors and helium superfluids. 
More recently optically cooled alkali atom clouds have 
been added to this list. Here we limit our discussion to 
the fermionic 3 He superfluids. 

The liquid 3 He phases provide attractive advantages 
as model systems for the study of various general con- 
cepts in quantum field theory: The liquid is composed 
of neutral particles in an inherently isotropic environ- 
ment — no complications arise from electrical charges 
or a symmetry constrained by a crystalline lattice. This 
allows one to concentrate on the consequences from a 
most complex symmetry breaking, which gives rise to 
a multi-dimensional order parameter space, but is well 
described by a detailed microscopic theory. Experimen- 
tally supcrfluid 3 He is devoid of extrinsic imperfections. 
In fact, with respect to impurities and dirt it is one of 
the purest of all condensed matter systems, excepting 
optically cooled atom clouds. The superfluid coherence 
length is large such that even surface roughness can be 
reduced sufficiently to transform the container walls to 
almost ideally behaving boundaries. There are topo- 
logically stable defects of different dimensionality and 
type, which often can be detected with NMR meth- 
ods with single-defect sensitivity. The study of the dif- 
ferent intrinsic mechanisms, by which these defects are 
formed, provide one of the important parallels to other 
systems. Phase transitions of both first and second order 
exist, which can be utilized in the investigation of de- 
fect formation. Another example is the use of the zero- 
temperature 3 He systems to model the complicated phys- 
ical vacuum of quantum field theory — the modern ether: 



The bosonic and fermionic excitations in 3 He (in the di- 
lute limit) are in many respects similar to the excita- 
tions of the physical vacuum — the elementary particles. 
This approach has been quite successful in constructing 
a physical picture of the interactions of elementary par- 
ticles with strings and domain walls, as described in a 
recent monograph by one of the present authors (Volovik 
2003). 

This review describes some examples from superfluid 
3 He research where strong connections exist to quantum 
field theory. In fact, some of these studies were only 
taken up because they can be used as laboratory models, 
to answer the question whether a proposed principle is 
physically correct or not. Our discussion has been split 
in two parts: The first part Chap. 2 deals with a de- 
tailed description of defect formation and evolution in 
a rapid non-equilibrium phase transition of second or- 
der, caused by the localized heating from an absorption 
event in ionizing radiation. This phenomenon originally 
became interesting because it can be thought to model 
phase transitions in the early Universe. Defect formation 
in these transitions has been suggested as the origin for 
the inhomogeneous mass distribution in the present Uni- 
verse. The second part Chap. 3 discusses concepts from 
vortex dynamics which have been used in Chap. 2 to an- 
alyze vortex formation in rapid phase transitions. This 
discussion proceeds in more general terms and continues 
to highlight connections to quantum field theory. 

2. DEFECT FORMATION IN QUENCH-COOLED 
SUPERFLUID TRANSITION 

2.1. Non-equilibrium phase transitions 

A rapid phase transition is generally associated with a 
large degree of inhomogencity. After all, this is the pro- 
cess by which materials like steel or amorphous solids are 
prepared. This disorder we attribute to heterogeneous 
extrinsic influence which is usually present in any system 
which we study: impurities, grain boundaries, and other 
defects depending on the particular system. To avoid dis- 
order and domain formation, we generally examine phase 
transitions in the adiabatic limit, as close to equilibrium 
as possible. 

But suppose we would have an ideally homogeneous 
infinite system with no boundaries. It is rapidly cool- 
ing from a symmetric high temperature state to a low 
temperature phase of lower symmetry, which we call the 
broken-symmetry phase, for example by uniform expan- 
sion. What would happen in such a transition? Are 
defects also formed in this case? Such a measurement 
of a homogeneous transition as a function of the transi- 
tion speed, is difficult. As we shall see, ultimately it also 
raises the question whether the source of the precipitated 
inhomogencity has been reliably identified. 
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2.2. Cosmic large-scale structure 

According to the standard theory of cosmology the 
Universe started off in the "Big Bang" in a homogeneous 
state. It then rapidly cooled through a sequence of phase 
transitions, in which the four fundamental forces of na- 
ture separated out. Today the Universe exists in a state 
with inhomogeneous large-scale structure. The clumped 
distribution of visible mass has become most evident from 
galaxy surveys — maps which show that galaxies form 
clusters and these in turn super clusters, such as the 
"Great Wall" , which are the largest structures discovered 
to date (Geller and Huchra 1989). The clustering takes 
the form of long chains or filaments, which are separated 
by large voids, regions empty of visible mass. Recent ex- 
tended galaxy surveys indicate that the length scale of 
large-scale structure is of order 100 Mpc 1 (Peacock et al. 
2001). 

Another image of large-scale structure has been pre- 
served in the cosmic microwave background radiation, 
from the time when the Universe had cooled to a few eV, 
when nuclei and electrons combined to form atoms and 
the Universe became transparent to photons. Since then 
the background radiation has cooled in the expanding 
Universe. It now matches to within 3 parts in 10 the 
spectrum of a black body at 2.728 K, as measured eg. 
with the spectrometers on board of the satellite Cosmic 
Background Explorer (COBE) in 1990. Later balloon- 
borne measurements examined the spatial distribution 
of the residual anisotropy with high angular resolution 
and amplitude sensitivity. The recent satellite Wilkinson 
Microwave Anisotropy Probe (WMAP) has extended this 
work, mapping all-sky surveys with an angular resolution 
better than 1° on different frequency bands, to separate 
out the interfering signal from our galaxy (Bennett et al. 
2003; Spergel et al. 2003). Combined with other infor- 
mation, this anisotropy of the cosmic microwave back- 
ground radiation, 30 /xK in amplitude, which represents 
the density fluctuations in the structure of the Universe 
when it was only 300 000 years old, is expected to ex- 
plain the mechanisms which originally seeded and led to 
the formation of the large-scale structure, as we observe 
it today. 



2.3. Kibble-Zurek mechanism 

One of the early explanations for the origin of large- 
scale structure was offered by Tom Kibble (1976). He 
suggested that the inhomogeneity was created in rapid 
phase transitions of the early expanding Universe. Even 



in a perfectly homogeneous transition of second order de- 
fects can be expected to form, if the transition proceeds 
faster than the order parameter of the broken-symmetry 
phase is able to relax. In such a non-equilibrium transi- 
tion the new low-temperature phase starts to form, due 
to fluctuations of the order parameter, simultaneously 
and independently in many parts of the system. Sub- 
sequently during further cooling, these regions grow to- 
gether to form the new broken-symmetry phase. At the 
boundaries, where different causally disconnected regions 
meet, the order parameter does not necessarily match 
and a domain structure is formed. 

If the broken symmetry is the U(l) gauge symmetry, 
then domains with different value of the order parame- 
ter phase are formed. Such a random domain structure 
reduces to a network of linear defects, which are vor- 
tex lines in superfluids and superconductors or perhaps 
cosmic strings in the Early Universe. If the symmetry 
break is more complicated, as is the case in 3 He superflu- 
ids, then defects of different dimensionality and structure 
may be formed. 

Subsequent numerical simulations of rapidly cooled 
second order phase transitions confirmed that defects 
were indeed forming (Vachaspati and Vilenkin 1984). 
In 1985 Wojciech Zurek proposed a conceptually pow- 
erful phcnomcnological approach how to understand a 
phase transition far out of equilibrium, when the out- 
come from the transition becomes time dependent (Zurek 
1985, 1996). He characterizes the transition speed with 
a quench time 



T Q = 



dT 



df 



(1) 



T=T C 



which allows him to approximate temperature and time 
with a linear dependence during the thermal quench. 

The quench time tq is compared to the order param- 
eter relaxation time r, which in a Ginzburg-Landau sys- 
tem at a second order phase transition is of the general 
form 



r(T) = t (1 - T/T c y 



(2) 



1 1 Mpc = 10 6 parsec = 3.262-10 6 light years. 



In superfluid 3 He, to is on the order of tq ~ £o/vf, 
where £o is the zero temperature limiting value of the 
temperature (T) and pressure (P) dependent superfluid 
coherence length £(T,P). Close to T c , it is of the form 
£(T,P) = £o(-P)(l - T/Tc)- 1 / 2 . The second quantity, 
up, is the velocity of the thermal quasiparticle excita- 
tions which are excited above the superfluid energy gap. 

As sketched in Fig. 1, this means that below T c the 
order parameter coherence spreads out with the veloc- 
ity c(T) ~ £/t - Co (1 - T/Tc) 1 / 2 /^. The freeze-out of 
defects occurs at t — izurck, when the causally discon- 
nected regions have grown together and superfluid co- 
herence becomes established in the whole volume. At the 
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FIG. 1. Principle of KZ mechanism. Rapid thermal quench 
through T c : (Left) Temperature T(t) and its linear approxi- 
mation T — T c (l — i/ r Q) during the quench, as a function of 
time t. (Right) Superfluid coherence length £(T) and order 
parameter relaxation time r(T) diverge at T c . At the freeze- 
out point tzurok, when phase equilibrium is achieved, the edge 
of the correlated region, the causal horizon, has moved out to 
a distance £H(izurok), which has to equal the coherence length 

£(tZurck)- 

freeze-out temperature T(tzurck) < T c , the causal horizon 
has travelled the distance ^H^Zurck) = J Zurek c(T)dt = 
CoT"q(1 — 7z ure k/T c ) 3 / 2 /to which has to be equal to the 
coherence length £(tzurck)- This condition establishes the 
freeze-out temperature Tz ure k/T c = 1 — ^tq/tq at the 
freeze-out time tzurck = ^/toTq , when the domain size 
has reached the value 

ev-e H (tZurck)=?0 (tq/tq) 1 / 4 . (3) 

In superfluid 3 He we may have £o ~ 20 nm, To <~ 1 ns, 
and in the best cases a cool-down time of tq <~ 1 ^s 
can be reached. From these values we expect the domain 
structure to display a characteristic length scale of or- 
der £ v ~ 0.1 /um. Assuming a U(l)-symmetry-breaking 
transition, where vortex lines are formed at the domain 
boundaries, the average inter-vortex distance and radius 
of curvature in the randomly organized vortex network 
is on the order of the domain size £ v . In general, we can 
expect a rapid quench to lead to an initial defect density 
(defined as vortex length per unit volume) 



where the numerical factor a; <~ 1 - 100 depends on the 
details of the model system. Numerical simulation calcu- 
lations (Ruutu et al. 1998a) give ai w 2 for 3 He-B. In this 
system with global symmetry and negligible thermal fluc- 
tuations the Kibblc-Zurek (KZ) model of defect forma- 
tion can be expected to work well. Moreover, Ginzburg- 
Landau theory provides a good description close to T c 
and might lead via microscopic calculations to more rig- 
orous predictions even on a time-dependent phase tran- 
sition. 



2.4. Experimental verification of KZ mechanism 

The predictions of the KZ model have been tested in 
numerous numerical experiments under varying condi- 
tions (cf. Bettencourt et al. 2000). Laboratory experi- 
ments are fewer and their results are less conclusive. 

The first experiments were performed on liquid crys- 
tals (Chuang et al. 1991; Bowick et al. 1994), where the 
transition from the isotropic to the nematic phase was ex- 
amined. This is a weakly first order transition of rod-like 
molecules from a disordered to an ordered state, where 
defects of different dimensionalities can be formed. The 
results were generally found to be consistent with the 
KZ interpretation. However, a liquid-crystal transition 
at room temperature is rather different from the ideal 
KZ case. 

The first experiments on coherent quantum systems 
were started with liquid 4 He (Hendry et al. 1994), by 
releasing the pressure on a convoluted phosphor-bronze 
bellows so that the sample expanded from the normal 
liquid through the A transition into the superfluid phase. 
For a mechanical pressure system the quench time is of 
order tq ~ 10 ms. The depressurization is followed by a 
dead time of similar length, which is required for damping 
down the vibrations after the mechanical shock. The 
vortex density after the transition is determined from the 
attenuation in second sound propagation through bulk 
liquid. Expansion of 4 He across the T\(P) line (which 
has negative slope in the (P, T) phase diagram) requires 
careful elimination of all residual flow which might result 
from the expansion itself, foremost the flow out of the 
filling capillary of the sample cell, but also that around 
the convolutions of the bellows. Such flow will inevitably 
exceed the critical velocities, which approach zero at T\, 
and leads to vortex nucleation (Sec. 2.6.1). After the 
initial attempts, where vortices were detected but a fill 
line was still present, later much improved measurements 
failed to yield any evidence of vortex line production, 
in apparent contradiction with the KZ prediction (Dodd 
et al. 1998). 

The early 4 He expansion experiments were analyzed 
by Gill and Kibble (1996) who concluded that the ini- 
tially detected vortex densities were unreasonably large 
and must have originated from extrinsic effects, most 
likely created by the flow out the fill line. The absence 
of vorticity in the later results was shown not to be evi- 
dence against the KZ mechanism (Karra and Rivers 1998; 
Rivers 2000): Although the pressure drop in expansion 
cooling is large, 20 - 30 bar, the change in relative tem- 
perature e = 1 - T/T\ is small, |Ae| < 0.1. The final 
temperature £/ < 0.06 is below T\ but still above the 
Ginzburg temperature Tq, ie. in the regime where ther- 
mal fluctuations might be sufficiently effective to wipe 
out domain formation. Here Tq is defined as the tem- 
perature at which U^Tq equals the energy of a vortex 
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loop of the size £ v > which is expected after the quench 
according to Eq. (3). Thus the current view about the 
4 He expansion experiments holds that the vortices decay 
away before the observation switches on (Rivers 2000; 
Hendry et al. 2000). 

Simultaneously with the 4 He expansion experiments, 
measurements of quite different nature were carried out 
in 3 He-B (Ruutu et al. 1996a; Baucrle et al. 1996). 
Here a detailed comparison with the KZ predictions was 
achieved. These efforts were followed by measurements 
of thermal quenches through the superconducting tran- 
sition in different superconducting geometries. 

The most conclusive results have been obtained with 
high-T c YBCO films (Maniv et al. 2003). The film is 
heated above its transition temperature of 90 K with a 
laser pulse. The hot film is then cooled from below back 
to the ambient temperature of 77 K. This happens lat- 
erally uniformly through the contact with an optically 
transparent substrate of much larger thickness and heat 
capacity. The maximum cooling rate was measured to 
be 10 8 K/s, by monitoring the electrical resistance of the 
film. The total net flux, which remains frozen in the 
pinning sites of the YBCO film after the transition, was 
measured with a sensing coil below the substrate and con- 
nected to a SQUID magnetometer. This flux was found 
to display very weak dependence on the cooling rate, as 
predicted by the KZ scaling theory. 

As discussed by Kibble and Rajantie (2003), the im- 
portant information to obtain from such measurements 
would be the spatial distribution of the flux quanta in the 
film. Spatial correlations between the fluxoids would sup- 
port the KZ mechanism, if the vortices were found to dis- 
play negative correlations, ie. neighboring vortices would 
prefer to have opposite signs. A competing mechanism 
involves thermal fluctuations of the magnetic field which 
generates vortices during the quench cooling through T c 
where the critical field vanishes: H C 2 — > 0. This second 
mechanism is predicted to display positive correlations, 
ic. vortices of same sign would tend to cluster. An analo- 
gous situation arises in the 3 He-B experiments which will 
be described below. Here at temperatures below 3 mK 
thermal fluctuations are not important since any energy 
barriers will generally automatically be orders of mag- 
nitude higher. Nevertheless, even in these conditions a 
phenomenon can be identified which competes with the 
KZ mechanism (which is a volume effect). This involves 
the phase boundary between the normal liquid and 3 Hc-B 
(a surface effect) if a phase boundary exists in the pres- 
ence of flow. The reason is again that at T c the critical 
velocity for vortex formation vanishes: v c — > 0. 

In superconductors a measurement with spatial reso- 
lution and single-fluxoid sensitivity has been recently ac- 
complished (Kirtley et al. 2003). Here a scanning SQUID 
microprobe is used to map the flux quanta trapped in an 
amorphous MosSi film which had been lithographically 



patterned in ring-like elements. Unfortunately, in this 
measurement the cooling rate was limited to <~ 20 K/s 
at the T c = 7.8 K of the film and the detected quan- 
tized flux was concluded to originate from a thermally 
activated process. 

These examples illustrate the difficulties in attempting 
to find experimental proof for the KZ process. The KZ 
mechanism is intuitively easy to accept, but its exper- 
imental verification is complicated and involves careful 
considerations. In a phase transition defect formation is 
more usual than its absence. In most experiments both 
extrinsic and intrinsic sources of vortex formation are 
present and compete in actual vortex generation. For 
instance, even slow cooling of liquid 4 He across the A 
transition is known to produce primordial vortices which 
remain pinned as remanent vorticity on the walls of the 
container (Awschalom and Schwarz 1984). There are no 
claims in the literature yet that an experiment would 
have produced a perfectly vortex-free bulk sample of su- 
pcrfluid 4 He. 

A second point to note is that the KZ model describes 
a second order transition where the energy barrier, which 
separates the symmetric high-temperature phase from 
the broken-symmetry states at low-temperatures, van- 
ishes at T c . Thus the transition becomes an instability. 
A first order transition is different: Here the barrier re- 
mains finite and the low temperature phase has to be 
nucleated, usually by overcoming the barrier via thermal 
activation. Nucleation may occur in different parts of the 
system nearly simultaneously, depending on the proper- 
ties of the system, the fluctuations, the quench time tq, 
and the nucleation mechanism. Such a situation, which 
often is called bubble nucleation, also leads to domain 
formation and to a final state which is qualitatively simi- 
lar to that expected after the second order transition and 
the KZ mechanism. 

In all of the experimental examples listed above, other 
mechanisms are also possible which might replace or op- 
erate in parallel to the KZ process. The differences are 
subtle, as we shall see in Sec. 2.8 for 3 He-B or as has been 
emphasized for superconductors by Rajantie (2001). In 
practice measurements differ from the original KZ model, 
ie. from a homogeneous transition, which occurs simul- 
taneously in the whole system and where there are no 
gradients in density or temperature. In practical mea- 
surements on bulk material a rapid quench through the 
transition is often forced by strong gradients. On some 
level gradients will always appear in any laboratory sys- 
tem and their influence has to be investigated. 

Superfluid 3 He measurements have so far established 
the most quantitative comparison to the KZ scaling the- 
ory (Ruutu et al. 1996a; Bauerle et al. 1996). A happy 
coincidence of many valuable features has made this pos- 
sible: A reduced likelihood of defect formation from ex- 
trinsic sources, the virtual absence of remanent vorticity, 
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the possibility to perform thermal quenches in microsec- 
onds, the presence of a pure second order phase transi- 
tion, and the availability of measuring methods to detect 
different types of defects, often with a resolution of one 
single defect. A critical discussion of these questions is 
presented below. 

The search for firm proof of the KZ mechanism has 
grown to an important topic in condensed matter physics. 
The reason is that it provides (at least superficially) easy 
interpretation of the fastest non-equilibrium phase tran- 
sitions. Amusingly this question does not enjoy simi- 
lar interest any more in cosmology where the problem 
was born. Measurements of the anisotropy in the cos- 
mic microwave background radiation with the WMAP 
satellite (with an angular resolution <~ 0.25° and an am- 
plitude sensitivity <~ 1 /iK) have concluded that the an- 
gular distribution of the anisotropy expanded in its har- 
monic components gives a multipole expansion with a 
sharp dominating peak at £ ~ 200, which corresponds 
to an angular scale <~ 1° (Bennett ct al. 2003; Spcrgcl 
et al. 2003). This distribution can be fit with models 
which are based on inflationary expansion of the Early 
Universe, i.e. with models with an early period of accel- 
erated expansion. In contrast, topological defects, such 
as cosmic strings, are expected to generate a clearly dif- 
ferent broader peak. These two different scenarios have 
been the ruling contenders as the explanation for the ori- 
gin of large-scale structure. Thus at present time the KZ 
mechanism appears to survive only in condensed matter 
physics. 

2.5. Principle of superfluid 3 He experiments 

To study reliably defect formation in a quench-cooled 
transition, two basic requirements are the following: 
First, the transition can be repeated reproducibly, and 
second, a measurement is required to detect the defects 
after their formation, either before they annihilate or by 
stabilizing their presence with a bias field. Both of these 
requirements can be satisfied in 3 He experiments. 

2.5.1. Outline of experimental method 

A schematic illustration of one set of experiments (Ru- 
utu et al. 1996a) is depicted in Fig. 2. The top part dis- 
plays the events which are controlled by the observer in 
the laboratory. A cylindrical container filled with super- 
fluid 3 He-B is rotated at constant velocity il. The veloc- 
ity of rotation is maintained below the critical value at 
which quantized vortex lines are spontaneously formed. 
In other words, the initial state is one of metastable 
vortex- free counterflow of the superfluid and normal com- 
ponents in constant conditions. 




Neutron absorption Heated bubble Vortex network 

escape 

FIG. 2. Principle of the quench-cooled rotating experiment in 
superfluid 3 He-B: (Top) A cylindrical sample container with 
superfluid He-B is rotated at constant angular velocity Q 
and temperature T, while the NMR absorption is monitored 
continuously. When the sample is irradiated with neutrons, 
vortex lines are observed to form. (Bottom) Interpretation 
of the processes following a neutron absorption event in bulk 
superflow; see the text for details. (From Ruutu et al. 1996a). 

Next a weak source of thermal neutrons is placed in 
the vicinity of the 3 He sample. If the rotation velocity is 
sufficiently high, then vortices start to appear at a rate 
which is proportional to the neutron flux. The neutron 
source is positioned at a convenient distance d from the 
cryostat so that vortex lines are observed to form in well- 
resolved individual events. The experimental signal for 
the appearance of a new vortex line is an abrupt jump in 
NMR absorption. An example, measured at low rotation, 
is shown by the lower recorder trace in Fig. 3. 

The lower part of Fig. 2 shows in more detail what is 
thought to happen within the superfluid. Liquid 3 He can 
be locally heated with the absorption reaction of a ther- 
mal neutron: n + ^Hc — > p + fH + E , where E = 764 
keV. The reaction products, a proton and a triton, arc 
stopped by the liquid and produce two collinear ioniza- 
tion tracks (Meyer and Sloan 1997). The ionized parti- 
cles, electrons and 3 He ions, diffuse in the liquid to recom- 
bine such that 80 % or more of E is spent to heat a small 
volume with a radius R b ~ 50 ^m from the superfluid 
into the normal phase. The rest of the reaction energy es- 
capes in the form of ultraviolet emission (Stockton et al. 
1971; Adams et al. 1995) and, possibly, in the form of 
long-living molecular 3 He2 excitations (Keto et al. 1974; 
Kafanov et al. 2000), which partially relax at the walls 
of the container. 
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FIG. 3. Peak height of the NMR signal which monitors the 
velocity of the bias flow from rotation (the so-called counter- 
flow peak in the NMR absorption spectrum of He-B), shown 
as a function of time t during neutron irradiation. The lower 
trace is for a low and the upper trace for a high value of rota- 
tion velocity Q. The initial state in both cases is metastable 
vortcx-free counterflow. The vertical arrows indicate when 
the neutron source was turned on/off. Each step corresponds 
to one neutron absorption event and its height, when com- 
pared to the adjacent vertical calibration bar, gives the num- 
ber of vortex lines formed per event. The distance of the 
neutron source from the sample is denoted with d. (From 
Ruutu ct al. 1996b). 

For measurements in rotating 3 He, neutron absorption 
is an ideal heating process. For slow neutrons from a 
room temperature source with a Maxwellian distribution 
(at 20.4 °C or 0.0253 eV) the velocity at the peak of the 
distribution is 2200 m/s. The capture cross section of the 
3 He nucleus for these neutrons is huge: 5327 ± 10 barn = 
5.3- 10~ 21 cm 2 (Bcckurts and Wirtz 1964). This means a 
short mean free path of about 100 /im. Thus most neu- 
trons are absorbed within a close distance from the wall, 
although still within the bulk liquid, but such that the ve- 
locity of the applied counterflow is accurately specified. 
In other parts of the refrigerator neutron absorption is 
negligible and thus the temperature of the 3 He sample 
can be maintained stable during the irradiation. Finally, 
the cooling is fast, as the small volume heated by the 
neutron absorption event is embedded in the cold bulk 
system. 

Subsequently, the normal liquid in the heated neutron 
bubble cools back through T c in microseconds. The mea- 
surements in Fig. 3 demonstrate that vortex lines are 
indeed created in this process. Vortex loops, which form 
within a cooling neutron bubble in the bulk superfluid, 
would normally contract and disappear in the absence of 
rotation. This decay is brought about by the inter-vortex 
interactions between the vortex filaments in the presence 
of dissipation from mutual friction between the normal 
and superfluid components. In rotation the externally 



applied counterflow provides a bias which causes suffi- 
ciently large loops to expand to rectilinear vortex lines 
and thus maintains them for later detection. In the rotat- 
ing experiments in Helsinki these rectilinear vortex lines 
are counted with NMR methods (Ruutu et al. 1996a, 
1998a; Finnc ct al. 2004a). 

In a second series of 3 He experiments, performed in 
Grenoble (Bauerle et al. 1996, 1998a) and Lancaster 
(Bradley et al. 1998), the vortices formed in a neutron ab- 
sorption event are detected calorimctricly with very-low- 
temperature techniques. In the zero temperature limit 
mutual friction becomes exponentially small and the life 
time of the vorticity might be very long, even if there is 
no flow. In this situation the existence of the vorticity 
can be resolved as a deficit in the energy balance of the 
neutron absorption reaction. 

2.5.2. Interpretation of 3 He experiments 

Several alternative suggestions can be offered on how 
to explain these observations. The KZ model is one of 
them. Whatever the mechanism, the superfluid transi- 
tion in a localized bubble represents a new process for 
creating vortex rings. This phenomenon is associated 
with one of the fastest 2nd order phase transition, which 
has been probed. It takes place in the bulk superfluid and 
not in immediate contact with a solid wall, where most 
other vortex formation processes occur in the presence of 
applied counterflow. 

The KZ interpretation of the later stages of the exper- 
iment is contained in the two illustrations marked as (c) 
and (d) in Fig. 2. First a random vortex network (c) 
is created while cooling through T c . Next the network 
starts to evolve under the influence of inter-vortex inter- 
actions and the applied flow. At high temperatures the 
process is highly dissipative, owing to the large mutual 
friction force which connects the superfluid and normal 
fractions and acts on a moving vortex. The average inter- 
vortex distance increases, small loops are smoothed out, 
and the network becomes rarefied or "coarse grained". 
Also reconnections, which take place when two lines 
cross, contribute to the rarefaction. The applied flow 
favors the growth of loops with the right winding direc- 
tion and orientation. It causes sufficiently large loops to 
expand, while others contract or reorient themselves with 
respect to the flow. The final outcome is that correctly 
oriented large loops (d) , which exceed a critical threshold 
size, start expanding spontaneously as vortex rings, un- 
til they meet the chamber walls. There the superfluous 
sections of the ring will annihilate and only a rectilinear 
vortex line will finally be left over. It is pulled to the 
center of the container, where it remains in stationary 
state parallel to the rotation axis, stretched between the 
top and bottom surfaces. 
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This picture applies at high temperatures in 3 Hc-B 
where Kelvin-wave excitations on vortices are exponen- 
tially clamped and spontaneous loop formation on exist- 
ing vortices does therefore not occur. In this situation the 
number of vortex rings, which are extracted from the neu- 
tron bubble (Fig. 2d), remains conserved during the later 
expansion. Thus the number of rectilinear vortex lines, 
the end result from the neutron absorption event, charac- 
terizes the extraction process as a function of the applied 
bias velocity. This situation is very different from that in 
4 Hc-II where mutual friction damping never reaches this 
high values (except extremely close to T\). 

Very different time scales are at work in this process 
of vortex formation. Except for the retarded molecular 
excitations, the heating from the initial ionization and 
subsequent recombination is limited by the diffusion of 
charges in the liquid and takes place much faster than 
the thermal recovery. The initial vortex network forms 
during cooling through T c , for which the relevant time 
scale is microseconds. The later evolution of the network 
and the loop escape happen again on a much slower time 
scale, namely from milliseconds up to seconds, since here 
the vortex motion is governed by the mutual-friction- 
dependent superfluid hydrodynamics. 

2.6. Measurement of vortex lines in 3 He-B 

There are two major phases of superfluid 3 He, the 
A and B phases (Fig. 14). The neutron measurements 
have been performed in the quasi-isotropic 3 He-B. In the 
present context we may think of vortices in 3 He-B as 
being similar to those in superfluid 4 Hc-II, where only 
the U(l) symmetry is broken and the order parameter 
is of the general form * = |*(T)| e i$(r) : The vortices 
are topologically stable, have a singular core, inside of 
which \^f\ deviates from its bulk value, while outside the 
phase <I> changes by 2-kv on a closed path which encir- 
cles the core. A persistent superfluid current is trapped 
as a single-quantum circulation [y = 1) around the core, 
with the circulation quantum k = h/{2mz) — 0.0661 
mm 2 /s and a superflow velocity v S!Vortox = «/ {2ir) V$ = 
«/(27rr) <fi. 

2.6.1. Critical velocity of vortex formation 

The energy of a single rectilinear vortex line, aligned 
along the symmetry axis of a rotating cylindrical con- 
tainer of radius R and height Z, consists primarily of 
the hydrodynamic kinetic energy stored in the superfluid 
circulation trapped around the core, 

P 1 f 2 , T/ PsK 2 Z R 

v = 2j psV ^°^ dv = ln j ■ ( 5 > 



Here the logarithmic ultraviolet divergence has been cut 
off with the core radius, which has been approximated 
with the coherence length £. If the container rotates 
at an angular velocity ft, the state with the first vortex 
line becomes energetically preferred when the free energy 
E v — QL Z becomes negative. The hydrodynamic angu- 
lar momentum from the superflow circulating around the 
vortex core is given by 

L z = J p s rv^ VOItcx dV = p s nR 2 Z . (6) 

One finds that E v — Q,L Z < 0, when the velocity HR of 
the cylinder wall exceeds the Feynman critical velocity, 
v cl = k/(2ttR) ln(i?/£). With a container radius R of 
a few mm, the Feynman velocity is only ~ 10~ 2 mm/s. 
However, if we can exclude remanent vorticity and other 
extrinsic mechanisms of vortex formation (which is gen- 
erally the case in 3 He-B at temperatures T > 0.6 T c , but 
not in 4 He-II), then vortex- free superflow will persist as 
a meta-stable state to much higher velocities because of 
the existence of a finite nucleation energy barrier. 

The height of the nucleation barrier decreases with in- 
creasing velocity. In 4 Hc-II intrinsic vortex formation 
is observed in superflow through submicron-size orifices. 
Here the barrier is ultimately at sufficiently high veloc- 
ities overcome by thermal activation or at low tempera- 
tures by quantum tunnelling (Packard 1998; Varoquaux 
et al. 1998). The other possibility is that the thermo- 
dynamic stability limit of the superfluid is reached (An- 
dreev and Melnikovsky 2003). In 3 He-B the barrier is 
unpenetrably high in practically any situation and the 
flow velocity has to be increased to the point where the 
barrier approaches zero. Vortex formation then occurs 
in the form of an instability which takes place locally at 
some sharp asperity on the cylindrical wall of the rotat- 
ing sample where the stability limit of superflow is first 
reached (Parts et al. 1995; Ruutu ct al. 1997a). The 
average velocity at the cylindrical container wall we call 
the spontaneous critical velocity v c (T, P) (Fig. 13), which 
now depends on surface roughness and is therefore con- 
tainer specific. When the instability occurs and a vortex 
is formed, the actual local velocity at the critical dom- 
inating asperity equals the intrinsic bulk instability ve- 
locity v ch (T, P). 

When the flow reaches the instability value v c b at the 
dominating sharpest protrusion, a segment of a vortex 
ring is formed there. The energy of the smallest possible 
vortex ring, which has a radius comparable to the coher- 
ence length, is of order E v <~ p s K 2 £. This energy con- 
stitutes the nucleation barrier. On dimensional grounds 
we may estimate that E v /(kBT) ~ (£/a) (Tp/T), where 
Tp = h 2 /{2mza 2 ) <~ IK is the degeneracy temperature 
of the quantum fluid and a the interatomic spacing. In 
4 He-II the coherence length is of order £ ~ a and the 
temperature T <~ Tp, which means that E v /(k-sT) ~ 1 
and thermally activated nucleation becomes possible. 
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FIG. 4. (Top) Metastable state with a central vortex cluster 
in the rotating cylinder. This axially symmetric arrangement 
consists of rectilinear vortex lines surrounded by vortex-free 
counterflow. (Bottom) Radial distribution of the velocities of 
the superfiuid component, v s , and the normal component, v n - 
Because of its large viscosity, the normal fraction corotates 
with the container. (From Krusius et al. 1994). 

In contrast in 3 He-B the coherence length is more than 
10 nm and we find E v /(k^,T) <~ 10 5 at temperatures of 
a few mK or less. The consequence from such an enor- 
mous barrier is that both thermal activation and quan- 
tum tunnelling arc out of question as nucleation mech- 
anisms. The only remaining intrinsic mechanism is a 
hydrodynamic instability, which develops at rather high 
velocity. This fact has made possible the measurement 
of neutron-induced vortex formation as a function of the 
applied flow velocity v < v c . 

2.6.2. Rotating states of the superfiuid 

When a container with superfiuid is set into rotation 
and the formation of vortex lines is inhibited by a high 
energy barrier, then the superfiuid component remains 
at rest in the laboratory frame. Its velocity is zero in 
the whole container: v s = 0. The normal component, 
in contrast, corotate with the container and v n = ilr. 
This state of vortex-free counterflow is called the Landau 
state. It corresponds to the Meissner state in supercon- 
ductors, with complete flux expulsion. 

In 4 He-II vortex-free flow with sufficiently high veloc- 
ity is out of reach in most situations, because even at 
very low velocities remanent vorticity leads to efficient 
vortex formation. In 3 He-B the superfiuid coherence 
length is several orders of magnitude larger and remanent 
vorticity can be avoided if the container walls are suffi- 
ciently smooth. Coupled with a high nucleation barrier, 
metastable rotating states then become possible. These 
can include an arbitrary number TV of rectilinear vortex 



lines, where N < iV max . In most experimental situations 
the maximum possible vortex number, N max , equals that 
in the equilibrium vortex state, N cq 2-k£IR 2 / k. The 
equilibrium vortex state is obtained by cooling the con- 
tainer slowly at constant through T c . 

If the number of lines is smaller than that in the equi- 
librium state, then the existing vortex lines are confined 
within a central vortex cluster, as shown in Fig. 4. The 
confinement comes about through the Magnus force from 
the normal-superfluid counterflow, v = v n — v s , which cir- 
culates around the central cluster in the vortex-free re- 
gion. Within the cluster mutual repulsion keeps the lines 
apart such that they form a triangular array at the equi- 
librium density, n v = 20/ 'k, which is constant over the 
transverse cross section of the cluster. At this line den- 
sity the cluster corotates with the container at constant 
fi, like a solid body. The total number of lines is given 
by N — irn v Rl, where R v is the radius of the cluster. 

With N vortex lines in the cluster, the supcrfiow veloc- 
ity outside has the value v s = kN/ (2irr) , which is equiv- 
alent to that around a giant vortex with N circulation 
quanta k (Fig. 4). The counterflow velocity v vanishes 
inside the cluster (measured on length scales which ex- 
ceed the inter- vortex distance) , while outside it increases 
from zero at r = R v to its maximum value at the cylinder 
wall, r — R: 

In the absence of the neutron source, new vortex lines 
are only formed if v(R) is increased to its container- 
dependent spontaneous critical value v c (T, P). 

2.6.3. Experimental setup 

Fig. 5 shows a measuring setup for rotating NMR mea- 
surements in 3 He-B. The cylindrical sample is connected 
via an orifice to a long liquid 3 He column which is needed 
to fill the sample volume, to provide the thermal contact 
with the nuclear refrigeration stage, and to isolate the 
sample in a metal-free environment for NMR measure- 
ment. The major resistances in the thermal path be- 
tween the sample and the nuclear cooling stage are the 
orifice, the liquid column, and the Kapitza surface resis- 
tance of the porous sintered heat exchanger in a large 
liquid volume below the column. The most important 
characteristic of the sample container itself is its spon- 
taneous critical velocity v c {T,P), which may depend in 
complicated ways on temperature and pressure, since dif- 
ferent sources contribute to vortex formation. These arise 
from (i) the surface roughness of the cylindrical walls in 
the volume above the orifice, (ii) sites in which remanent 
vortices can be trapped, and (iii) the leakage of vortices 
through the orifice (since the volume in contact with the 
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FIG. 5. Liquid 3 He sample container with coils for NMR mea- 
surement. The cylindrical sample on the top is connected via 
a narrow aperture with a long liquid He column which serves 
as a thermal path to the sintered heat exchanger on the nu- 
clear refrigeration stage. The nuclear stage is located below 
this structure while the mixing chamber of the pre-cooling di- 
lution refrigerator is above. The excitation coil and the outer 
magnet have not been drawn to scale. (From Ruutu et al. 
1997a). 

rough heat exchanger surfaces will be filled with the equi- 
librium number of vortices already at low fi). In practice 
for a clean quartz cylinder the most important source is 
some isolated surface protrusion, a localized defect or a 
piece of dirt, which increases the flow velocity at a sharp 
asperity by up to an order of magnitude (Parts et al. 
1995; Ruutu et al. 1997a). To achieve a high critical ve- 
locity, the geometry and surface quality of the sample 
container are of utmost importance. 

The sample container structure in Fig. 5 is prepared 
from fused quartz. The sample cylinder itself has a ra- 
dius R = 2.5 mm, length Z = 7 mm, and wall thick- 
ness 0.5 mm. An orifice of 0.5 mm diameter links it with 
the 60 mm long connecting column. The long length is 
needed to place the sample in the middle of the NMR 
polarization magnet. The container is fused together 
from tubular and flat plate parts in an oxygen-acetylene 
flame. After assembly the structure is annealed in an 
oven, etched, and cleaned with solvents. All indications 
point to an average surface roughness well below 1 fim, 
such that it is dirt particles and perhaps isolated localized 
defects on the glass surface which ultimately control the 
critical velocity. For the quartz glass cylinder in Fig. 5 
the spontaneous critical velocity v c was found to be a 



factor of three or more larger than the lowest velocity 
v cn (T, P) at which vortices start to appear in neutron 
irradiation (cf. Fig. 13). 

The non-invasive continuous-wave NMR measurement 
is conducted from the outside with a system of three 
orthogonal superconducting coils. The innermost is a 
saddle-shaped detector coil, wound from 25 /im solid Nb 
wire on a thin epoxy shell and mounted directly on the 
sample cylinder. It is thermally anchored with a Cu strip 
to the nuclear cooling stage. The outermost coil is a 
magnet which produces the axially oriented homogeneous 
polarization field. This end-compensated solenoid on a 
bronze body is fixed mechanically and thermally to the 
mixing chamber. The coil in the middle is used for rf 
excitation. It consists of two turns in each half and is 
fixed inside the polarization magnet. All parts shown 
here are located inside a superconducting Nb shield, to 
avoid interference from the demagnetization field which is 
required for cooling and temperature stabilization. The 
Nb jacket is part of the heat shield which is fixed to the 
mixing chamber. 

The neutron source is a weak Am/Be specimen 
which is sandwiched between two paraffin moderator tiles 
so that a cubic box results with sides of 25 cm length. For 
irradiation the box is carried to the cryostat and placed 
at room temperature at a desired distance from the 3 He 
sample. By varying this distance the neutron flux can be 
adjusted to the required low value, as determined from 
the vortex formation rate N in Fig. 8. Between measure- 
ments the source is moved a fair distance away from the 
cryostat (~ 25 m). 

2.6.4. NMR measurement 

The unusual NMR properties of the 3 He superfluids arc 
a direct consequence from the Cooper pairing in states 
with spin S = 1 and orbital momentum L = 1. In 3 He-B 
the NMR absorption spectrum, recorded at low rf exci- 
tation with traditional transverse cw methods, is related 
to the spatial distribution of the order parameter orien- 
tation in the rotating cylinder. Spin-orbit coupling, al- 
though weak, is responsible for the large frequency shifts 
which give rise to absorption peaks shifted from the Lar- 
mor frequency (Vollhardt and Wolfle 1990). The height 
of the so-called counterflow (CF) peak, which depends on 
the CF velocity around the vortex cluster, can be used 
to determine the number of rectilinear vortex lines in the 
cluster. 

In the left panel of Fig. 6 NMR absorption spectra are 
shown, which have been recorded at the same f2, but 
with a different number N of vortex lines in the vortex 
cluster. In the vortex-free state the large CF peak on the 
right increases rapidly with f2 when the orienting effect 
from the CF grows. The absorption intensity for this 
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FIG. 6. NMR measurement of vortex lines in superfluid 3 He-B, when rotated in a long cylinder with the magnetic field oriented 
axially: (Left) NMR spectra recorded at a reference velocity Q Te f = 1.2 rad/s with different number of vortex lines N . The 
Larmor frequency is at the left and the counterflow peak at the right vertical arrow. (Right) Two calibration measurements of 
the number of vortex lines N as a function of the ratio of the Larmor and counterflow peak heights, /iLar//icf , measured at two 
different rotation velocities Q ie f. (From Xu et al. 1996). 



growth in peak height is shifted from the asymmetric 
peak on the left close to the Larmor frequency. This 
second peak reflects the order parameter orientation in 
the center of the container in the region of the vortex 
cluster. At fixed fl, if more vortex lines are added, the 
central vortex cluster expands, the CF velocity outside 
the cluster is reduced, and the absorption intensity from 
the CF peak is shifted back to the Larmor peak (Kopu 
et al. 2000). 

The frequency shift of the CF peak increases monotoni- 
cally with decreasing temperature and is used to measure 
the temperature. The total integrated absorption in the 
NMR spectrum is proportional to the B phase suscepti- 
bility \b(T,P) and decreases rapidly with temperature. 
However, by maintaining the CF peak at fixed frequency 
shift, the temperature is kept stable, and the relative 
amount of absorption in the Larmor and CF peaks can 
then be used to determine the number of vortex lines. 
By monitoring the peak height of either of the two ab- 
sorption maxima, one can detect a change in N, as was 
illustrated in Fig. 3. 

In the right panel of Fig. 6 the ratio of the two peak 
heights has been calibrated to give the vortex number 
at two different angular velocities f2 re f. In both cases 
only a relatively small number of vortex lines N <C N cq 
is present and the CF peak is large (Xu et al. 1996). 
These calibration plots were measured by starting with 
an initially vortex- free sample rotating at fi re f , into which 
a given number of vortex lines N was introduced with 
neutron absorption reactions, as in Fig. 3. After that the 
irradiation was stopped and the entire NMR spectrum 
was recorded (by sweeping the magnetic polarization field 
H). From the spectrum the ratio of the peak heights of 
the two absorption maxima is worked out and plotted 
versus N. This gives the linear relationships in Fig. 6. 



The reduction in the CF peak height from the addition 
of one single vortex line can be discerned with good res- 
olution in favorable conditions (Fig. 3). The optimiza- 
tion of this measurement has been analyzed by Kopu 
et al. (2000). With a small-size superconducting mag- 
net and its modest field homogeneity (AH/H ~ 10~ 4 ), 
the best conditions are usually achieved at low magnetic 
field (H ~ 10 - 20 mT) close to T c (T > 0.8 T c ) and 
at relatively high CF velocity (£1 > 0.6 rad/s), where 
the CF peak is well developed (Korhonen et al. 1990). 
The single- vortex signal, i.e. the change in peak height 
per vortex line, decreases with increasing f2 and CF peak 
height, as shown in Fig. 7. In neutron irradiation mea- 
surements, the rate of vortex line creation N is deter- 
mined directly from records like the two traces in Fig. 3. 
The plot in Fig. 7 provides a yard stick for the single- 
vortex signal, to estimate the number of new lines if mul- 
tiple lines are created in one neutron absorption event. 

2.7. Vortex formation in neutron irradiation 

A measurement of vortex-line formation in neutron ir- 
radiated normal-superHuid counterflow is performed at 
constant ambient conditions. The externally controlled 
variables include the rotation velocity Q, temperature T, 
pressure P, magnetic field H, and neutron flux <fi n . The 
initial state is one of vortex- free counterflow (N = 0). 
When stable conditions have been reached, a weak neu- 
tron source is placed at a distance d from the 3 He-B 
sample and the output from the NMR spectrometer is 
monitored, as shown in Fig. 3. From this NMR absorp- 
tion record as a function of time, the vortex lines can 
be counted which are formed during a given irradiation 
period. 
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FIG. 7. Absolute magnitude of the single-vortex signal (as 
in Fig. 3), measured in the vortex- free counterflow state at 
different £1. The reduction in the height of the counterflow 
peak is given here in millivolts at the output of the cooled 
preamplifier which has a gain of £s 10 and is located inside 
the vacuum jacket of the cryostat. The scatter in the data 
arises from residual differences in the measuring temperature, 
perhaps from small differences in the order-parameter texture 
from one cool down to the next, and from an electronic in- 
terference signal which is generated by the rotation of the 
cryostat. More details on the NMR spectrometer are given 
by Parts et al. (1995) and Ruutu et al. (1997a). 

The process evidently exhibits stochastic variation. To 
measure the vortex- formation rate dN/dt — N, records 
with a sufficiently large number of detected neutron ab- 
sorption events are analyzed to obtain a representative 
result. After that one of the experimental parameters is 
changed and a new run is performed. In this way the 
dependence of the vortex formation process on the exter- 
nal variables can be studied. It was found that N varies 
as a function of all the external parameters, i.e. neutron 
flux, rotation velocity, temperature, pressure, and mag- 
netic field. In the following we shall first describe these 
empirically established dependences. 

The rate of vortex-line formation is proportional to 
the neutron flux. The latter is varied by changing the 
distance of the source from the cryostat. In this way 
individual absorption events can be studied, which are 
well separated in time. A calibration plot of the measured 
rate N, as a function of the distance between sample and 
source, is shown in Fig. 8. By means of this plot the 
results can be scaled to correspond to the same incident 
neutron flux, for example with the source at the minimum 
distance d = 22 cm, which is given by the outer radius of 
the liquid He dewar. 

A most informative feature is the dependence of N on 
rotation velocity (Ruutu et al. 1996a, b, 1998a). Rotation 
produces the applied bias velocity at the neutron absorp- 
tion site. Since the absorption happens close to the wall 
of the container, the velocity of this bias flow is approxi- 
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FIG. 8. Rate of vortex-line formation in neutron irradiated 
counterflow, N, as a function of the distance d between neu- 
tron source and 3 He-B sample. The initial state for all data 
points is vortex-free counterflow at a rotation velocity of 0.91 
rad/s. The neutron irradiation time is 30 min. The fitted 
curve is of the form In [1 + (R B /d) 2 ], where R B is the radius of 
the front surface of the paraffin moderator box (perpendicular 
to d), in which the Am-Be source is embedded. 

mately equal to v(R) (Eq. (7)), the counterflow velocity 
at the cylinder wall. We shall denote this velocity sim- 
ply with v. It depends both on the angular velocity ft 
and the number of vortex lines N(t) which are already 
present in a central cluster. The bias flow provides the 
force which allows vortex rings to escape from the heated 
neutron bubble and to expand to rectilinear vortex lines, 
which are then preserved in the central cluster. 

These features can be seen in Fig. 9. It demonstrates 
that it is the counterflow velocity v(N), which is the im- 
portant variable governing the rate N and not, for in- 
stance, that of the normal component, v n (R) = QR. Here 
the total vortex number N (t) is recorded during neutron 
irradiation at constant external conditions over a time 
span of six hours. The irradiation is started from the 
vortex-free state: N(t = 0) = 0. The rate N(t), at which 
vortex lines accumulate in the center of the container, is 
not constant: Initially both the counterflow velocity and 
the rate of vortex formation are the highest. When more 
vortex lines collect in the central cluster, both the veloc- 
ity v and the rate N fall off. Finally the vortex number 
N approaches a saturation value, beyond which no more 
lines are formed. This means that there exists a lower 
limit, a threshold value for the counterflow velocity v, 
below which the neutron absorption events produce no 
vortex lines at all. 

The dependence of A^ on the counterflow velocity v is 
obtained from experiments like that in Fig. 9. However, 
more efficient is a measurement of only the initial slope 
of the N(t) record. Typically 15-30 min accumulation 
periods are then enough. Only a small number of vor- 
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FIG. 9. Cumulative number of vortex lines N(t) as a function 
of time t, after turning on the neutron irradiation at constant 
flux on an initially vortex-free He-B sample rotating at con- 
stant £2. The dashed curve represents a fit to Eq. (10) with 
7 = 1.1 min^ 1 and v cn — 1.9 mm/s. The result demon- 
strates that the rate of vortex-line formation is controlled by 
the counterflow velocity v(N) and not by that of the normal 
component v n (R) = QR, which is constant during this entire 
measurement. (From Ruutu et al. 1996b). 
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FIG. 10. Initial rate N of vortex line formation at different 
rotation velocities £2. N has been measured as the average 
of all vortex lines formed during an accumulation period of 
15 - 30 min. The fitted curves are of the form 7 [(£2/£2 cn ) 3 — 
1] (Eq. (8)). Their low velocity end points determine the 
threshold velocity v cn = £2 cn .R. 
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FIG. 11. Vortex-formation rate TV plotted versus the cube 
of the bias velocity v. This plot is used to determine the 
threshold velocity v cn . A line has been fit to each set of data, 
measured under different external conditions, to identify the 
horizontal intercept v^ n and the common vertical intercept 
-7 = -1.4 min" 1 . (From Ruutu et al. 1996a). 



tices form during such a run compared to the number of 
vortices in the equilibrium state. Thus the total decrease 
in the counterflow velocity is small and can be accounted 
for by assigning the average of the velocities before and 
after irradiation as the appropriate value of the applied 
bias velocity. An example of the measured rates is shown 
in Fig. 10. 

These measurements reveal a vortex-formation rate N 
as a function of v which has an onset at a threshold ve- 
locity v cn , followed by a rapid non-linear increase. As 
shown in Fig. 3, close above the threshold a successful 
neutron absorption event produces one vortex line, but 
at high flow rates many lines may result. Therefore the 
nonlinear dependence N(SY) arises in the following man- 
ner: Initially the fraction of successful absorption events 
increases, ie. of those which produce one new vortex line. 
Eventually this effect is limited by the neutron flux, when 
essentially all absorbed neutrons produce at least one new 
vortex. The second part of the increase is brought about 
by the fact that more and more lines are produced, on 
an average, in each absorption event. 

As demonstrated in Fig. 11 the dependence of the rate 
N on the counterflow velocity v can be approximated by 
the empirical expression 



N(v) = 7 



v 

Vcn 



(8) 



with 7 and v cn as parameters. This equation also de- 
scribes the results in Fig. 9: If initially v in the vortex- 
free state is only slightly larger than the critical velocity 
v cn , then Eq. (8) can be linearized to N ~ 3j(v/v cn — 1)- 
Thus N(t) is obtained by integration from 
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FIG. 12. Neutron irradiation and vortex-line yield at high 
rotation velocity. The total number of lines, which are cre- 
ated in one irradiation session, can be determined in several 
ways: 1) In this example the direct count of the steps in 
signal amplitude yields 88 lines (cf. Fig. 3). The number 
next to each step denotes its equivalent in rectilinear vortex 
lines. 2) From Fig. 7 one finds that the drop in signal am- 
plitude per vortex line is 11 mV in the present measuring 
conditions. The total amplitude drop of 960 mV corresponds 
thus to 87 vortex lines. 3) A measurement of the annihila- 
tion threshold, by successive deceleration to lower and lower 
fl (Ruutu et al. 1997a), gives the total number of lines in 
the vortex cluster. Here this measurement yields f2 v = 0.245 
rad/s. The annihilation threshold, which at £l v > 0.2 rad/s 
coincides with the equilibriu m vortex sta te, corresponds to 
N = itR 2 2CI v k~ 1 (1 - 0.18\Arads _1 /fiv), as measured by 
Ruutu et al. (1998b). With R = 2.45 mm, this gives 88 vor- 
tex lines. 4) The rate equation (8) , with the measured values 
of 7 and v cn from Fig. 11, gives 88 lines for an irradiation 
time of 8.5 min (with fi cn = 0.75 rad/s). Thus all of the re- 
duction in the height of the counterflow peak (cf. Fig. 6) can 
be reliably accounted for in many different ways. 



j y( t ) = 37 l^-^ 2 ^-i 



2t:Rv ci 



with the solution 

N(t) = — K0) - v cn ] 



1 — exp — 



(9) 



• (10) 



This equation has been fit to the measurements in Fig. 9, 
to give the two parameters 7 and v cn . The resulting 
values agree with those of the horizontal and vertical zero 
intercepts in Fig. 11. 

At high rotation velocities a large number of vortex 
lines is produced in each neutron absorption event. The 
accumulation record for one neutron irradiation session 
(similar to the upper trace in Fig. 3) then consists of 
many large steps where, in the worst case, two absorption 



events might even be overlapping. Nevertheless, the final 
total number of lines, which have been collected into the 
central cluster, can be determined in a few different ways 
and these give consistent answers. These tests have been 
described in Fig. 12. One of them includes the use of the 
rate equation (8) where the multiplier 7 has the value 
from Fig. 11. 

In fact, it can be guessed from Fig. 11 that Eq. (8) 
has wider applicability than is apparent from the present 
examples: The same equation holds universally under dif- 
ferent externally applied conditions with all dependence 
on external variables contained in v cn (T,P,H), while 
7 is proportional to the neutron flux but does not de- 
pend on temperature, pressure, or magnetic field. Hence 
N(v/v cn ) appears to be a universal function for all mea- 
surements which have been performed in the temperature 
regime T > 0.8 T c . 

To summarize, the experimental result displays two 
distinguishing features: 1) the cubic dependence on the 
bias v and 2) the universality that all dependence of 
the vortex formation properties on the experimental vari- 
ables T, P, and H is contained in the threshold velocity 

fen- 



2.8. Volume or surface mechanism? 

In the following sections we examine the various ex- 
perimental properties of the neutron-induced vortex for- 
mation. The analysis demonstrates that most features 
are consistent with the KZ mechanism - none have so far 
been found contradictory. But does this fact constitute 
solid proof of the mechanism? 

Neutron-induced vortex formation in a rotating super- 
fluid does not exactly fit the ideal KZ model — a rapid 
second order transition in an infinite homogeneous sys- 
tem. In the heated neutron bubble there is a strong ther- 
mal gradient and a strict boundary condition applies at 
its exterior, imposed by the bulk supcrfluid state outside. 
The comparison of experiment and model is further com- 
plicated by the fact that the cool down occurs so fast that 
an extrapolation from the equilibrium state theories be- 
comes uncertain, whether it concerns the hydrodynamics 
or even the superfiuid state itself. (The applicability of 
the KZ mechanism to inhomogeneous transitions is con- 
sidered in more detail in Sec. 2.15.) Perhaps, other pro- 
cesses of hydrodynamic origin can be suggested which 
account for the experimental observations? 

The most viable alternative arises from the boundary 
condition in the presence of the applied bias flow: It is 
the instability of superflow along the normal-supcrfluid 
interface in the outer peripheral regions of the neutron 
bubble. In the Ginzburg-Landau temperature regime 
the container-dependent spontaneous critical velocity 
v c (T, P) decreases with increasing temperature and van- 
ishes at T c (dashed curves in Fig. 13). It is related (via 
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surface roughness enhancement) to the intrinsic instabil- 
ity velocity v c b(T,P) of the bulk superfluid (Parts et al. 
1995; Ruutu et al. 1997a), which has qualitatively a sim- 
ilar temperature and pressure dependence, but is larger 
in magnitude. In the outermost region of the neutron 
bubble the fluid remains in the B phase, but is heated 
above the surrounding bulk temperature To. Conse- 
quently, if no other process intervenes, the superflow in- 
stability has to occur within a peripheral shell surround- 
ing the hot neutron bubble where v c \,(T, P) drops below 
the applied bias v. This process was discussed by Ruutu 
et al. (1998a) and compared to rotating measurements, 
but was found incompatible with the measured results. 
As explanation it was suggested that the KZ mechanism 
is inherently the fastest process for creating vortices. 

The interplay between the interior volume of the neu- 
tron bubble and its boundary was examined by Aran- 
son et al. (1999) in a rapidly cooling model system with 
a scalar order parameter, using the thermal diffusion 
equation to account for cooling and the time-dependent 
Ginzburg-Landau equation for order-parameter relax- 
ation. Rigorous analytic calculations showed that the 
normal-superfluid interface becomes unstable in the pres- 
ence of superflow along the interface. Numerical calcula- 
tions (Aranson et al. 1999, 2001) further confirmed that 
the instability develops such that vortex rings are formed. 
The rings encircle the bubble perpendicular to the ap- 
plied flow and screen the superflow, so that the super- 
fluid velocity (in the rotating frame) is zero inside the 
bubble. The number of rings depends on the shape and 
size of the bubble and on the bias velocity. In this pro- 
cess it does not matter what the state of the liquid is 
inside the bubble: The rings are formed within an exter- 
nal cooler layer, which remains throughout the process in 
the B phase. After their formation, the rings start to ex- 
pand and might eventually be pulled away by the Magnus 
force. Thus these rings could be the source for the vortex 
lines which are observed in the rotating measurements. 

However, the calculations also confirmed that in addi- 
tion to the superflow instability at the bubble boundary, 
the KZ mechanism in the bubble interior is also present 
- thus in this calculation both processes produce vor- 
tices. In the simulation the loops, which eventually man- 
age to escape into the bulk bias flow, originate from the 
boundary while the random vortex network in the inte- 
rior collapses and decays away. Thus Aranson et al. con- 
firmed the presence of the KZ mechanism, but concluded 
that the vortices observed in the rotating measurements 
are produced by the boundary instability. The time- 
dependent Ginzburg-Landau treatment with one com- 
mon relaxation time for both the phase and the mag- 
nitude of the order parameter is rigorously obeyed only 
within a narrow temperature interval from T c down to 
1 — T/T c ^ 10 ~ 3 . The bulk liquid temperature in the ro- 
tating measurements is at least an order of magnitude be- 



low this range (see Fig. 13). Another difference between 
experiment and calculation lies in the size of the neutron 
bubble: The numerical calculation is carried out on a 
grid with unit length smaller than the coherence length 
£(Tb). To make this tractable, the reaction energy Eq is 
reduced by one or two orders of magnitude. This scales 
down the neutron bubble radius i?b via the thermal diffu- 
sion equation, which is used to model the neutron bubble 
(see Sec. 2.9.1). A smaller bubble (see Eq. (14)) means 
a faster velocity of the phase front (oc 1/Rb). At present 
time it is not known how reliably the calculation with its 
differences extrapolates into the range of the measure- 
ments. 

How do we, in fact, distinguish from measurements 
whether it is the KZ mechanism or the boundary insta- 
bility which is responsible for the observed vortex lines? 
Clearly the critical velocity v cn does not discriminate be- 
tween them. In both cases v cn is the velocity of the bias 
flow at which the largest vortex ring which fits in the 
bubble is stable, as will be shown in Sec. 2.9. However, 
the different nature of the two processes has important 
consequences for the vortex formation rate N(v). The 
row of vortex rings, which is produced by the boundary 
instability to screen the interior of the neutron bubble 
from the bias flow, resembles a vortex sheet. The sheet 
is formed in the heated peripheral shell of the neutron 
bubble, where the outside superfluid moves with velocity 
v along the bubble interface and the superfluid inside is 
stationary. The density of vorticity in the sheet is v/k 
and the number of loops produced by one neutron ab- 
sorption event is N oc vR^/k, where i?b is the size of 
the bubble along the flow direction. Thus N(v) grows 
linearly with the applied bias flow velocity, N(v) oc v, 
which agrees with the result of Aranson et al. (1999). In 
contrast, the KZ mechanism is a volume effect, which re- 
sults in the cubic dependence expressed in Eq. (8) (see 
also Sec. 2.12.2). The measured result in the figures of 
the previous section is not linear, but is consistent with 
the cubic dependence. 

Are there other arguments emerging from the mea- 
surements which would support the KZ origin of the ob- 
served vortices? Several more can be listed in addition to 
the rate dependence N(v), which favor a volume effect: 
(i) The boundary instability is expected to be a deter- 
ministic process which in every neutron event produces 
the same surface density of vortex rings. Variations in 
the number of rings arise only owing to variations in the 
shape of the neutron bubble and its orientation with re- 
spect to the bias flow. (Unfortunately, the distribution 
of the number of rings was not studied in detail in the 
simulations by Aranson et al. 2001.) In contrast, the 
KZ mechanism produces a random vortex network from 
which the number of loops extracted by the bias flow 
varies and follows a stochastic distribution (Fig. 30) . (ii) 
The boundary instability should not be particularly sen- 
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FIG. 13. Critical values of the applied bias velocity v vs. 
temperature at different pressures: (Solid lines) Threshold 
velocity v cn oc (1 — T/Tc) 1 / 3 for vortex formation in neutron 
irradiation (Ruutu et al. 1996a). (Dashed lines) Critical ve- 
locity v c oc (1 — T/Tc) 1 / 4 of spontaneous vortex formation in 
the same quartz glass container in the absence of the neutron 
source. (From Parts et al. 1995; Ruutu et al. 1997a). 



sitive to the processes inside the heated neutron bubble 
and only regular mass-flow vortices are expected to form 
in such an instability. The KZ mechanism, on the other 
hand, can be expected to produce all possible defects. 
Their presence might be either directly observed in the 
final state or via their influence on the evolution of the 
vortex network inside the neutron bubble. Such observa- 
tions are discussed in Sec. 2.10. 



2.9. Threshold velocity for vortex loop escape 

The dependence on the counterflow bias v can be stud- 
ied from the threshold v cn up to the critical limit v c at 
which a vortex is spontaneously nucleated at the cylin- 
drical wall in the absence of the neutron flux (Parts 
et al. 1995; Ruutu et al. 1997a). The threshold velocity 
v cn (T, P 1 H) is one of the features which can be examined 
to learn more about neutron-induced vortex formation. 
Measurements of both critical velocities, v cn and v c , are 
shown in Fig. 13 as a function of temperature. 



2.9.1. Properties of threshold velocity 

By definition the threshold v cn represents the smallest 
bias velocity at which a vortex ring can escape from the 
neutron bubble after the absorption event. It can be 
connected with the bubble size in the following manner. 
A vortex ring of radius r Q is in equilibrium in the applied 
bias flow at v if it satisfies the equation 



As explained in more detail in Sec. 2.12.2, a ring with 
a radius larger than r a will expand in the flow while a 
smaller one will contract. Thus the threshold or mini- 
mum velocity at which a vortex ring can start to expand 
towards a rectilinear vortex line corresponds to the maxi- 
mum possible vortex-ring size. This must be comparable 
to the diameter of the heated bubble. For a simple es- 
timate we set the vortex ring radius equal to that of a 
spherical neutron bubble: r (u C n) ~ Rb- (In fact, the 
numerical simulations to be described in Sec. 2.14.3 sug- 
gest that r (v cn ) ~ 2i?b because of the complex convo- 
luted shape of the largest rings in the random vortex 
network.) 

A simple thermal diffusion model can be used to yield 
an order of magnitude estimate for the radius i?b of the 
bubble which originally was heated above T c . In the tem- 
perature range close to T c the cooling occurs via diffusion 
of quasiparticle excitations out into the surrounding su- 
pcrfluid with a diffusion constant D w vpl, where vp is 
their Fermi velocity and I their mean free path. The dif- 
ference from the surrounding bulk temperature T as a 
function of the radial distance r from the centre of the 
bubble can be calculated from the diffusion equation 



dT(r,t) _ (d 2 T 2dT 
dt \ dr 2 r dr 



(12) 



With the assumption that at t = the reaction energy 
E is deposited at r = 0, the solution is given by 



T(r,t)-T «^ 



1 




77TT ' (13) 



C v (AirDt) 3 / 2 



where C' v is the specific heat. For now, we may assume 
that all the energy of the neutron absorption reaction is 
deposited as heat. The bubble of normal fluid, T(r) > T c , 
first expands and reaches a maximum radius 



C V T C 



1/3 



(1 - T /T c ) 



-1/3 



(14) 



(11) 



It then starts cooling and rapidly shrinks with the char- 
acteristic time tq <~ R±,/D <~ 1/xs. Since v cn is inversely 
proportional to r <~ i?b, it has the temperature depen- 
dence v cn oc (1 — Xb/Tc) 1 / 3 . This is in agreement with 
the solid curves in Fig. 13 which have been fitted to mea- 
surements on v cn . The prefactor of these curves is in 
agreement with that from Eqs. (14) and (11) within a 
factor of ~ 2, and its increase with increasing pressure 
is well described by the decrease in bubble size accord- 
ing to Eq. (14), where C v and T c increase with pressure 
(Greywall 1986). 



2.9.2. Influence of 3 He-A on the threshold velocity v cn 

Measurements on the threshold velocity v cn offer inter- 
esting possibilities to study whether the interior of the 
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FIG. 14. Threshold velocity v cn for the onset of vortex for- 
mation during neutron irradiation: (Left) Phase diagram of 
3 He superfluids in the pressure vs. temperature plane, with 
the A^B transition at Tab(0) in zero field (solid line) and 
at Tab(H) in nonzero field (dashed line). Two quench tra- 
jectories, distinguishing different types of measurements on 
the right, are marked with (a) and (b). (Right top) The pres- 
sure dependence of v cn displays a steep change at the pressure 
Ppcp of the polycritical point. (Right bottom) The magnetic 
field dependence of v cn is parabolic, similar to that of the equi- 
librium state A^B transition Tab(H). (From Ruutu et al. 
1998a). 



heated bubble participates in vortex formation or not. 
In the left panel of Fig. 14 the phase diagram of 3 He is 
shown at temperatures below 2.5 mK. In the measure- 
ments which we discussed so far, the liquid pressure has 
been below 21.2 bar, so that the heated bubble cools from 
the normal phase directly into the B phase (along trajec- 
tory b). On the top right in Fig. 14 the measured v cn is 
plotted as a function of pressure P at constant reduced 
temperature T/T c . Here the pressure dependence dis- 
plays an abrupt increase at about 21.2 bar, the pressure 
Ppcp of the polycritical point: It thus makes an unex- 
pectedly large difference whether the quench trajectory 
follows a path denoted with a or with b! The two cool- 
ing trajectories differ in that above Ppcp a new phase, 
3 He-A, is a stable intermediate phase between the nor- 
mal and B phases. Consequently, although the bulk liq- 
uid is well in the B phase in all of the measurements of 
Fig. 14, vortex formation is less than expected on the 
basis of extrapolations from lower pressures, when the 
quench trajectory crosses the stable A-phase regime. 

The fitted curves in Fig. 14 represent v cn (P) at pres- 
sures below Ppcp 



v cn = AK/(iTTR h ) In (fib/0 , 



(15) 



curves has the value A — 2.1. This is the only fitting pa- 
rameter in Figs. 13 and 14 where the same fit is compared 
to measurements as a function of both temperature and 
pressure. The agreement is reasonable and suggests that 
the spherical thermal diffusion model is not too far off. 
However, as seen in Fig. 14, v cn is clearly above the fit 
at P > Ppcp- Also this offset is largest at the highest 
temperature, which in Fig. 14 is 0.96 T c . Both features 
suggest that vortex formation is reduced when the rela- 
tive range of A-phase stability increases over the quench 
trajectory below T c . 

The lower right panel of Fig. 14 shows the dependence 
of v cn on the applied magnetic field H. Below Ppcp the 
magnetic field acts to stabilize 3 He-A in a narrow in- 
terval from T c down to the first order A-^B transition 
at Tab(P, H). This result confirms our previous con- 
clusion: v cn again increases when the range of A-phase 
stability increases. The parabolic magnetic field depen- 
dence of v cn (H) in the lower right panel is reminiscent of 
that of the equilibrium Tab(-£0 transition temperature. 2 
It should be noted that the measurements in the upper 
right panel of Fig. 14 have not been carried out in zero 
magnetic field but in 11.7mT. However, as shown in the 
lower panel, the field dependence in the range 0-12mT 
is not visible within the experimental precision. 

Thus by increasing the pressure above Ppcp or by in- 
creasing the magnetic field we reduce the A-phase energy 
minimum relative to that of the B phase and both oper- 
ations act to increase v cn . If we go back to Figs. 10 and 
11, we note that both operations seem to leave the rate 
equation (8) and its rate constant 7 unchanged. To be 
more careful, we conclude that in first order the changes 
from increased A-phase stability appear to affect only 
v cn . This conclusion is compatible with the KZ model, 
as will be discussed in the next sections. What about the 
competing model, the superflow instability at the bound- 
ary of the heated neutron bubble, can it also account for 
the observations in Fig. 14? At first glance this does not 
appear to be the case. 

The superflow instability occurs in B phase in the outer 
periphery of the heated bubble in a shell where the tem- 
perature is close to T c . In principle this phenomenon 
should not be influenced by the appearance of A phase in 
the interior of the bubble during its cool down to the bath 
temperature To. However, it can be affected by the pres- 
ence of an A phase shell around the hot bubble, bounded 
by an inner surface T = T C and an outer T = Tab, i-e- a 
region where the temperature exceeds the AB transition 
temperature Tab(P, H). The B-phase bulk critical veloc- 



where Rb is obtained from Eq. (14). If we assume that 
all of the reaction energy is transformed to heat (Eq = 
764 keV) , then the common scaling factor A of the three 



2 According to Tang et al. (1991), below Ppcp the equilib- 
rium A^B transition temperature is in first order of the form 



T AB (P,H) 
(mT)" 2 . 



T C (P) (1 - bH 2 ), where b(P) ~ (0.5 - 10) ■ 10" 
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ity is approximately equal to the pair-breaking velocity 
and in the Ginzburg-Landau temperature regime it is of 
the form (Vollhardt et al. 1980; Kleinert 1980) 

v ch {T,P) « v c0 (P) v 7 ! - T/T c , 

where v c0 (P) = 1.61 (l+Ff/3) k B T c /p F . This velocity is 
a smooth function of pressure and cannot account for the 
steep rise of v cn at Ppcp in Fig. 14 (upper right panel). 
The B-phase boundary instability occurs within an outer 
shell where the inner surface is defined by T = T c and 
the outer by the condition v c b(T,P) = v. For the data 
in Fig. 14 the latter condition puts the outer surface of 
this shell at 1 - T/T c ~ 10~ 3 - 10" 5 , which falls inside 
the shell within which A-phase is stable. Thus here it 
is the A-phase shell within which the boundary instabil- 
ity has to occur. Such a situation is complicated (see 
Blaauwgeers et al. 2002; Parts et al. 1993) and possibly 
one where the escape of vorticity from the A-phase shell 
into B-phase is reduced. In first order these considera- 
tions suggest that v CD increases with increasing pressure 
and magnetic field, when the range of the stable A-phase 
regime increases over the length of the quench trajec- 
tory as a function of temperature. However, it is unclear 
whether in a 3-phase problem — where an outer region 
of warm supcrfluid surrounds a hot interior with normal 
liquid, while both are embedded in a cold bath of B- 
phase — coherence is established such that we can apply 
shell structure considerations. Therefore, in the absence 
of more quantitative estimates, we must conclude that 
the measurements in Fig. 14 do not yet allow a prefer- 
ence between the KZ model or the superflow instability 
at the neutron bubble boundary. Clearly more exten- 
sive measurements are needed to study rapid thermal 
quenches in the 3-phase situation, as will be emphasized 
in Sec. 2.13.2. 

2.10. Other defect structures formed in neutron irradiation 

2.10.1. Radiation-induced supercooled A — > B transition 

The first order transition of supercooled 3 He-A to 
3 He-B is another experimentally confirmed phenomenon 
which is catalyzed by the interior of the neutron bub- 
ble. Usual arguments about homogeneous first-order 
phase transitions show that the A^B transition is for- 
bidden (Leggett 1984). The normal phase symmetry 
SO(3) x SO (3) x [7(1) of liquid 3 He can be broken to the 
[7(1) x U(l) symmetry of the A phase or to the SO (3) 
symmetry of the B-phase, with only a small energy dif- 
ference between these two states, but separated by an 
extremely high energy barrier (~ 10 6 fcBT c ). Neverthe- 
less, so far it has not been reported that the transition 
would not have taken place ultimately, if only one cools 
to sufficiently low temperature. 



The A-^B transition has not been found to be affected 
by rotation (Hakonen et al. 1985) and in the absence of 
ionizing radiation it is generally believed to be caused by 
extrinsic sources, perhaps at solid surfaces. A dramatic 
increase in the transition probability has been shown to 
take place in the bulk liquid in the presence of 7 radi- 
ation or thermal neutrons (Schiffer et al. 1992; Schiffcr 
and Osheroff 1995). A generally accepted understand- 
ing of this feature is still missing, because two quite dif- 
ferent models, the "baked Alaska" (Leggett 1992; Schif- 
fer et al. 1995; Leggett 2002) and the KZ mechanism 
(Volovik 1996; Bunkov and Timofecvskaya 1998a, b), are 
discussed in this context. Since only the latter is appro- 
priate in radiation-induced vortex formation, we briefly 
describe it in the context of the supercooled A— >B tran- 
sition. 

When 3 He-A is supercooled to low temperatures, below 
the thermodynamic A^B transition, the sample stays in 
the A phase for long times, perhaps even indefinitely 
depending on temperature and other largely unknown 
requirements. The deeper the supercooling, the faster 
the transition into the B phase follows. With increas- 
ing magnetic field the transition probability decreases, 
since a magnetic field acts to stabilize the A phase. How- 
ever, when neutron irradiation is turned on, the transi- 
tion probability increases sharply. In such measurements 
the initial state is supercooled 3 He-A at varying tem- 
peratures and magnetic fields. Consequently here the A- 
phase state is enforced by the boundary condition around 
the hot neutron bubble: The boundary condition is not 
the source for the A— >B transition! Rather the new B- 
phase state has to emerge from the interior of the bubble. 

According to the KZ model, in different parts of the 
rapidly cooling neutron bubble the order parameter may 
fall initially either in the A or B phase energy minima. 
Thus a domain structure of size £ v (Eq. 3) is laid down, 
where A and B-phase blobs form a patch work. When 
these blobs later grow together, then AB interfaces are 
formed. The lower the temperature, the lower is the B- 
phase energy minimum relative to that of the A phase, 
the larger is the proportion of B-phase blobs in the patch 
work, and the more likely it is that many of them manage 
to merge together to form one large bubble, where the 
AB interface exceeds the critical diameter of about 1 /im, 
needed for spontaneous A-phase expansion to start. This 
one seed is then sufficient to initiate the A— >B transition 
in the whole sample. 

The time dependent Ginzburg-Landau calculations of 
Bunkov and Timofeevskaya (1998a, b) are consistent with 
this model. In the calculations there is no spatial depen- 
dence, an initial fluctuation is imposed at T c into a ran- 
dom direction of the phase space, and the evolution to- 
wards the final state is followed. As seen in in Fig. 15, the 
energy difference between the A and B phases is small, 
but pressure dependent. During rapid cool down through 
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FIG. 15. Normalized energy difference, AE(b-a) / E(P), be- 
tween A and B phases (right vertical axis) and probabil- 
ity of A-phase formation in a rapid cool-down (left verti- 
cal axis), plotted as a function of pressure at temperatures 
close to T c in the Ginzburg-Landau regime. The results 
have been calculated with the weak coupling (3 parameters 
{(3i — —1, 2, 2, 2, —2), except for the data point marked with 
a cross (x) on the upper curve, for which strong-coupling 
parameter values have been used according to the spin fluc- 
tuation model. (From Bunkov and Timofeevskaya 1998a, b). 



T c the order parameter may settle into either of these 
two energy minima with a pressure-dependent probabil- 
ity, which is also shown in Fig. 15. 

Thus the KZ mechanism provides a simple explanation 
of the A— >B transition probability as a function tempera- 
ture and magnetic field in supercooled 3 He-A. This model 
is independent of the boundary condition, which dictates 
3 He-A. Furthermore, in ambient conditions, where only 
B phase is stable, A-phase blobs ultimately shrink away 
and only B phase remains. Before that, however, a net- 
work of AB interfaces is created which interact with the 
formation of other topological defects. This feature can 
be used to explain the sharp increase in v cn at Ppcp in 
Fig. 14 (top right panel). 



2.10.2. Vortex formation, AB interfaces, and KZ mechanism 

The KZ model predicts that in different parts of the 
cooling neutron bubble the order parameter may fall ini- 
tially either in A or B phase energy minima, and a patch 
work of AB interfaces is laid down. The relative num- 
ber of A and B-phase blobs of size £ v is not only deter- 
mined by the difference in the A and B-phase energies: 
It also depends on the trajectory from the normal phase 
to the new energy minimum in the phase space spanned 
by the order parameter components. The latter aspect 
gives more preference to the higher energy A state: Al- 
though the A-phase is energetically not favorable in a 
wide pressure range, the probability of its formation has 



a value close to 50 % even at lower pressures, as shown 
in Fig. 15. The AB interfaces interfere with the simul- 
taneous formation of a random vortex network, which is 
laid down separately in both the A and B phases. It 
is known from experiments on both a moving AB inter- 
face (Krusius et al. 1994) and a stationary AB interface 
(Blaauwgeers et al. 2002) that the penetration of vortex 
lines through the AB phase boundary is suppressed. We 
thus expect that the presence of A-phase blobs reduces 
the combined volume which remains available for B phase 
blobs and their vortex network inside a neutron bubble. 
Consequently, when B-phase blobs merge to larger units, 
the overall volume, into which the B-phase vortex net- 
work is confined, is smaller, and B-phase vortex forma- 
tion becomes impeded. This process increases the value 
of v cn . 

To conclude, the KZ mechanism provides a unified 
view on the problems of pressure and magnetic field de- 
pendences of v cn (Fig. 14) and also the A^B transition in 
supercooled 3 He-A. In both cases AB interfaces appear 
as a new type of defect in the rapidly cooling neutron 
bubble. Normally, to create a sizeable bubble of B phase 
within bulk A liquid requires that an AB interface of 
large size is formed. This is a slow and energy consuming 
process. In quench cooling through T c the freeze-out of 
disconnected blobs with A and B-phase like fluctuations 
happens first and these are initially separated by super- 
cooled normal liquid. The AB interfaces appear later as 
metastable defects when the blobs grow together. Such 
an interpretation suggests that the KZ mechanism is the 
fastest process by which defects are created, before other 
effects manage to switch on. 



2.10.3. Spin-mass vortex 

Unlike vortex lines, which are produced in neutron ir- 
radiation and are then partly stabilized by the applied 
flow, no bias is provided for maintaining AB interfaces in 
Fig. 14. Thus in these rotating experiments, the evidence 
for AB interfaces remains indirect. However, there exists 
another topologically stable defect, which is formed in 
a neutron absorption event and is preserved by the bias 
flow for later examination. This is a combined object 
called the spin-mass vortex (Kondo et al. 1992; Korho- 
nen et al. 1993). 

The signature from the spin-mass vortex is explained in 
Fig. 16, where the counterflow peak height is plotted as a 
function of time during neutron irradiation (as in Figs. 3 
and 12). This accumulation record shows one oversize 
step in the absorption amplitude, which coincides with a 
similar discontinuous jump in the out-of-phase dispersion 
signal. As explained in the context of Fig. 12, the total 
number of vortex lines extracted from an irradiation ses- 
sion can be determined in several different ways. In the 
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FIG. 16. Signal from spin-mass vortex in neutron irradiation. 
This irradiation session includes one large jump in signal am- 
plitude. Like in Fig. 12, the various methods for extracting 
the total number of lines give: 1) A direct count of the steps 
in signal amplitude yields 45 lines plus the large jump. The 
combined amplitude drop for these 45 lines amounts to 716 
mV which translates to 16 mV/vortex, in agreement with 
Fig. 7. With this step size the large jump of 390 mV would 
correspond to 25 additional lines. 2) A measurement of the 
annihilation threshold gives £l v — 0.155 rad/s, which cor- 
responds to 48 lines. 3) The rate Eq. (8) gives N = 6.2 
lines/min or a total of 48 lines for an irradiation of 7.5 min 
duration with the measured threshold value of Q cn = 0.75 
rad/s. For these three estimates to be consistent, the large 
jump cannot be 25 lines. Instead it is interpreted to represent 
one spin-mass vortex line. In this case the size of the jump 
is not related to the number of lines, but to the counterflow 
velocity v(R) (Eq. (7)). It controls the length v(R)/2Q of the 
soliton sheet (Fig. 18). The dispersion signal also displays 
a simultaneous large jump, while it behaves smoothly when 
usual mass-current vortex lines are formed. Outside the dis- 
continuity the dispersion signal is slowly drifting, because it is 
very sensitive to a residual temperature creep. (From Eltsov 
et al. 2000). 



case of Fig. 16 a comparison of the different line counts 
shows that, within the uncertainty limits, the large jump 
cannot represent more than a few vortex lines. As will be 
seen below, this is the signature from a spin-mass vortex, 
which in Fig. 16 should be counted as one line. 

The identification of the spin-mass vortex is based on 
its peculiar NMR spectrum (Kopu and Krusius 2001). 
In an axially oriented magnetic field the NMR spectrum 
has a textural cut off frequency, beyond which no ab- 
sorption is allowed (excepting line broadening effects). 
The distinguishing feature of the spin-mass vortex is the 
absorption which has been pushed to the maximum pos- 
sible frequency shift, beyond that of the counterflow (CF) 
peak (which corresponds to 80 % of the maximum value, 
see Fig. 17). 



Because of this shifted absorption, the CF peak is 
slightly shifted both in frequency and in height in the 
presence of a spin-mass vortex, compared to its loca- 
tion without the spin-mass vortex. This difference is 
displayed in Fig. 17, where the spectrum plotted with 
the solid curve includes one spin-mass vortex while the 
dashed spectrum is without it. The small frequency shift 
between these two peaks is the reason why also the NMR 
dispersion signal is discontinuous when the spin-mass 
vortex is formed during the neutron irradiation in Fig. 16. 
Since both the NMR absorption and dispersion signals 
are continuously monitored at the lock-in-amplifier out- 
put, a jump in the dispersion channel provides the im- 
mediate alert for a spin-mass vortex. In contrast, no 
discontinuity in the dispersion signal is observed in first 
order when a small number of vortex lines is added to the 
cluster at high bias velocity (see spectra in Fig. 6). In 
this case some absorption intensity is shifted towards the 
Larmor edge to small frequency shifts, and the change in 
the dispersion signal at the site of the CF maximum is 
negligible. 

The large-scale configuration, in which the spin-mass 
vortex appears, is shown in the top part of Fig. 18. It 
consists of a linear object and a planar domain-wall-like 
soliton sheet. The line defect provides one termination 
for the sheet while its second end is anchored on the 
cylindrical wall of the container. The line defect consists 
of a singular core, which supports a circulating supcrfluid 
mass current, but has additionally trapped a disclination 
line in spin-orbit coupling. The latter feature can also 
be described in terms of a spin vortex which supports a 
persistent spin current. Thus the composite object, the 
spin-mass vortex, has the combined properties of a usual 
mass-current vortex and a spin-current vortex. A gain 
in the core energy of the combined structure provides 
the energy barrier against their dissociation. The soliton 
sheet is the continuation into the bulk supcrfluid from 
the breaking of the spin-orbit coupling: It separates two 
degenerate domains with anti-parallel orientations of the 
B-phase anisotropy axis n (lower part of Fig. 18). The 
spin-mass vortex feels the Magnus force from the applied 
bias flow, but this is partly compensated by the surface 
tension of the soliton sheet. The equilibrium position 
of the spin-mass vortex is therefore at the edge of the 
vortex cluster (top part of Fig. 18) from where it can be 
selectively removed by annihilation, as was done in the 
experiment of Fig. 17. 

The origin of the two defects in the spin-mass vortex 
can be seen from the form which the B-phase order pa- 
rameter takes (Vollhardt and Wolfle 1990): 



A B (T,P) R aj (n,6) . 



(16) 



It includes the isotropic energy gap Ab (T, P) , the phase 
factor e 1 ^, and the rotation matrix R a j(h,0). The lat- 
ter defines the rotation by which the 50(3)l orbital 
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FIG. 17. NMR signal characteristics of the spin-mass vor- 
tex. (Insert) This 20 min irradiation session includes one 
large jump, in addition to vortex-line formation. The lat- 
ter amounts to 64 lines, as determined from the measured 
annihilation threshold Q v = 0.19 rad/s. The rate Eq. (8) 
gives N = 3.35 lines/min or a total of 67 lines (with the 
measured Q cn — 1.10 rad/s). The combined drop in signal 
amplitude from these lines is 206 mV, which corresponds to 
3.2 mV/line. Here at a lower temperature the step size per 
vortex is much reduced from that in Fig. 12. The single large 
jump corresponds to 36 additional lines, but like in Fig. 16 it 
is interpreted to represent one spin- mass vortex. (Main panel) 
NMR spectra of the accumulated vortex sample, recorded at a 
lower reference velocity Q re f : The solid curve with the smaller 
counterflow peak is the original state after irradiation and de- 
celeration from 2.29 to 1.61 rad/s. The dashed curve with 
the higher counterflow peak was traced later after cycling Q 
from 1.61 to 0.20 rad/s and back, ie. after decelerating close 
to, but still above the annihilation threshold of the regular 
vortex cluster. In the former spectrum (solid curve) the soli- 
ton sheet shifts a sizeable fraction of the absorption to the 
maximum possible frequency shift (at right vertical arrow). 
After cycling fi the outermost spin-mass vortex (see Fig. 18 
top panel) is annihilated (dashed curve) and the CF peak is 
shifted to the value (at left vertical arrow) without the soli- 
ton sheet. This shift in the frequency of the counterflow peak 
creates the jump in the dispersion signal in Fig. 16. Since the 
remaining regular cluster was not decelerated at any stage be- 
low its annihilation threshold of 0.19 rad/s, no usual vortex 
lines were annihilated from the dashed-curve spectrum. 

and SO(3)s spin spaces are rotated with respect to each 
other around the axis n by the angle 9. To minimize 
the weak spin-orbit interaction, it is required that 9 = 
arccos(— i) w 104° homogeneously everywhere within 
the bulk supcrfluid. The mass-current vortex is asso- 
ciated with a 2ir circulation in the phase factor while 
the spin current involves a disclination line in R a j (n, 9) 
field. The minimum energy configuration becomes then 
the structure depicted in the lower part of Fig. 18. Here 




FIG. 18. (Top) The spin-mass vortex (SMV) is a composite 
object, which carries a trapped spin current as well as a mass 
current. In the rotating container it is confined by the Mag- 
nus force to the edge of the vortex cluster, which also includes 
usual mass-current vortex lines (MV). The extension of the 
spin vortex into the bulk liquid is a soliton tail (grey strip). 
The spin-mass vortex itself provides one termination line of 
the soliton while its second end is anchored on the cylindrical 
wall. (Bottom) Inside the soliton the angle 6 changes through 
7T while outside it is homogeneously at 104°. At the soli- 
ton surface the rotation axis n is oriented perpendicular to 
the surface while further away it goes smoothly over into the 
flare-out texture of an ideal cylinder in an axially oriented 
magnetic field. The distortion of the n texture from the flare- 
out configuration in the region around the soliton (Kopu and 
Krusius 2001) gives rise to the NMR signature of the spin- 
mass vortex in Fig. 17. 

the volume in the bulk, where the spin-orbit interaction is 
not minimized and 9 traverses through it, is concentrated 
within a soliton tail. 

In 3 Hc-B weak anisotropy energies arise in an exter- 
nal magnetic field. They produce an extended orien- 
tational distribution, or texture, of the anisotropy axis 
n. The characteristic length scale of this texture in 
the conditions of the rotating NMR measurements is 
£h oc 1/H <~ 1 mm. Since the NMR frequency is con- 
trolled locally by the orientation of n with respect to 
H, the NMR spectrum is an image of the n texture in 
the cylindrical container. In an axially oriented magnetic 
field and an infinitely long cylinder, the texture is known 
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to be of the "flare-out" form (Vollhardt and Wolfle 1990). 
In this texture there are no regions where n would lie in 
the transverse plane perpendicular to the field H. In 
contrast, at the soliton sheet n is oriented strictly per- 
pendicular to the sheet and thus to H. This part of the 
texture is responsible for shifting NMR absorption into 
the region of maximum possible frequency shift and be- 
comes thus the signature of the spin-mass vortex in the 
NMR spectrum (Fig. 17). 

The presence of the soliton sheet has many important 
experimental implications: It determines the location of 
the spin-mass vortex at the edge of the cluster and it 
leaves a clear signature in the NMR spectrum at high 
bias velocities when the sheet becomes stretched. It also 
affects the threshold velocity at which a spin-mass vortex 
might be expected to escape from the neutron bubble. 
This was only observed to happen well above v cn , the 
threshold velocity for usual mass- vortex rings (Fig. 13). 
It is explained by the fact that for spontaneous expansion 
a spin-mass vortex loop has to exceed, in addition to the 
usual barrier formed by the line tension of the small loop, 
also the surface tension from the soliton sheet. 

The spin-mass vortex was originally discovered as a de- 
fect which is formed when an A^B transition front moves 
slowly close to adiabatic conditions through a rotating 
container (Kondo et al. 1992; Korhonen et al. 1993). In 
this experiment the initial A-phase state is one with the 
equilibrium number of doubly quantized singularity-free 
vortex lines while the final state is found to contain less 
than the equilibrium number of singly quantized B-phasc 
vortex lines plus some number of spin-mass vortices. This 
means that A-phase vortex lines interact with the mov- 
ing AB interface, they are not easily converted to B-phasc 
vorticity, a critical value of bias flow velocity has to be ex- 
ceeded before the conversion becomes possible, and even 
then some fraction of the conversion leads to vorticity 
with the additional defect in R a j(h, 9) (Parts et al. 1993; 
Krusius et al. 1994). In neutron measurements the spin- 
mass vortex was not observed at 2.0 bar, but at 18.0 bar. 
In the KZ scenario this is the pressure regime where the 
probability for the formation of A phase blobs in the neu- 
tron bubble increases (Fig. 15). Interestingly one might 
therefore surmise that even in the neutron bubble the 
spin-mass vortex could result from vorticity which orig- 
inally was contained within A-phase blobs and which is 
partly transferred to the B phase when the A phase blobs 
shrink away. 

The two NMR spectra in Fig. 17 illustrate the easy 
identification of the spin-mass vortex after neutron irra- 
diation. Although the spin-mass vortex is a rare event in 
neutron irradiation, compared to the yield of vortex lines, 
its presence demonstrates that in addition to usual mass- 
current vortices also other order-parameter defects are 
created. This lends support to the discussion in Sec. 2.9.2 
that also the AB interface could be among such defects. 



The presence of the spin-mass vortex also limits the pos- 
sible mechanisms of defect formation. An explanation 
in terms of the superflow instability at the neutron bub- 
ble boundary (Sec. 2.8) is not a probable alternative: The 
spin-mass vortex is not normally created at the bulk crit- 
ical velocity limit v c b(T, P) as a response to superflow. 

2.11. Vortex formation in gamma radiation 

To study the dependence of the vortex formation rate 
N on the dimensions of the heated neutron bubble in 
the 3 He-B bath, it would be useful to investigate alterna- 
tive heating methods. In addition to thermal neutrons, 
beta and gamma radiation have frequently been used as 
a source of heat. In the rotating experiments a 60 Co 
gamma source was experimented with, which emits 7 rays 
at 1.17 and 1.33 MeV energy. 

A gamma ray, which scatters in the container wall or 
in the liquid 3 He bath, knocks off a secondary electron 
which gives rise to a heat release in the liquid 3 He bath. 
The interaction probability is roughly proportional to the 
density of the material and thus a large fraction of the 
secondary electrons originate from the quartz glass wall. 
Although the initial energies of the secondary electrons 
cover a broad spectrum, the variation in the effective 
heat release is narrower. Energetic electrons knock off 
additional electrons which produce their own ionization 
tracks. However, once the energy has dropped in the 
regime of several keV, the remaining ionization track is 
short: roughly the rate of energy loss per unit distance 
travelled increases inversely proportional to the electron's 
decreasing energy. Therefore the heating becomes con- 
centrated within a distance < 10 ^m of the final stopping 
point. 

As a source for localized heating in liquid 3 He exper- 
iments, gamma rays suffer from several shortcomings, 
compared to thermal neutrons: 1) It is not possible to 
concentrate on single heat-release events, since many sec- 
ondary electrons can be produced by a single 7 ray. Thus 
there is more variation in these events. For instance, the 
endpoints of different ionization tracks may fall close to 
each other and may merge to produce a large hot bubble 
of variable size. 2) The events are not localized in a well- 
determined location within the container. 3) It is not 
only the liquid 3 He bath which preferentially absorbs the 
heat, but even more so the structural materials of the 
refrigerator. This interferes with the temperature con- 
trol of the measurements: The heating produced by the 
7 radiation in the metal parts of the cryogenic equipment 
may rise to excessive levels, and radiation shielding and 
collimation need to be built into the experiment. 

It was found that qualitatively vortex formation in 
gamma and neutron radiation are similar. But the mea- 
surements also showed that quantitative differences are 
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large: The threshold velocity v cn is substantially in- 
creased and thus the yield N(Q) at a given rotation ve- 
locity f2 in the vortex-free state is smaller. So far no 
careful quantitative measurements have been carried out 
with 7 radiation. 

2.12. Bias dependence of loop extraction 

The dependence of the vortex formation rate N on the 
externally applied bias v is the most tangible quantita- 
tive result from the rotating measurements. It will be 
described and analyzed below (Ruutu et al. 1998a). The 
empirical rate equation (8) is shown to be consistent with 
the KZ model, including the measured value of the pref- 
actor 7. 

2.12.1. Experimental velocity dependence 

The frequency and number of vortex lines can be 
counted from the accumulation record (Figs. 3, 12), when 
the NMR absorption is monitored as a function of time 
during neutron irradiation. Fig. 19 shows the vortex for- 
mation rate N (a), the rate of those neutron absorption 
events N c which produce at least one line and thus be- 
come observable (b), and the average number of lines 
produced by each absorption event (c) . All three quanti- 
ties are determined independently and directly from the 
absorption records. In Fig. 19 the results are plotted as 
a function of the normalized bias v/v cn . To construct 
the plot, the horizontal axis was divided into equal bins 
from which all individual measurements were averaged to 
yield the evenly distributed data points displayed in the 
graphs. 

The rates in Fig. 19a,b are proportional to the intensity 
of the neutron flux, which is contained in the prefactors of 
the expressions given in the different panels of the figure. 
For counting the rates from the NMR absorption records 
it is vital that absorption events do not overlap, when the 
applied counterflow velocity is increased. Therefore three 
different source positions were used which were scaled to 
the same distance using the measured graph in Fig. 8. 

The rates in Fig. 19 increase rapidly with the applied 
bias velocity v: At v/v cn w 4.5, close to the maximum ve- 
locity limit imposed by the spontaneous instability limit 
(Fig. 13), there are almost no unsuccessful (or unob- 
served) absorption events left: N e (oo) w iV e (4.5 v cn ). 
The panel in the middle shows that this saturation cor- 
responds to a flux of 20 neutrons/min. The same flux 
of neutrons was estimated to be absorbed in the sam- 
ple based on independent measurements with commer- 
cial 3 He counters. Also we note that in agreement with 
the scaling properties expressed by the rate equation (8), 
measurements at the two pressures of 2 and 18 bar fall 
on the same universal curves. 




FIG. 19. Rates of vortex line formation plotted vs normal- 
ized bias velocity x = v/v cn : (a) total number of lines formed 
per minute (N), (b) number of observed neutron absorption 
events per minute (N c ), and (c) average number of lines per 
observed event (N ve ). All three quantities have been deter- 
mined independently from discontinuities in NMR absorption, 
ie. each data point is based on one or several accumulation 
plots like those in Fig. 3. The data in the two upper plots are 
proportional to the neutron flux, while the data in the bot- 
tom plot are neutron- flux independent. Solid curves are fits 
to the data, given by the expressions in each panel. Triangles 
are results from numerical simulations described in Sec. 2.14. 
(From Ruutu et al. 1998a). 

Fig. 19 demonstrates that the neutron-induced vortex 
formation process involves a strong stochastic element: 
Close to the critical threshold at v = l.lv cn only one 
neutron capture event from 40 manages to produce a 
sufficiently large vortex loop for spontaneous expansion. 
On increasing the bias flow by a factor of four, almost all 
neutron capture events give rise to at least one escaping 
vortex loop. Secondly, the nonlinear increase of TV as a 
function of v/v cn in the top panel has to be attributed 
to two factors: (i) the increase in the event rate N c in 
the middle panel and (ii) the rapid rise in the number 
of lines produced per event in the bottom panel. These 
conclusions fit with the KZ predictions, as we shall see 
in the next section. 

The most detailed information from the rate measure- 
ments is the dispersion into events in which a given num- 
ber of lines is formed. Fig. 20 displays the rates N ci 
of events which produce i = 1-5 lines. This data ex- 
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FIG. 20. Rates N c i of vortex formation events, grouped ac- 
cording to the number of lines i formed per absorption event, 
plotted vs v/v cn . Triangles are result of numerical simulations 
described in Sec. 2.14. Solid curves are spline fits to simula- 
tion data. (Inset) Threshold velocity v cni /v cn i for the onset 
of an event with i lines, plotted vs the number of lines i. The 
solid curve represents the fit v cn i/v cn i = [2.0 (i — 1) + 
(From Ruutu ct al. 1998a). 

hibits larger scatter, owing to its statistical variation, but 
again after averaging one gets for each value of i a curve, 
which peaks at a maximum, and then trails off. With in- 
creasing value of i the curves shift to successively higher 
velocities. The starting points of each curve, the thresh- 
old velocities v cn i, are plotted in the inset. Their val- 
ues increase monotonically with each consecutive value 
of i. This means that at and immediately above the first 
threshold, v cn = v cn \, only single- vortex events occur. 

A surprising finding from the measurements is the to- 
tal absence of a background contribution to the measured 
rates. A number of tests were performed to look for spon- 
taneous events in the absence of the neutron source, to 
check whether a contribution from the background radia- 
tion level should be subtracted from the measured rates. 
For instance, a vortex-free sample was rotated for 90 min 
at different velocities (0.9, 1.3, and 2.1 rad/s at 2.0 bar 
and 0.94 T c ), but not a single event was noted. 



2.12.2. Analytic model of vortex loop escape 

An analytic calculation can be constructed for indepen- 
dent vortex loops in the applied flow, which explains the 
cubic velocity dependence of the escape rate in Fig. 19. 
Vortex rings, which are part of the random vortex net- 
work after the quench, have a critical radius (Eq. (11)) 
below which they contract and above which they expand 
spontaneously in homogeneous superflow. By comparing 
this radius to the typical curvature in the random net- 
work, we arrive at the desired result (Ruutu et al. 1996a). 

During the quench through the superfiuid transition a 
random vortex network is formed with a characteristic 
length scale of the order of £ v . The later evolution of the 
network leads to a gradual increase in this length. The 
average inter-vortex distance or, cquivalently, the typical 
radius of curvature of the loops increases with time. We 
shall call this length £(i). This "coarse-graining" process 
preserves the random character of the network, in other 
words the network remains self similar or scale invariant. 
Only later a change occurs in this respect, when the loops 
become sufficiently large to interact with the externally 
applied bias flow: This causes large loops to expand, if 
they are oriented transverse to the flow with the correct 
winding direction, while small loops contract, and loops 
with the wrong sign of circulation tend to both contract 
and to change orientation. The time scale of these pro- 
cesses is determined by the magnitude of mutual friction 
between the normal and superfiuid components. In the 
Ginzburg-Landau temperature range the evolution of the 
network occurs in milliseconds and the expansion of the 
extracted loops into rectilinear lines in a fraction of a 
second. 

In superfiuid 3 He the large viscosity of the normal com- 
ponent clamps it to corotation with the container. In 
the rotating frame of reference we may write v n = and 
v s = v. The energy of a vortex loop, which is stationary 
with respect to the walls, is given by (Donnelly 1991) 

£ = E kin + pv , (17) 

where v is the velocity of the applied superflow. In sim- 
ple configurations, the hydrodynamic kinetic energy or 
self-energy of the loop arises from the trapped superfiuid 
circulation with the velocity v S:Vort = (k/2tt) V$, 

Skin = \j p s < vort dV = eL , (18) 

and is proportional to the length L of the loop and its 
line tension, 

Here we neglect the small contribution from the core en- 
ergy, use intervortex spacing £(t) as the upper cutoff, and 
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the superfluid coherence length £(T, P) as the lower cut- 
off for the integration in Eq. (18). This equation is valid 
in the logarithmic approximation, when >• £(T, P). 
While the first term in Eq. (17) is proportional to the 
length L of the loop, the second term involves its linear 
momentum, 



> = J p s v s ,vort dV = p s n J X7$dV = p s nS , 



(20) 



where the last step follows from Gauss's theorem and 
involves the area S of the loop in the direction of the 
normal S/S to the plane of the loop. Thus we write for 
the energy of a loop 



£(L, S, t) = p s K 



4tt £ 



(21) 



where S is the algebraic area, perpendicular to the flow 
and of proper winding direction. This equation expresses 
the balance between a contracting loop due to its own line 
tension, which dominates at small applied velocities, and 
expansion by the Magnus force from the external super- 
flow, which dominates at high applied velocities. The di- 
vide is the equilibrium condition, which was expressed in 
Eq. (11) and corresponds to the situation when the height 
of the energy barrier, which resists loop expansion, van- 
ishes. In this extremal configuration p is antiparallel to 
v, the loop moves with the velocity — v in a frame of the 
superfluid component, but is stationary in the rotating 
frame. In a more general sense, if we consider loops in the 
random network which still deviate from circular shape, 
the extremal case degenerates to a saddle point. This is 
because the extremum requires also minimization with 
respect to deviations from circular shape, ie. the total 
energy is invariant under small variations of the radius 
of the ring, or 8£ = 5E^ n + v • <5p = 0. 

The expansion of the vortex loop should be calculated 
by including the mutual friction forces. In our analytic 
description of vortex loop escape we shall neglect such 
complexity. Instead we shall make use of three scaling 
relations which apply to Brownian networks and are de- 
scribed in more detail in Sec. 2.14. These expressions 
relate the mean values in the statistical distributions of 
the loop diameter T>, area S, and density n to the length 
L of the loop: 

V = AL 5 f 1_<5 , (Aw 0.93, 5 w 0.47), (22) 
\S\ = BV 2 -^ C , (BwO.14, C«0), (23) 



n = CL-P , (C w 0.29, (3 w 2.3) . (24) 

For a Brownian random walk in infinite space the values 
of 6, (3 and C are 1/2, 5/2 and 0. 



The important simplification here is that these rela- 
tions are assumed valid during the entire evolution of 
the network, until sufficiently large rings are extracted 
by the bias flow into the bulk. Using Eqs. (22) and (23), 
we may write Eq. (21) for the energy of a loop in the 
form 



£{V,t) = p s nV 2 



4<(t)A 2 



In 



vB 



■ (25) 



When the scale £(t) exceeds a critical size £ c (v), which 
depends on the particular value of the applied superflow 
velocity v, 



1 



?c(u) ~ A 2 B Attv 



(26) 



the energy in Eq. (25) becomes negative and the loops 
in the network start expanding spontaneously. The to- 
tal number of loops Nb, which will be extracted from 
one neutron bubble, can then be obtained by integrating 
their density from the smallest size £ c to the upper cutoff, 
which is provided by the diameter of the entire network, 
or that of the heated bubble, 2Rb- 



= V h \ 

JL 



dV n{V) 



(27) 



Here the density distribution n(L) = C £~ 3/2 L~ 5/2 , 
combined with that for the average diameter T>(L) = 
A (Li f/ 2 , gives n(V) dV = 2A 3 CV- 4 dV. On inserting 
this into the integral (27) we obtain 



iV b = -irA 3 C 



(28) 



From this equation we see that the requirement 
Nb{v cn ) = returns us the definition of the threshold 
velocity v cn : £ c (t> = v cn ) = 2i?b- This in turn gives us 
from Eq. (26) for the radius of the heated bubble 



Rb = 



1 



A 2 B 8irv c 



In 



2i?b 



Z(T,P) 



(29) 



which we used in Sec. 2.9.1 to derive the temperature de- 
pendence of the threshold velocity: v cn oc (1 — T/Tc) 1 / 3 . 

Eqs. (26) and (29) thus show that ( c oc l/« and i?b oc 
l/t> cn , so that we may write for the vortex- formation rate 
N = (j) n N h from Eq. (28) 



V 
Vcn 



(30) 



where (f> n is the neutron flux. This is the form of the mea- 
sured cubic velocity dependence in the empiric Eq. (8). 
By inserting A w 0.93, C w 0.29 from Eqs. (22) and (24) 
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respectively and <j) n ~ 20 neutrons/min, as determined 
from the saturation of the event rate N c in Fig. 19b, we 
obtain for the rate factor in Eq. (30) 7 w 1.6 min -1 , 
which agrees with the experimental value in Figs. 11 and 
19a. We note that the cubic dependence on the applied 
bias flow in Eq. (30) comes only from the assumption 
that the whole volume of the heated bubble contributes 
equally to the production of vortices, while the values 
of the prefactor 7 and of the constant term in Eq. (30) 
depend on the scaling relations (22-24). 

The definition of the threshold velocity i> cn i, which 
applies for an event in which i rings are formed simul- 
taneously, is roughly consistent with the requirement 
Nb(v = v cni ) w i. This gives v cni /v cn — i 1/3 , which 
agrees with the measured result in Fig. 20. 

To summarize, we note that the KZ model, combined 
with the simplest possible interpretation for the loop es- 
cape from a random vortex network, reproduces both the 
measured cubic dependence on the normalized velocity, 
(i>/v cn ) 3 , and the magnitude of the extraction rate. 



2.13. Neutron-induced vortex formation at low 
temperatures 

The low temperature properties of superfluid 3 He could 
not have been extrapolated from measurements restricted 
to temperatures above 0.80 T c , even within the frame- 
work of the general theory. Similarly, to understand the 
neutron-induced vortex formation process, it is impor- 
tant to extend the measurements to lower temperatures. 
One such study has been performed with NMR down to 
0.4 T c in a rotating sample (Finne et al. 2004a). This 
measurement is a continuation of the NMR work, which 
we have discussed so far above 0.80 T c . Another series of 
measurements has been performed calorimetrically below 
0.2 T c in a quiescent bath (Bauerle et al. 1996, 1998a). At 
these very low temperatures quantized vortex lines have 
a long life time in zero flow. Thus there is sufficient time 
after the neutron absorption event to detect the vortices. 

NMR measurement in 3 He-B loses rapidly in ampli- 
tude resolution on cooling to lower temperatures: the 
susceptibility drops with decreasing temperature and the 
width of the NMR spectrum increases. Thus the change 
of NMR absorption per one rectilinear vortex line quickly 
diminishes and single- vortex resolution is lost. Also be- 
low 0.6 T c a new phenomenon enters, namely superfluid 
turbulence (Finne et al. 2003). If here a few seed vor- 
tex loops are injected into rotating vortex-free superflow, 
they become easily unstable and form a turbulent vor- 
tex tangle. In the turbulent state the number of vortices 
rapidly multiplies and the final state becomes an equilib- 
rium array of rectilinear lines. 

In neutron-irradiated superflow this is exactly what 
happens: First vortex loops are extracted from the neu- 
tron bubble and injected into the bias flow. Their number 
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FIG. 21. 3 He sample with NMR measuring setup. The sam- 
ple is contained in a fused quartz tube which has a diameter 
of 6 mm and length 110 mm. This space is separated from the 
rest of the liquid 3 He volume with a partition disc. In the 
center of the disc an orifice of 0.75 mm diameter provides the 
thermal contact to the liquid column which connects to the 
sintered heat exchanger on the nuclear cooling stage. Two su- 
perconducting solenoidal coil systems with end-compensation 
sections produce two independent homogeneous field regions 
with axially oriented magnetic fields. An exterior niobium 
cylinder provides shielding from external fields and additional 
homogenization of the NMR fields. The NMR magnets and 
the Nb shield are thermally connected to the mixing cham- 
ber of the pre-cooling dilution refrigerator and have no solid 
connection to the sample container in the center. The two 
split-half detection coils are fixed directly on the sample con- 
tainer. (From Finne et al. 2004a). 



can be tuned to some extent (in the range 1-5), by choos- 
ing the value of the bias velocity (Fig. 19c). Then with 
some probability, which increases towards lower temper- 
atures and higher bias velocities, turbulence starts. As 
a result the sample is suddenly filled with the equilib- 
rium number of vortex lines, which means that instead 
of a few rectilinear lines the final state contains of order 
~ 10 3 lines. In these conditions the final state has lit- 
tle connection with the loop extraction process and the 
details of the neutron-induced vortex formation become 
lost. But we can use this method to start turbulence in 
a rotating superfluid column. 
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2.13.1. Experimental techniques 

The rotating NMR measurements at low temperatures 
are performed in the setup of Fig. 21. For the study 
of turbulence it proved fortunate that this arrangement 
includes two NMR detection coils at both ends of the 
long sample cylinder. A second bonus is the relatively 
high critical velocity of the container so that vortex-free 
superflow could be maintained to above 3.5rad/s below 
0.8 T c at 29.0 bar. 

On cooling to lower temperatures vortex-free superflow 
generally becomes more and more difficult to generate 
and maintain. The lower section of the 3 He volume be- 
low the orifice is directly connected with the sintered heat 
exchanger. This region is flooded with remanent vortices 
from the porous sinter already at low rotation. Generally 
one finds that these vortices start to leak through the 
orifice into the sample volume on cooling below 0.6 T c . 
Surprisingly the sample tube, with which the neutron 
measurements were performed, seemed immune to this 
problem in spite of the fact that the orifice had a rela- 
tively large diameter of 0.75 mm. 

Below 0.5 T c remanent vortices become a problem. 
Here successive accelerations to rotation have to be sepa- 
rated by extensive waiting periods at stand still, to allow 
the slow annihilation of remanent vortices. If rotation is 
started too soon, remanent vortices are still present and 
become easily unstable in rotation. As a result vortex- 
free rotation at velocities above lrad/s is then not pos- 
sible to achieve. Just below 0.60 T c the waiting time is 
around 5min in the container of Fig. 21, but at 0.40 T c 
it is found to be 30min or more. No procedure using 
rotation in different directions or with different decelera- 
tion rates seemed to help in cutting down on the waiting 
period. 

The experimental setup in Fig. 21 is equipped with two 
independent continuous- wave NMR spectrometers. Each 
spectrometer includes a split- half excitation/detection 
coil, wound from thin superconducting wire. The two 
coils are installed at both ends of the sample cylinder. 
Each coil is part of a high-Q tank circuit, with a reso- 
nance frequency approaching 1 MHz and Q~ 10 4 . Fig. 22 
shows a number of absorption spectra measured with the 
bottom spectrometer, similar to the high-temperature 
examples in the left panel of Fig. 6. The spectra illus- 
trate the relevant NMR line shapes and their large width 
at low temperatures, with varying numbers of rectilin- 
ear vortex lines in the sample. Again the peak height of 
the CF maximum can be calibrated to yield the number 
of lines N in the linear limit, when N <C iV cq . Often 
enough, to determine N, one has to increase fi to some 
higher reference value where the sensitivity of the CF 
peak is restored. 

Fig. 23 shows an example of a calibration where the CF 
peak height is measured as a function of N (or equiva- 




Magnetic field sweep in the bottom spectrometer (mT) 
FIG. 22. NMR absorption spectra of 3 He-B in rotation at low 
temperatures. As in the left panel of Fig. 6, changes in the 
line shape can be calibrated to give the number of vortex lines 
N. Here N is characterized by the rotation velocity Q V (N) at 
which a given number of lines N is in the equilibrium state. 
The different spectra have been measured with the RF exci- 
tation at constant frequency /, using a linear sweep of the 
axially oriented polarization field H. The spectra have been 
recorded at constant temperature and thus all have the same 
integrated total absorption. The sharp absorption maximum 
at low field is the counterflow peak (CF). Its shift from the 
Larmor field (at 21.02 mT) is used for temperature measure- 
ment. When a central cluster of rectilinear vortex lines is 
formed, the height of the CF peak is reduced. In the equilib- 
rium vortex state (fi = f2 v ), where the number of vortex lines 
reaches its maximum, the spectrum looks very different: it 
has appreciable absorption at high fields and borders promi- 
nently to the Larmor edge. This spectrum is more similar to 
that of the non-rotating state (Q = 0). As shown in Fig. 23, 
when the vortex number is small, Q v <C Q, the reduction in 
the CF peak height can be calibrated to give fi v and thus N. 
(From Finne et al. 2004a). 

lently Q v ). The measurement has been performed for 
small vortex clusters (N < 100), by measuring each cali- 
bration point in a four step process: 1) First the spectrum 
is recorded in the vortex-free state in the reference con- 
ditions and the CF peak height A (Q rc f) is obtained. 2) 
Next a large number of vortex lines is created by irradi- 
ating with neutrons (with the rotation at O re f or higher) . 
3) The sample is then decelerated to the low rotation ve- 
locity fiiow ^ref so that part of the vortex lines are 
observed to annihilate. Since the long cylinder is ori- 
ented along the rotation axis only within a precision of 
~ 0.5°, this means that any annihilation barrier must be 
negligible or that at fli ow the sample is in the equilibrium 
vortex state f2i ow = Sl v (N), with a known number of vor- 
tex lines N. 4) Finally the rotation is increased back to 
the reference value fi rc f and the spectrum is recorded in 
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Number of vortex lines in the cluster N 




Cluster size £l v (rad/s) 
FIG. 23. Calibration of CF peak height A versus vortex line 
number N. The CF peak height in the vortex-free state, 
j4o(fi r ef), is compared to the peak height A(fl re f,Q v (N)), 
which is measured at the same rotation velocity f2 re f for a 
small vortex cluster of size fi v , which is prepared as de- 
scribed in the text. The quantity plotted on the vertical 
scale is the relative reduction in peak heights, AA/Aq = 
[Ao(firof) — A(fl TC {,£l v (N))]/Ao(£l TC f), measured in constant 
conditions. The solid line is a fit, given by the expres- 
sion in the figure. The conversion from O v (bottom axis) 
to N (top axis) in the continuum picture is (Ruutu et al. 
1998b) N = JV (1 - d c JR) 2 , where N = % R 2 (2Q, V / n) and 
d cq = [K/(8Tvn v ) ln(K/27rO v r^)] 1/2 . Here r c ~ £(T,P) is the 
radius of the vortex core. Using this conversion it is found 
that AA/Aq is a linear function of N, similar to the right 
panel in Fig. 6. (From Finne et al. 2004a) . 

order to measure the CF peak height A(£l ie f, Q V (N)). In 
Fig. 23 the reduction AA = A (Q Te {) — A(Q ie {, Ct v ) in the 
CF peak height of the two spectra is plotted as a function 
of f2 v . To reduce the dependence on drift and other irreg- 
ularities it is normalized to the CF peak height A (£l Te f) 
of the vortex-free reference state. In contrast to the cal- 
ibration in the right panel of Fig. 6, this procedure does 
not rely on single- vortex resolution. However like always, 
it does require that the order parameter texture is stable 
and reproducible and that new vortices are not formed 
during rotational acceleration. 

The result in Fig. 23 is a smooth parabola, which be- 
comes a straight line if AA/Ao is plotted as a function 
of the number of rectilinear lines N. Any sample with 
an unknown number of lines, which is less than the max- 
imum calibrated number, can then be measured in the 
same reference conditions (f2 re f , T, and P) and compared 
to this plot, to determine N . In Fig. 24 a measurement 
of vortex formation in neutron irradiation is shown for 
which the calibration was used. After each irradiation 
session at different bias flow velocity v = HR the rotation 
is changed to O re f and the NMR absorption spectrum is 
recorded. From this spectrum the reduction in CF peak 
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FIG. 24. Rate of vortex formation in neutron irradiation. 
The average number of rectilinear vortex lines created per 
unit time during the irradiation period is shown as a function 
of the rotation velocity f2. The data are fitted with the ex- 
pression N = 0.336 [(fi/1.43) 3 - 1] 

min (with Q, in rad/s). 
Depending on the rate N, the irradiation time varies here from 
30 min to 4.5 h, so that the number of accumulated vortices 
remains within the range of the calibration in Fig. 23. This 
calibration is used to determine the number of vortices from 
the relative reduction in the CF peak height, AA/Ao- The 
distance of the neutron source from the sample was d = 34 cm. 
The range of the bias flow is limited in these measurements 
between the critical velocity £l cn — 1.43 rad/s and the upper 
limit fi = 3.0 rad/s, where the turbulent events start to occur. 
(From Finne et al. 2004a). 

height is determined, by comparing to the spectra of the 
vortex-free state which are measured regularly between 
neutron irradiation sessions. 

2.13.2. Measurement of vortex formation rate 

In Fig. 24 the rate of vortex formation TV is measured 
at 0.53 T c as a function of the bias velocity v — QR, dur- 
ing neutron irradiation. The result supports the cubic 
rate equation (30). The critical velocity v cn — Q cn R — 
4.3 mm/s is in the same range as the data in Fig. 13. This 
is plausible if the critical velocity is determined only by 
the size of the neutron bubble which does not change ap- 
preciably with decreasing temperature, Eq. (14). How- 
ever, the rate factor 7 is 40 times smaller than in the 
measurements above 0.80 T c and below 21.2 bar (when 
scaled to the same neutron flux). A second similar mea- 
surement at 10.2 bar and 0.57 T c gives a critical velocity 
v cn = 3.6 mm/s and a rate constant 7 which is 3 times 
larger than the measurement in Fig. 24. 

These results raise the question how do v cn and 7 vary 
over a wider temperature and pressure range. In partic- 
ular this becomes a problem at high pressures, where a 
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wide range of stable A phase exists between T c and the 
ambient B-phase bath temperature To, where the irradi- 
ation is performed. In the 10 bar measurement A phase 
is not stable. In both of these measurements the avail- 
able range of bias velocities is limited (to approximately 
v cn < v < 2v cn ), by the appearance of turbulent events 
at higher velocities. Both of these low temperature mea- 
surements are consistent with the cubic form of the rate 
equation, but do not conclusively prove it owing to the 
limited velocity range (unlike the high-temperature re- 
sult in Fig. 19). Clearly more extensive measurements as 
a function of pressure and temperature are called for. 



2.13.3. Superfluid turbulence in neutron irradiation 

At high rotation and low temperatures below 0.60 T c 
neutron irradiation events may become turbulent, simi- 
lar to vortex formation from other sources (Finne et al. 
2003). This means that vortex loops, which have been 
extracted from the neutron bubble and are injected into 
the bias flow, may start to interact, to produce a vortex 
network of large scale. This tangle then blows up and fills 
the rotating sample with the equilibrium number of vor- 
tex lines. Some NMR characteristics of such a neutron 
capture event are illustrated in Fig. 25. The basic exper- 
imental feature here is that the sample is suddenly filled 
with the equilibrium number of rectilinear vortex lines, 
apparently as a result of one neutron capture event: The 
NMR absorption spectrum jumps (via a brief transitory 
period) from a line shape with a large CF peak to that 
of the rotating equilibrium state of totally different form 
(Fig. 22). 

Turbulent events are observed in neutron irradiation 
only if (i) the rotation velocity exceeds fi cn and (ii) the 
temperature is sufficiently low so that vortex motion is 
not heavily damped by mutual friction. Even then these 
processes are stochastic such that the vortex loops ex- 
tracted from the neutron bubble only rarely achieve the 
proper initial conditions in which turbulent loop expan- 
sion starts to evolve. With decreasing temperature and 
increasing rotation the probability of turbulent events 
increases. This is understandable for the following rea- 
sons: (i) With decreasing temperature mutual friction 
damping is reduced, Kelvin wave excitations on the ex- 
isting vortices grow in amplitude, more new loops are 
formed, and via reconnection processes the intersecting 
loops multiply to a turbulent cascade, (ii) With increas- 
ing bias velocity the number of vortex loops, which are 
injected into the bias flow, rapidly increases (Fig. 19, bot- 
tom) and the probability of their intersections increases. 
Here we shall not discuss general features of turbulent 
superfluid dynamics, only describe the main signatures 
of the neutron-induced turbulent events. 

Fig. 25 illustrates how the height of the CF peak sud- 
denly drops to zero in the NMR spectrum when a tur- 
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FIG. 25. NMR signatures of turbulent vortex formation in 
neutron irradiation. The CF peak height (Fig. 22) is moni- 
tored as a function of time in the two NMR detection coils. 
Here a vortex tangle suddenly starts to expand, first in the 
bottom coil (at t = 0), when its CF peak rapidly drops to 
zero. After a delay of 70s, which the turbulent front needs 
to travel to the lower edge of the top coil 90 mm higher along 
the column, a similar collapse is recorded by the top spec- 
trometer. The flight time r = z/(aflR) corresponds to the 
situation where z = 95 mm and turbulence is first formed in 
the middle of the bottom coil. This sequence of events is 
schematically displayed here at fi = 1.61rad/s. The sudden 
collapse of the CF peak means that the turbulent tangle is 
rapidly polarized at these large flow velocities and that the 
global counterflow between the normal and superfluid com- 
ponents is thereby removed. Measurements of this type show 
that from the initial injection site, where the extracted vor- 
tex loops first start to intersect in the bias flow, turbulence 
expands by forming two fronts which move at constant veloc- 
ity towards the top and bottom ends of the rotating column. 
(From Finne et al. 2004a). 



bulcnt event starts to evolve. In this schematic example 
the turbulent tangle first appears inside the bottom coil. 
The collapse of the CF peak height is the first feature 
in the NMR absorption spectrum which signals the tur- 
bulence: Its rapid decay shows that the turbulent state 
becomes polarized, to mimic on an average solid-body 
rotation. Simultaneously the NMR absorption intensity 
from the CF peak is transferred close to the Larmor edge 
of the spectrum where a new sharp peak rapidly grows 
in intensity. This new peak then slowly decays to the 
line shape of the equilibrium state (Fig. 22). The inten- 
sity in the Larmor region reflects how the vortex density 
evolves within the coil: It first rapidly overshoots to a 
value which is about twice that in equilibrium and then 
slowly rarefies to the equilibrium value. 

The decaying CF signals from the two coils do not 
overlap in Fig. 25. Even the more slowly relaxing over- 
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FIG. 26. Turbulent expansion of quantized vorticity in a ro- 
tating superfluid column. Initially the column rotates in the 
vortex-free state. If sufficiently many loops are extracted from 
a neutron bubble, they may interact to develop a superfluid 
tangle. The tangle forms with a sharp advancing front and a 
trailing tail which decays into the equilibrium rotating state 
with rectilinear lines. Two turbulent fronts thus propagate 
along the column in opposite directions, while in the mid- 
dle there exists already an ordered array of rectilinear vortex 
lines, if the column is sufficiently long. The axial velocity of 
the fronts is controlled by the dissipative mutual friction a. 
The ordered array in the middle remains fixed in the rotat- 
ing frame, while the two fronts rotate, as determined by the 
reactive mutual friction a! . 



shoots in the Larmor region do not overlap in time when 
T < 0.5T C and turbulence is initiated at one end of the 
sample tube. This means that turbulence propagates 
along the column as a stratified layer: By the time the 
NMR absorption in the top coil gives the first indica- 
tion of the approaching turbulent front, the bottom coil 
has already settled into its stable equilibrium line shape. 
The turbulent front is formed from a relatively thin layer 
of disordered tangle. Here the polarization of the circu- 
lation reaches its final equilibrium value almost immedi- 
ately, even before the vortex density has reached its peak, 
and well before the vortex configuration has stabilized to 
its equilibrium state of rectilinear lines. These observa- 
tions suggest that vorticity propagates in the rotating 
column in the peculiar configuration shown in Fig. 26. 
Here turbulence does not fill the entire column, only short 
sections, before it straightens into rectilinear lines. Such 
a configuration has not been described before - simply 
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FIG. 27. Turbulent vortex formation in neutron irradiation at 
29.0 bar pressure. The sample is irradiated in constant con- 
ditions until a turbulent vortex expansion event takes place. 
The irradiation time required to achieve the turbulent event 
is measured. Results from measurements at two different con- 
stant conditions are shown in this plot in the form of cumu- 
lative probability distributions. The solid curve is a fit of 
the data measured at 0.53 T c to a distribution function of the 
form P(t) = f — (I -pn)^ 1 , where the probability of a single 
neutron absorption event to start turbulence is p n — 0.092. 
The rate iV c of the neutron-induced vortex injection events is 
taken as 0.55 N from Fig. 24 as appropriate for v/v cn — 2.32 
(see Fig. 19). There were no cases among the two sets of mea- 
surements where a turbulent event would not have been ob- 
served after an irradiation period of 20 min. The inset shows 
the axial location of the initial turbulent seed, the site of the 
neutron capture, measured from the orifice upward with a 
technique explained in Fig. 25. The corresponding irradia- 
tion time needed to achieve the event is shown on the vertical 
axis. As expected, the sites are randomly distributed along 
the sample. 



because it has not been possible to set up this situation 
in He-II, ie. one where vortices are injected in high- 
velocity rotating superflow. 

The turbulent front moves along the column with 
a fixed velocity which was measured by Finne et al. 
(2004b). This velocity has the same value as that of 
a single short vortex filament moving along the cylinder 
wall in the initial vortex-free bias flow: v z = aflR. Here 
a is the dissipative mutual friction coefficient which was 
measured by Bevan et al. (1997a). This means that from 
the delay between the signals of the two detector coils (as 
marked in Fig. 25), we may calculate the axial location 
z where the turbulent event started. This location has 
been plotted for the data in Fig. 27 in the inset of this 
figure. 

The probability of a neutron capture to trigger a tur- 
bulent process is studied in Fig. 27. Two measurements 
are shown, at different temperatures and rotation veloc- 
ities. The cumulative probability distribution has been 
plotted for the irradiation time needed to achieve the 
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turbulent event. The sample is irradiated at constant 
conditions until the CF peak is observed to collapse sud- 
denly. The irradiation time is plotted on the horizontal 
axis. On the vertical axis the number of turbulent events 
observed within this time is shown, normalized to the to- 
tal number of events: 13 events in one case (at 0.53 T c ) 
and 9 in the second (at 0.45 T c ). As seen from the plot, 
the number of events is insufficient to produce smooth 
probability distributions, but the irradiation times are 
observed to be distributed over the same range in the 
two cases, ie. their distributions have similar average 
values and widths. This in spite of the fact that the mea- 
surements at 0.45 T c were performed at half the rotation 
of those at 0.53 T c . At high temperatures a reduction 
by two in velocity results in a significant decrease in the 
yield of vortex loops from one neutron absorption event 
(Fig. 19). The fact that the two distributions in Fig. 27 
do not differ significantly indicates that with decreasing 
temperature the transition to turbulence becomes more 
probable and less sensitive to the initial configuration of 
the injected loops. 

Let us consider the measurement at 0.53 T c in more 
detail. The state of the sample changes during the ir- 
radiation at constant O: (i) the initial state is vortex- 
free, (ii) during the irradiation rectilinear vortex lines 
are formed at a rate which can be extrapolated from 
Fig. 24, (iii) until finally all vortex-free CF is completely 
terminated in a turbulent event. In the final step the 
state of the sample changes from one with only a small 
central vortex cluster to one with the equilibrium num- 
ber of vortex lines (N eq <~ 2600). Note that to observe 
a new turbulent event the existing vortices have to be 
annihilated, by stopping rotation. Then the vortex-free 
state can be prepared again and a new irradiation ses- 
sion can be started. The longest irradiation time is 
here ~ 20min. During this period the cluster grows 
at the rate N = 3.9 vortices/min, so that it contains 
about 80 vortices when the turbulent event finally starts. 
At this point the CF velocity at the sample boundary 
v = v n — v s = QR — kN/ (2-kR) has been reduced by 
2.7% from the initial state. Since the mean irradiation 
time in the measured distribution is only ~ 250 s, the 
reduction in CF velocity by vortices formed before a tur- 
bulent event is minor. We may thus view the result in 
Fig. 27 as representative of these particular values of ro- 
tation and temperature. 

One may wonder whether a turbulent process results 
from a single neutron capture event or from the coin- 
cidence of two or more events. In the latter case the 
simultaneous events need to be sufficiently close not only 
in time, but also in space, so that expanding loops, which 
have been extracted from the two neutron bubbles, have 
a possibility to intersect. (The probability of the bubbles 
themselves to intersect is very small.) The intersection 
of the loops expanding from the two random positions 
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FIG. 28. Bolometer box, which is used to measure calorimet- 
rically the energy balance of the neutron capture reaction in 
the low temperature limit of ballistic quasiparticle motion. 
Two vibrating superconducting wire loops are included, one 
for measuring the temperature of the liquid within the box 
and the other as a heater to calibrate the calorimetric mea- 
surement. (From Bauerle et al. 1998b). 

along the height of the cylindrical sample is more likely 
to occur closer to the middle than at the top or bottom 
end. However, the events listed in the inset of Fig. 27 
occur randomly along the sample. In all cases the NMR 
signatures from the turbulent events are similar, there is 
no prominent variation in their appearance depending on 
where the event starts. Additionally, at lower tempera- 
tures, the turbulent events occur close above Q cn . Here 
successful neutron absorption events, which lead to the 
extraction of loops to the bulk, are rare but, neverthe- 
less, turbulence becomes more probable. From this we 
presume that a single neutron capture event in suitable 
conditions must be able to start the turbulence. For a 
more careful proof of this point the measurements should 
be repeated as a function of the neutron flux. 

To conclude, with deceasing temperature turbulence 
moves closer to the critical threshold v cn . This means 
that the range of bias velocities, in which the neutron- 
induced vortex formation process can be studied, de- 
creases and ultimately below 0.45 T c a measurement of 
the bias dependence becomes obsolete. Nevertheless, 
here and at lower temperatures it is a practical technique 
to start turbulence. 

2.13.4. Calorimetry of vortex network 

Another type of low temperature experiments are the 
calorimetric measurements of Bauerle et al. (1996, 1998a) 
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which determine the amount of heat which is dumped 
into the liquid 3 He-B bath in each individual neutron 
absorption event. They allow quantitative comparison 
to the KZ model. These measurements are performed in 
the very low temperature limit of ballistic quasiparticlc 
motion in a quiescent bath, with zero bias flow. 

The measuring probe is a superconducting wire loop. 
It is oscillated in the liquid with a frequency of a few 
hundred Hz and a high Q value, by driving the loop at 
resonance with an ac current in a dc magnetic field ori- 
ented perpendicular to the wire. The damping of the wire 
oscillations measures the density of quasiparticle excita- 
tions in the liquid and can be calibrated to give the tem- 
perature. For this a second vibrating wire loop is used 
as a heater (Fig. 28) . The latter is excited with a known 
current pulse to heat up the liquid, by driving the wire at 
supercritical velocities where the break-down of Cooper 
pairs gives rise to a shower of quasiparticles. To contain 
the ballistically moving quasiparticles, both wire loops 
are placed inside a box from which the particles leak out 
through a pin hole at a well-known rate (Bauerle et al. 
1998b). 

The thermal connection between the liquid in the 
calorimeter box and that in the surrounding 3 He-B bath 
is only via the small pin hole. A neutron absorption event 
inside the box heats up the liquid and a thermal pulse 
is recorded with the wire resonator. The rise time is de- 
termined by the resonator properties, while the trailing 
edge (with a time constant of about 1 min) monitors the 
much slower leakage of quasiparticle excitations from the 
box through the orifice. If energy is released into the box 
on a still longer time scale it is not recorded in the form 
of pulses. 

It is found that neutron absorption events amount to 
roughly 100 keV smaller thermal pulses than the 764 keV, 
which a slow neutron is expected to yield for the absorp- 
tion reaction with a 3 He nucleus. Since the various re- 
combination channels of the ionized charge and the sub- 
sequent thermalization processes in liquid 3 He are poorly 
known, it is not quite clear how large an energy deficit 
one should expect. Obvious losses include the ultraviolet 
radiation absorbed in the walls of the bolometer box and 
the retarded relaxation of excited electronic states of he- 
lium atoms and molecular complexes. However, Bauerle 
et al. (1996, 1998a) expect that these contributions are 
not strongly pressure dependent and on the order of 7% 
of the reaction energy. They then ascribe the remain- 
ing energy deficit to the random vortex network which is 
created in the neutron bubble and which in the low tem- 
perature limit should have a long life time, when mutual 
friction approaches zero and no flow is applied. 

In Table I the measured missing energy at three dif- 
ferent liquid pressures has been compared to that esti- 
mated from the KZ model. The measured result A.E exp 
is recorded on the top most line while the calculated com- 



parison proceeds stepwise from one line to the next in the 
downward direction of the table. The final result Ai? t hoor 
can be found on the lower most line. The thermal diffu- 
sion constant, D = 3k^/C v = v^tt , is derived from the 
conductivity fcx, which has been tabulated by Whcatlcy 
(1975) while all other liquid 3 He values are taken from 
Greywall (1986). The agreement between the top and 
bottom lines is within a factor of 2, which is surprising, 
if we remember the uncertainties which are built into this 
comparison. Table I has been included in this context to 
draw attention on the magnitudes of the different quan- 
tities. 

The significance of the calorimetric measurement rests 
entirely on the quantitative analysis of the measured data 
- whether an energy deficit is present and can be ascribed 
to vortices. The quantitative analysis is hampered by 
three kinds of uncertainties: 

(1) The proportions of the radiative and retarded con- 
tributions in the thermalization after neutron absorption 
are not known sufficiently well (Leggett 2002). 

(2) To calibrate the bolometer box, the heater wire is 
vibrated at a velocity which exceeds the pair-breaking 
limit. It is now known that when the vibrating wire loop 
reaches this limit, then a beam of quasiparticles and vor- 
tex loops is created. The fraction of energy spent on gen- 
erating vortices is temperature dependent and appears to 
decrease rapidly with increasing temperature. With in- 
creasing drive level more and more vortices are formed at 
different points along the wire loop (Fisher et al. 2001; 
Bradley et al. 2004). The calibration of the bolometer 
box then becomes uncertain, if a sizeable fraction of the 
energy of the calibrating electrical pulse ends up in a long 
lived turbulent vortex network. 

(3) The decay time of a vortex tangle is unknown in 
the zero temperature limit. The calorimetric measure- 
ment of the energy deficit in the neutron absorption event 
requires a time constant of tens of seconds so that the de- 
cay of the tangle would not contribute to the measured 
thermal pulse. On the other hand, the measurements of 
Fisher et al. (2001) on vibrating-wire-generated tangles 
suggest that their time constant of decay is seconds. The 
simulation calculations of Barenghi and Samuels (2002) 
suggest even faster decay. However, the decay time is 
expected to depend on the topology and line density of 
the vortex tangle since it proceeds via reconnection into 
loops which fly away in the quiescent bath and finally 
annihilate on solid surfaces (where part of the energy 
may go directly to the wall, escaping the liquid). In 
this respect the random vortex network created in the 
neutron absorption event (Sec. 2.14.1) is expected to be 
quite different from the vortex tangle generated by the vi- 
brating wire: According to the measurements of Bradley 
et al. (2004), during stationary state vortex generation 
and decay by a vibrating wire the line density reaches 
only a very low value of 20mm~ 2 (an inter- vortex dis- 
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TABLE I. Comparison of measured and estimated total vortex line energies at different pressures in a random vortex network, 
as generated by a neutron absorption event according to the KZ model. The ambient temperature of the 3 He-B bath is taken 
to be To = 0.16 mK. For the value of tq we use the time it takes for the normal phase bubble to disappear. The coefficient 
ai = 2.1, for estimating the total vortex-line length L v , is taken from the simulation results in Sec. 2.14.1, when the value of 
the dimensionless bubble diameter is D b — 200. 



tance <~ 0.2mm). In comparison, in Table I the initial 
line density expected in the neutron bubble, as precipi- 
tated by the KZ process, is of order 10 7 mm~ 2 . 

The original interpretation of the calorimetric mea- 
surements supports the KZ model of vortex formation. 
However, to keep this statement valid in view of recent 
developments on vortex dynamics at ultra low tempera- 
tures, more work is required. In principle a calorimetric 
measurement of the energy in the random vortex net- 
work, which is generated in a neutron absorption event, 
provides more direct proof of the KZ mechanism, than 
rotating measurements: Here in the absence of applied 
flow, the discussion about the superflow instability at the 
neutron bubble boundary is irrelevant. Also the calori- 
metric measurement gives directly the domain size £ v of 
the inhomogencity in the order parameter distribution 
(Eq. (3)) right after the quench, which can otherwise be 
inferred only indirectly from other types of measurement. 

It is interesting to note that the understanding gained 
from the calorimetric work has turned 3 He-B into an at- 
tractive absorber material in the search for dark matter 
particles. In such a dark matter detector the collision en- 
ergy will be measured with an array of bolometer boxes 
equipped with micro-mechanically fabricated vibrating 



resonators (Collin et al. 2004). In the next sections we 
shall analyze further the initial loop formation in the neu- 
tron bubble and its later evolution, which we so far have 
omitted. 

2.14. Simulation of loop extraction 

In spite of the agreement, which we have listed so far 
between the measured characteristics of neutron-induced 
vortex formation and the KZ model, a deeper under- 
standing of the processes involved would be important. 
Owing to its phcnomenological content, the KZ model 
is based on general concepts and contains few media- 
dependent parameters (such as the supcrfiuid coher- 
ence length £ or the order-parameter relaxation time r) . 
These are known in the case of supcrfiuid 3 He, hence the 
predictions of the model for the initial state of the vortex 
network can be calculated. In the rotating experiments 
the initial state is connected with observable quantities 
only through the complex evolution of the network, gov- 
erned by superfiuid hydrodynamics. This means that a 
more rigorous comparison requires numerical simulation. 

A number of numerical simulations exist on the cvolu- 
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tion of a network of linear defects. These apply to cosmic 
strings, liquid crystals, and vortices in supcrfluid 4 Hc. 
Such results cannot be directly transferred to neutron- 
induced vortex formation in 3 He-B. The differences with 
the cosmic string and liquid crystal calculations arise 
from the different boundary conditions and the equations 
which govern the evolution of the network in the applied 
external bias fields. In the case of superfluid 4 He, stud- 
ies of random vortex networks often concentrate on vor- 
tex flow driven by thermal counterflow between the nor- 
mal and superfluid components in a stationary situation 
(Tough 1982). In contrast, the vortex network, which is 
produced in 3 He-B in a neutron absorption event, is in 
a state of rapid non-stationary evolution. Nevertheless, 
standard techniques exist (see eg. Schwarz 1978, 1985, 
1988) and can be applied for solving also the transient 
problem in 3 He-B. Here we describe calculations which 
address the dependence of the vortex-formation rate on 
the normalized bias velocity x = v/v CD (Ruutu et al. 
1998a). 

2.14.1. Initial loop distribution 

Vachaspati and Vilenkin (1984) developed an approach 
for the simulation of the initial network of linear defects 
after a rapid phase transition, adapted to the case of 
cosmic strings. This technique has been used in most 
studies since then (see e.g. the reviews by Hindmarsh 
and Kibble 1995 and Bray 1994). 

We approximate the neutron bubble with a cubic vol- 
ume which is subdivided into smaller cubes, such that 
the size of these volume elements equals the length scale 
of the initial inhomogencity in the order-parameter dis- 
tribution at the moment of defect formation. This length 
is on the order of the coherence length and in the simula- 
tion it plays the role of a "unit" length: It is a parameter 
of the model, which has to be derived by other means. In 
the case of the superfluid transition this length is given 
by Eq. (3). 

An arbitrary phase is assigned to each vertex of the 
grid, to model the initial random inhomogeneity of the 
order parameter. It is usually assumed that the distribu- 
tion of the phase in each segment of the grid between two 
vertices corresponds to the shortest path in the phase cir- 
cle. Thus it is possible to determine whether a line defect 
pierces a given face. Then the centers of the correspond- 
ing faces are connected to form closed or open (depending 
upon specific boundary conditions) linear defects, strings 
in the cosmological case and vortices in the superfluid. 

Vachaspati and Vilenkin assigned to each vertex a 
value of the phase from the following set: {0, 27r/3, 47r/3} 
and then studied the statistical properties of the result- 
ing network of strings. They found that most (70%) of 
the strings were in the form of open segments, which ex- 



tended from one boundary of the system to another. For 
closed loops they found two scaling relations to hold: 

n = CL-P, (31) 

and 

V = AL S , (32) 

where w § , S w ^ , n is the density of loops with a given 
length L, and V is the average spatial size of a loop. It is 
usually defined as an average of straight-line dimensions 
in x, y and z directions. In this model both the charac- 
teristic inter-vortex distance and the radius of curvature 
are on the order of the length scale of the spatial inho- 
mogeneity (i.e. the size of the small cubes, which here 
has been set equal to unity). Later, other variations of 
this model have been studied, including other types of 
grids and other sets of allowed phases. However, it has 
been found that the scaling relations (31) and (32) hold 
universally in each case. 

The direct applicability of these results to the case of 
vortex formation in 3 He-B is not evident. In cosmology 
the open lines are the most significant ones: only these 
strings may survive during the later evolution if no bias 
field exists to prevent closed loops from contracting and 
annihilating. In the case of a normal-liquid bubble within 
the bulk superfluid there exists an obvious boundary con- 
dition: the phase is fixed at the boundary and there are 
no open lines at all. Thus it is of interest to find out 
the influence of the boundary conditions on the scaling 
relations (31) and (32). 

Another question, which arises in the case of superfluid 
vortices, concerns the interaction of the loops with the ex- 
ternal bias field due to the normal-superfluid counterflow. 
The energy of a tangled vortex loop in the counterflow 
is proportional to its algebraic area S = §ydz in the di- 
rection x of the counterflow velocity, Eq. (21). Thus the 
dependence of S on the length of the loop is of interest 
as well. 

The simulation is performed in a cubic lattice with a 
fixed (zero) phase at the outer boundary. The set of 
allowed phases is not limited. It is known that such lim- 
itations affect the number of open lines in cosmic-string 
simulations, but this should be irrelevant for the network 
in the neutron bubble. Several vortices may pass through 
one cell in the grid. In this case the corresponding faces 
of the cell are connected at random. For calculating the 
length of a loop, it is assumed that both straight and 
curved segments of the loop inside one cell have unit 
length. The size of a loop in the direction of a specific 
coordinate axis is measured as the number of cells in the 
projection of the loop along this axis. 

The size of the volume, which undergoes the super- 
fluid transition in the neutron absorption event, is about 
50 /im, while the characteristic inter- vortex distance in 
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FIG. 29. Scaling properties of vortex loops in the initial net- 
work: density n (top) and average spatial size T> (middle), 
plotted as a function of the length L of the loops. Below their 
algebraic area \S\ is shown (bottom), as a function of the 2- 
dimensional spatial size T>2 of the loop. The calculations are 
performed on a lattice of 100 x 100 x 100 cells. The solid lines 
are fits to scaling laws obeying equations of the form of (31), 
(32), and (33), respectively. The inset in each plot shows the 
dependence of the scaling parameters on the size of the cubic 
neutron "bubble". 



the initial network is of the order of 1 /an. Calculations 
have been performed for several "bubble" volumes: start- 
ing from 6 x 6 x 6 up to 200 x 200 x 200. For each bubble 



size the loop distributions obtained from a large number 
of (up to 1000) initial distributions of random phases are 
averaged. The resulting distributions of n, T>, and S are 
shown in Fig. 29. One can see that despite differences 
in the boundary conditions the Vachaspati-Vilenkin re- 
lations (31) and (32) hold in the case of vortices in su- 
pcrfluid helium, but the exponents and prefactors are 
slightly different and depend on the size of the bubble. 

For the algebraic area S of a loop as a function of the 
corresponding 2-dimensional diameter T>2 the additional 
scaling law is 

-C 



\S\ = BT>2 



(33) 



where C ~ 0. Here T>2 is the average of the straight- 
line dimensions of a loop in y and z-directions. Thus 
the oriented area of a tangled loop is proportional to the 
area of a circle of the same straight-line size. The scaling 
relation (33), as well as (31) and (32), are of the form 
expected for a Brownian particle, for which the square of 
the average displacement on the i-th step is proportional 
to i and therefore the mean value of the square of the 
oriented area is given by 
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These scaling relations were used in the analytic model of 
vortex-loop extraction from the cooling neutron bubble 
in the bias flow [ie. Eqs. (22) - (24)]. 



2.14.2. Network evolution under scaling assumptions 

After its formation, the vortex network evolves un- 
der the influence of the inter-vortex interactions and the 
normal-superfluid counterflow v. Its characteristic length 
scale £(i) increases with time. Vortices start to escape 
from the bubble when the energy gain due to the exter- 
nal counterflow becomes larger than the energy of the 
superflow associated with the vortex itself. This corre- 
sponds to the critical value of £, expressed by Eqs. (11) 
or (26). 

Tangled vortex flow in superfluid He has been stud- 
ied numerically by Schwarz (1978, 1985, 1988), Samuels 
(1992), Aarts and de Waele (1994), Nemirovskii and Fisz- 
don (1994), Barcnghi et al. (1997), Tsubota and Yoneda 
(1995), and others. A large number of calculations has 
been devoted to the evolution of the initial network of 
cosmic strings (see reviews by Hindmarsh and Kibble 
1995; Bray 1994) and also of linear defects in liquid crys- 
tals (Toyoki 1994; Zapotocky et al. 1995). In all three 
cases the initial state is quite similar, but the equations 
controlling the evolution are different. The common fea- 
ture is that the interaction between the loops leads to re- 
connections when the loops cross each other. It has been 
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shown by Kagan and Svistunov (1994) that the scaling 
relations remain valid in a random network if the vortices 
are allowed to reconnect when they cross each other, but 
all other interactions are neglected for simplicity. In most 
simulation work, both in the case of cosmic strings and 
liquid crystals, it has been found that the scaling rela- 
tions are preserved during the evolution, even if inter- 
vortex interactions are included. 

To calculate the escape rate from the network, two 
crude assumptions are made, which are essentially the 
same as in the analytic treatment in Sec. 2.12.2: 1) We 
assume that the scaling laws remain valid until £ grows 
comparable in size to the critical value in Eq. (26). 2) 
At this point the influence of the external counterflow 
becomes suddenly so significant that all sufficiently large 
loops immediately escape by expanding to a rectilinear 
vortex line. In the numerical simulations the state before 
escape is modelled by the same method as was used to 
construct the initial state. It is assumed that not only 
the scaling relations but also the statistical properties of 
the vortex tangle remain the same during the evolution 
as in the "initial" state with the characteristic length £. 

At the moment of escape, £ ~ r oc v cn /v = l/x, and 
£ is used as the size of a cell in the Vachaspati-Vilenkin 
method. For integer values of x the bubble is represented 
by a grid with x x x x x vertices and a random phase is 
assigned to each vertex. To satisfy the boundary condi- 
tion, this grid is surrounded by a shell of vertices with 
fixed zero phase, representing the uniform bulk superfluid 
outside the heated bubble. Thus the whole grid contains 
[x + 1) X (x + 1) X (x + 1) cells. Such a correspondence 
between the size of the grid and the velocity is not too 
artificial even at small values of x: for example, at the 
critical velocity v = v cn one gets a grid with (x+ l) 3 = 8 
cells and 3x3x3 vertices, but the phase can be non-zero 
only at one vertex in the middle of the grid. In this case 
no vortices can appear in agreement with the definition 
of v cn - 

By counting the vortex lines produced from a number 
of random-phase distributions it is possible to calculate 
the probability distribution for the number of loops es- 
caping per absorption event. For these calculations it 
is assumed that each loop, which survives until the mo- 
ment when £ ~ r , will form an observable vortex line: 
In the case of a tangled loop the probability is high that 
at least one arc is oriented favorably with respect to the 
counterflow and will be extracted. 

In Fig. 30 the calculated probability distribution is 
plotted for observed neutron absorption events at dif- 
ferent values of x (x = 2, 3, and 4) expressed as the frac- 
tion of those neutron absorption events which produce a 
given number of vortices, normalized by the number of all 
events which give rise to at least one vortex. These results 
are in remarkably good agreement with the experimen- 
tal data (without any fitting parameters). However, the 
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FIG. 30. Distribution of the number of loops escaping per ab- 
sorption event, normalized to the total number of absorption 
events and given at three different values of the external bias 
field, v/vcn- The solid lines represent the simulations while 
the bars denote the experiment. 

agreement with the fraction of "zero" events, i.e. absorp- 
tion events which produce no vortices at all, is poorer, 
underestimating their number, especially at x = 2. Ex- 
perimentally the "zero" events can be extracted from the 
measured data by comparing the event rate at a given 
rotation velocity with the saturated event rate at the 
highest velocity. 

One reason for this discrepancy is that at low Rb/r 
ratio (ie. at low velocities) loops with a radius of cur- 
vature ~ r , which are capable of escape into the bulk 
superfluid, have a rather large probability to be oriented 
in such a way that they do not have sufficiently large seg- 
ments oriented favorably with respect to the counterflow, 
owing to space constraints. Hence they contract and do 
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not give a contribution to the observed signal. Taking 
this into account one could develop more elaborate tech- 
niques for counting vortices than the simple "count all" 
method described above. In this case better agreement 
with all experimental data could be achieved. Particu- 
larly, if only vortices with positive algebraic area with 
respect to the countcrflow (S = §y dz) are counted, then 
good agreement with the number of "zero" events would 
be obtained. However, such refinements do not change 
the results essentially, compared to uncertainties of the 
model. 



2.14.3. Direct simulation of network evolution 

There is no direct evidence for the validity of the 
assumptions, which were made in the previous section 
about the evolution of the network, ie. whether the net- 
work remains self-similar (or scale invariant) and the scal- 
ing laws Eqs. (31) - (33) can be applied during its later 
evolution. Preliminary calculations have been carried 
out using the techniques developed by Schwarz for the 
simulation of superfluid turbulence in 4 He. These meth- 
ods give good numerical agreement with experimental 
data in 4 He (Schwarz 1978, 1985, 1988) and have been 
used extensively by Schwarz and many others (Aarts and 
de Waele 1994; Nemirovskii and Fiszdon 1994; Samuels 
1992; Barenghi et al. 1997; Tsubota and Yoneda 1995). 

In such calculations vortices are generally considered to 
be 1-dimcnsional objects without internal structure. In 
4 He the diameter of the core is much smaller than other 
characteristic lengths, foremost the average radius of cur- 
vature for the loops or the inter- vortex distance, and even 
l n (£/£) can be treated as a large parameter. This is not 
the case in 3 He-B due to the large coherence length, es- 
pecially in the early stages of the evolution when the vor- 
tex density is largest. However, as mentioned in Sec. 2.3, 
the characteristic inter-vortex distance and the radius of 
curvature may still be several times larger than the co- 
herence length and the diameter of the vortex core. For 
now, we shall continue considering the vortices as linear 
objects. 

There are several forces acting on a vortex line in a 
superfluid (Donnelly 1991). The Magnus force appears 
in the presence of superflow, 



flv-I = p a KVs' x (v L - V s l), 



(34) 



where s is the radius vector of a point on the vortex line 
and the prime denotes the derivative with respect to the 
length of the line (ie. s' is a unit vector tangent to the 
vortex line at s) , Vl = s is the local velocity of the vortex 
line, and v s i is the local superfluid velocity at this point. 
For singular vortices in 3 Hc-B the number of circulation 
quanta is v = 1. The local superfluid velocity v s i is a 



sum of the superfluid velocity v s far from the network 
and the velocity induced by all the vortices in the tangle: 



v s i(r) = v s + 



(s — r) x ds 



Air /, all 
loops 



(35) 



The lordanskii force arises from the Aharonov-Bohm 
scattering of quasiparticles from the velocity field of the 
vortex, 



ftordanskii = Pn^s' X (v L - V„) , 



(36) 



where v n is the velocity of the normal component (ie. the 
heat bath of the fermionic quasiparticles). The Kopnin 
or spectral flow force has the same form, but originates 
from the spectral flow of the quasiparticle levels in the 
vortex core: 



f s f = m 3 C(T)nvs' x (vl - 



(37) 



where the temperature dependent parameter C(T) deter- 
mines the spectral flow in the core. All these three forces 
are of topological origin: they act in the transverse di- 
rection and are thus nondissipative. They are discussed 
in more details in Sec. 3. 

In contrast, the nontopological friction force ff r acts in 
the longitudinal direction, 



ffr = -d\\p s KVs' x [s' x (v„ - V L )] 



(38) 



Here the factor p s nv is the same as in the Magnus force. 

Neglecting the vortex mass, we may write the force 
balance equation for the vortex element: 



flvl + ftordanskii + fsf + ffr = . 



(39) 



It is convenient to rewrite the balance of forces in the 
following form: 

s' x [(v L - v s i) - a'(v n - v sl )] =as'x[s'x (v n - v sl )] , 

(40) 

where a and a' are the dimcnsionless mutual friction 
coefficients. The a parameters are actually the exper- 
imentally determined mutual friction quantities (see eg. 
Bevan et al. 1997b). The inverse coefficients are the d 
parameters: 



a + i(l-a') = 



1 



d|| - i(l - d_i) 



(41) 



where i = y/— 1. The transverse mutual friction parame- 
ter d± is expressed in terms of three temperature depen- 
dent functions which determine the Magnus, lordanskii 
and spectral-flow forces: 



d ± (T) = (m 3 C(T) - p n (T))/p s (T) . 



(42) 
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Eq. (40) is complicated because of the term (35) which 
contains the integral over the vortex network. To solve 
Eq. (40), one may follow Schwarz and neglect the influ- 
ence of all other vortex segments on v s i except the one 
containing the point of interest (the so called local self- 
induced approximation). In this case 

v sl w v s + /3s' x s" , where = — In ^ . (43) 

4-7T t, 

The leading correction to this simplification comes from 
the nearest neglected vortex segments and can be esti- 
mated to be of order £/[61n(£/£)], where b is the average 
inter-vortex distance. For a network in 3 He-B, £ ~ b and 
ln(£/£) ~ 1 in the early stages of the evolution and thus 
at this point the approximation is rather crude. How- 
ever, the main effect in the initial stages arises from the 
reconnection of vortex lines which happen to cross each 
other and this effect is taken into account below. 

A second simplification close to T c comes from a' <C a. 
We may rewrite the equation of vortex motion in the 
form: 

v L = v s + /3s' x s" + as' x [(v n - v s ) - /3s' x s"]. (44) 

In the rotating reference frame the normal component is 
at rest, v n = 0, and one has v s = v, where v = v n — v s 
is the counterflow velocity. We shall be comparing to 
experimental results measured at T = 0.96 T c and P = 
18 bar, where the most detailed data were collected. In 
these conditions ««10 (Bevan et al. 1995) and in the 
rotating frame Eq. (44) can be simplified further to the 
form 

s = a/3s" + as' x v , (45) 

The first term on the right hand side of Eq. (45) causes 
a loop to shrink while the second represents growth or 
shrinking, depending upon the orientation of the loop 
with respect to the counterflow. 

In the numerical calculation vortices are considered as 
lines in the same 3-dimcnsional lattice as before. The 
temporal and spatial coordinates are discrete and there- 
fore the network evolves in discrete steps: During each 
step in time a vertex of a loop can jump to one of the ad- 
jacent vertices of the lattice. The probability of a jump 
in any direction is proportional to the component of s 
in Eq. (45) in this direction, so that the average velocity 
equals s. Other more elaborate representations of vortex 
loops, based on splines for instance, have been used to 
model the dynamics of smooth vortex lines in the con- 
tinuous space. However, for the preliminary study of the 
network evolution the simplest lattice representation was 
used. The simple case of one circular loop in the coun- 
terflow, where the result is known analytically, can be 
checked separately and it is found to give correct results. 



It is well known that the interactions between neigh- 
boring vortex segments in the tangle play an important 
role in the evolution of a dense network. These interac- 
tions lead to the reconnection of loops which cross each 
other. The reconnections are the main source for the 
rarcfication of the network at the initial stages of evolu- 
tion, which process through formation and decay of small 
loops. The reconnection process has been studied numer- 
ically (Schwarz 1985; Koplik and Levine 1993) and was 
found to occur with a probability close to unity for vor- 
tex lines approaching each other. In the simulations the 
vertices are connected in such a way that closed loops are 
obtained after each time step during the evolution. If a 
cell is pierced by two vortices, the two segments coming 
in and the two segments going out are connected ran- 
domly. This corresponds to a reconnection probability 
of 0.5. It is reasonable to suppose that this convention 
leads to a slowing down in the evolution, in particular in 
the initial phase, but does not produce gross qualitative 
changes. 

The vortex tangle is initially produced by the proce- 
dure described in Sec. 2.14.2, in a lattice 40 x 40 x 40, 
and its evolution is followed until all vortices have dis- 
appeared or loops have expanded and formed circular 
planar rings far from their location of origin. These large 
rings are counted, since it is fair to assume that each one 
of them will produce a detectable rectilinear vortex line 
in the rotating container. The calculation is repeated 
for a large number of random initial distributions of the 
phase at any given value of v. The results are averaged 
to obtain the average number of vortices N\,, which are 
produced from one simulated bubble as a function of v, 
and are compared to experimental data in Fig. 31. It is 
also possible to study the validity of the scaling law (31) 
at different stages of the evolution in the course of these 
simulations. In the absence of counterflow [v = 0) it is 
found that the relation is valid for large loops, L > 4£, 
at least during much of the early evolution (in the late 
stages, when n approaches zero, the statistical noise ex- 
ceeds n). 

Fig. 31 illustrates the simulation results. No vortices 
are obtained at low counterflow velocity v < v CD , but 
when v > v cni their number N\> starts to increase rapidly 
with v. The value of the threshold velocity v cn corre- 
sponds to the situation when the largest radius of curva- 
ture r (i> C n) becomes equal to the diameter of the initial 
volume of the vortex network (ie. the diameter of the 
heated bubble). It then becomes possible for a loop to 
escape into the bulk superfluid if it consists of at least one 
arc with sufficiently large radius of curvature > r . The 
same calculations were also performed for a tangle with 
an initial volume of 20 x 20 x 20, but no differences were 
found in the dependence Nb(v/v cn ). This is additional 
circumstantial evidence for the fact that the scaling rela- 
tions approximately survive during the evolution of the 
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FIG. 31. The number of vortices Nh escaped from the bubble 
heated by one absorbed neutron as a function of v/v cn : (■) 
simulation of the network evolution in the local self-induced 
approximation (Eq. (43)) and (a) by including in addition 
approximately the polarization of the tangle by the superflow 
(Eq. (46)). Both calculations were performed on a 40 x 40 x 
40 lattice, but the results remain unchanged if the lattice is 
reduced in size to 20 x 20 x 20. The experimental data (o) 
represent measurements at P = 18 bar and T — 0.96 T c 
(Ruutu et al. 1998a). 



network also in the counterflow until the moment of es- 
cape and that the network can be approximated as scale 
invariant. 

The experimental data (denoted with (o) in Fig. 31) 
fit the universal dependence Nb(v/v cn ) — j[{v/v cn ) 3 — 1], 
where 7 equals 7 in Eq. (8) divided by the neutron flux. 
The numerical results (■) lie higher (note that no fitting 
parameters are involved) and do not display a cubic de- 
pendence in the experimental range 1 < v/v cn < 4. The 
differences could be explained by the approximations in 
the expression for v s i, which neglect significant contri- 
butions from the inter-vortex interactions. The polariza- 
tion by the external counterflow causes the loops with 
unfavourable orientation to contract and loops with the 
proper winding direction and orientation to grow. How- 
ever, as the polarization evolves, it also screens the vortex 
tangle from the external counterflow. To check whether 
the polarization has a significant effect on the results, 
the calculations were repeated by taking the screening 
approximately into account in the expression for v s i: 



v s i = v s + f3s' x s 



! " + ^ /other 
loops 



(s — r) x ds 



(46) 



Here the contribution from each loop to the superflow is 
averaged over the volume. The results (a) show that this 



effect is significant and should be taken into account more 
accurately, to reproduce the network evolution correctly. 

To summarize we note that the preliminary simulation 
work shows that it is possible to obtain numerical agree- 
ment between the KZ mechanism and the rotating exper- 
iments, if one assumes that the scaling relations of the 
initial vortex tangle are roughly obeyed also during its 
later evolution in the presence of the external superflow. 
This assumption should still be checked with sufficiently 
accurate simulation, with full account of the non-local 
inter-vortex interactions. Future simulation of transient 
networks should then answer the question how much in- 
formation about the initial state of the vortex tangle can 
be retrieved from the rotating experiment, where only 
the final stationary state result is measured. 



2.15. Superfluid transition as a moving phase front 

In any real laboratory system a rapid phase transition 
becomes inhomogeneous: The transition will be driven 
by a strong gradient in one or several of the externally 
controlled variables. In the case of the superfluid transi- 
tion in the aftermath of a neutron absorption event it is 
a steep and rapidly relaxing thermal gradient. In this sit- 
uation even a second order transition becomes spatially 
inhomogeneous and is localized into a phase front. In the 
limit of very fast transitions the order-parameter relax- 
ation slows down the propagation of the superfluid phase 
and the thermal front, where the temperature drops be- 
low T c , may escape ahead. This means that the mov- 
ing phase boundary breaks down into a leading thermal 
front and a trailing superfluid interface. The width of 
the region between these two zones is determined by 
the relative velocities of thermal and superfluid relax- 
ations. Thermal fluctuations in this region are amplified 
in the transition process and are carried over as order- 
parameter inhomogencity into the new symmetry-broken 
phase. This is the central claim of the KZ model. Be- 
low we shall briefly discuss the influence of the thermal 
gradient on defect formation, as analyzed by Kibble and 
Volovik (1997) and Kopnin and Thuneberg (1999). 



2.15.1. Neutron absorption and heating 

The decay products from the neutron absorption re- 
action generate ionization tracks which can be estimated 
from a standard calculation of stopping power (Meyer 
and Sloan 1997). This leads to a cigar-shaped volume 
of ionized particles, with the largest concentration at the 
end points of the two tracks. The probabilities and re- 
laxation times of the different recombination channels for 
the ionized charge are not well known in liquid 3 He. Also 
the thermalization of the reaction energy may not pro- 
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duce a heated region which preserves the shape of the 
original volume with the ionized charge. 

Initially the recombination processes are expected to 
lead to particles with large kinetic energies in the eV 
range, which are well outside the thermal distribution 
and for which the recoil velocities become more and more 
randomly oriented. Energetic particles suffer collisions 
with their nearest neighbors and the mean free path in- 
creases only slowly for atoms participating in these col- 
lisions, until all particles are slowed down and become 
thermalized to the ambient conditions (Bunkov and Tim- 
ofeevskaya 1998a, b). This means that the reaction energy 
remains initially localized. In the calorimetric experi- 
ments at the lowest temperatures (Bauerle et al. 1996, 
1998a) the thermal mean free path exceeds the container 
dimensions. Nevertheless, the energy is probably not im- 
mediately dispersed into the entire container volume, but 
remains centralized within a bubble of limited size during 
cooling through T c , when the vortex network is formed. 
This is the conclusion to be drawn from the comparison 
between experiment and the KZ mechanism in Table I, 
where it is assumed that the thermal diffusion mean free 
path is the same as in the normal fluid at T c . 

In an inhomogeneous initial state with large thermal 
gradients the second order phase transition is turned into 
one where a normal-to-superfluid phase front with finite 
width sweeps through the heated bubble. If the velocity 
of the phase front, vt ~ R^/tq ~ 6 m/s, is sufficiently 
high, comparable to a critical value «tc ~ v-p (to/tq)j, 
then the KZ mechanism is again expected to dominate, 
similarly as in the homogeneous case (Kibble and Volovik 
1997). 

In contrast, if the majority of the kinetic energy is 
assumed to be thermalized by quasiparticles with en- 
ergies comparable to the high-energy tail of the ther- 
mal Maxwcllian velocity distribution, then the mean free 
paths are long, the volume heated above the ambient be- 
comes large, and its temperature distribution may even 
become nonmonotonic like in a "Baked Alaska", as has 
been described by Leggett (1992). The Baked Alaska sce- 
nario is also popular in high energy physics, where the 
formation of the false vacuum with a chiral condensate af- 
ter a hadron-hadron collision is considered (Bjorkcn 1997; 
Amclino-Camelia et al. 1997). In both cases a rather thin 
shell of radiated high energy particles expands, with the 
speed of light in a relativistic system and with the Fermi 
velocity v-p in 3 Hc, leaving behind a region at reduced 
temperature. Since this interior region is separated from 
the exterior vacuum by the hot shell, the cooldown into 
the broken symmetry state in the center is not biased by 
the external order parameter state. The Baked Alaska 
mechanism thus can solve the problem of the neutron- 
mediated formation of B phase from supercooled bulk 
A liquid (Leggett 1992), while in high energy physics it 
opens the possibility for the formation of a bubble of chi- 



ral condensate in a high energy collision (Bjorkcn 1997; 
Amelino-Camelia et al. 1997). 

In such conditions, when the quasiparticle mean free 
path exceeds or is comparable to the dimensions of the 
heated bubble, temperature is not a useful quantity for 
the description of its cooling. Most of the analysis of the 
previous sections is applicable only if we assume that the 
reaction energy remains reasonably well localized while 
the hot bubble cools through T c . In this case there is no 
Bakcd-Alaska effect: No hot shell will be formed which 
would separate the interior from the exterior region. In 
this situation the exterior region could be imagined to fix 
the phase in the cooling bubble, while the phase front is 
moving inward, suppressing the formation of new order- 
parameter states, which arc different from that in the 
bulk supcrfluid outside, and in the same manner sup- 
pressing the formation of vortices. However, it appears 
that there also exists another alternative: The influence 
of the exterior region may not be decisive if the phase 
transition front moves sufficiently rapidly. Which of these 
alternatives is realized in a particular situation is still 
very much a subject of discussion. 

For the interpretation of the measurements in neutron 
irradiation, a sophisticated understanding of the shape 
and size of the constant temperature contours within the 
heated bubble is not vitally necessary. In the final results 
the bubble size does not enter, since the data can be nor- 
malized in terms of the measured threshold velocity v cn . 
Its measurement provides an estimate of the circumfer- 
ence of the bubble, since the largest possible vortex ring 
has to be comparable in size to the neutron bubble. The 
diameter of this ring is 1 - 2 orders of magnitude larger 
than the expected average inter- vortex distance £ v in the 
initial random network which is created by the KZ mech- 
anism. 

2.15.2. Thermal gradient and velocity of phase front 

For a rough understanding of the superfluid transition 
in a temperature gradient let us consider the time de- 
pendent Ginzburg-Landau (TDGL) equation for a one- 
component order parameter W = A/Aq: 

T ^=(l-^)*-*|*| 2 +e 2 V 2 * , (47) 

where the notations for r cx 1/ A and £ are the same as 
in Sec. 2.3. This equation is the so-called over-damped 
limit of the more general TDGL equation which has a 
time derivative of second order. The over-damped limit 
has been used in numerical simulations (Laguna and 
Zurek 1997; Antunes et al. 1999; Aranson et al. 1999) 
and analytical estimations (Dziarmaga 1998, 1999) of the 
density of topological defects in a homogeneous quench. 
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An extension of the above equation to superconductiv- 
ity was used by Ibaceta and Calzetta (1999) to study 
the formation of defects after a homogeneous quench in 
a 2-dimensional type II superconductor. 

If the quench occurs homogeneously in the whole space 
r, the temperature depends only on one parameter, the 
quench time tq: 



T(t) 



1 — I 



(48) 



In the presence of a temperature gradient, say, along x, 
a new parameter appears, which together with tq char- 
acterizes the evolution of the temperature: 



TV- , l /)=(l-^^ , |T, 



(49) 



Here is the velocity of the temperature front which is 
related to the temperature gradient (Kibble and Volovik 
1997) 



VTTQ 



(50) 



With this new parameter vt we may classify transi- 
tions to belong to one of two limiting regimes: slow or 
fast propagation of the transition front. At slow veloci- 
ties, vt — > 0, the order parameter moves in step with the 
temperature front and 



|*(*,t)| 2 = 1- 



T(x - v T t) 



e(l-T(x-v T t)/T c ] 



(51) 

Here 8 is the step function. In this case phase coherence 
is preserved behind the transition front and no defect 
formation is possible. 

The opposite limit of large velocities, vt — > oo, is more 
interesting. Here the phase transition front lags behind 
the temperature front (Kopnin and Thuneberg 1999). In 
the space between these two boundaries the temperature 
is already below T c , but the phase transition did not yet 
happen, and the order parameter has not yet formed, 
\I> = 0. This situation is unstable towards the formation 
of blobs of the new phase with ^ ^ 0. This occurs inde- 
pendently in different regions of space, leading to vortex 
formation via the KZ mechanism. At a given point r the 
development of the instability can be found from the lin- 
earized TDGL equation, since during the initial growth 
of the order parameter \& from zero the cubic term can 
be neglected: 



d# t 

T — = — * 



dt 



(52) 



This gives an exponentially growing order parameter, 
which starts from some seed ^fluc, caused by fluctua- 
tions: 



tf(r,i) = # fluc (r) exp 



2tqt 



(53) 



Because of exponential growth, even if the seed is small, 
the modulus of the order parameter reaches its equi- 
librium value 

^Zurck- 



cq 



= yjl - T/T c after the Zurek time 



iZurek = V T Q T 



(54) 



This occurs independently in different regions of space 
and thus their order-parameter phases are not corre- 
lated. The spatial correlation is lost over distances ex- 
ceeding £ v when the gradient term in Eq. (47) becomes 
comparable to the other terms at t = tzurck- Equating 
the gradient term ^V 2 f <~ (Co/Cv)^ to, say, the term 
Todfy / dt\ tzurck — ■\/to/tq5', one obtains the characteris- 
tic Zurek length which determines the initial length scale 
of defects, 



Cv = Co (tq/t ) 



1/4 



(55) 



which is the same as in Eq. (3). 

We can estimate the lower limit for the characteristic 
value of the fluctuations ^Pfluc = Afl uc /Ao, which serve as 
a seed for vortex formation. If there are no other sources 
of fluctuations, caused by external noise for example, the 
initial seed is provided by thermal fluctuations of the 
order parameter in the volume The energy of such 
fluctuations is ^A^^Np/Ep, where Ep ~ k^Tp is the 
Fermi energy and Np the fermionic density of states in 
normal Fermi liquid. Equating this energy to tempera- 
ture T w T c one obtains the magnitude of the thermal 
fluctuations 



I* 



cq 



1/8 



k B T c 



(56) 



The same small parameter ksTc/Ep ~ a/£o ~ 10~ 3 — 
10~ 2 enters here, which is responsible for the extremely 
narrow temperature region of the critical fluctuations in 
3 He near T c . But in our case it only slightly increases the 
Zurek time by the factor •^/ln[|^ , cq |/| 1 J r fiuc|] and thus does 
not influence vortex formation in a homogeneous quench 
or in an inhomogeneous quench at large velocities of the 
temperature front. 

Clearly there must exist a characteristic velocity v t c of 
the propagating temperature front, which separates the 
fast and slow regimes, or correspondingly transitions with 
and without defect formation. The criterion for defect 
formation is that the Zurek time tzurok = ^/ T Q T o should 
be shorter than the time t sw in which the phase transition 
front sweeps through the space between the two bound- 
aries. The latter time is t sw = xq/vt, where xo is the 
lag - the distance between the temperature front T = T c 
and the region where superfluid coherence starts (order 
parameter front). If izurok < t sw , instabilities have time 
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to develop. If izurck > t sw , both fronts move coherently 
and the phase is continuous. Let us consider the latter 
case, with "laminar" motion, and find how the lag Xq de- 
pends on vt- From the equation izm-ek = £o(vT c )/'fT c we 
find the criterion for the threshold where laminar motion 
becomes unstable and defect formation starts. 

In steady laminar motion the order parameter depends 
on x — v^t. Introducing a dimensionless variable z and a 
dimcnsionless parameter ft, 



z=(x- v T t)(v T T Q $) 



-1/3 



/ V T T \ 4/3 / TQ 
I Co J U 



the linearized TDGL equation becomes 

, 



dz 2 dz 



1/3 



(57) 



(58) 



or 



= const • e hz/2 x(z) 



fx 

dz 2 



- z 



x = o 



(59) 



This means that \& is an Airy function, x( z ~ z o), centered 
at z — z = — h 2 /A and attenuated by the exponential 
factor e~ hz / 2 . 

When ft > 1, it follows from Eq. (59) that *(z) quickly 
vanishes as z increases above — ft 2 /4. Thus there is a 
supercooled region — h 2 /A < z < 0, where T < T c , but 
the order parameter is not yet formed: the solution is 
essentially \J f = 0. The lag between the order parameter 
and temperature fronts is \z \ = h 2 /A or in conventional 
units 



xq = — 



1 VtTQ t 

4 e 



(60) 



Setting izurck = x o(vt c )/vt c one can estimate the value 
for the velocity of the temperature front where laminar 
propagation becomes unstable: 



VTC — _ 

T \TQ 



1/4 



(61) 



This result corresponds to ft ~ 1 and is in agreement 
with that obtained from scaling arguments by Kibble and 
Volovik (1997). Here an exact numerical value cannot be 
offered for the threshold where laminar motion ends, but 
this can in principle be done using the complete TDGL 
equation (Kopnin and Thuneberg 1999). For the neutron 
bubble we might take vt ~ Rb/Tq, which gives wt ~ 10 
m/ s. This value provides also an estimate for the limiting 
velocity v^ c ■ These considerations suggest that the ther- 
mal gradient should be sufficiently steep in the neutron 
bubble such that defect formation is to be expected. 




FIG. 32. Normal-to-superfluid transition, when described as 
a change-over from an infinite to a finite vortex network. 
(Top) The disordered phase (T > T c ) is an infinite network 
of the defects in the ordered phase - the vacuum of tangled 
closed loops. The coherence length £ is the mean distance be- 
tween the elements in this network. (Bottom) In the ordered 
phase (T < T c ) the defects also form closed loops, but the 
coherence length is now the upper limit in the distribution of 
loop sizes. 



2.16. Quench of infinite vortex tangle 

2.16.1. Vorticity on microscopic and macroscopic scales 

Onsager (1949) was the first to interpret the A- 
transition in liquid 4 He from the superfluid to the normal 
state in terms of quantized vortices: When the concentra- 
tion of the thermally activated quantized vortices reaches 
the point where they form a connected tangle through- 
out the liquid (and their line tension vanishes), the liquid 
goes normal (Fig. 32). This proliferation of vortices to 
an infinite network destroys superfluidity, since the phase 
slippage processes caused by the back reaction from the 
tangle to the superfluid current lead to the decay of this 
current. 

In Ginzburg-Landau theory vortices are identified as 
lines of zeroes in the scalar complex order parameter 

= l^le**, around which the phase winding number v is 
nonzero. Above the transition, in the symmetric phase, 
thermal fluctuations of the order parameter give rise to 
an infinite network of zeroes - to vortices which exist on 
the microscopic scale but are absent on the macroscopic 
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scale (Kleinert 1989). This is another way of describing 
the fact that there is no long range order in the symmet- 
ric phase. The properties of such microscopic vortices 
- topologically nontrivial zeroes - have been followed 
across the thermodynamic phase transition in numerical 
investigations (Antunes and Bettencourt 1998; Antunes 
et al. 1998; Rajantie et al. 1998; Kajantie et al. 1998; Ra- 
jantie 1998). The renormalization group description of 
the phase transition, based on the Ginzburg-Landau free 
energy functional, also contains the microscopic quan- 
tized vortices, but in an implicit form of zeroes (see e.g. 
Guillou and Zinn- Justin 1980 and Albert 1982). 

Some attempts have been made to reformulate the 3- 
dimensional phase transition in terms of this vortex pic- 
ture in the renormalization group approach (Williams 
1993a, b, 1999; Chattopadhyay et al. 1993), in a man- 
ner similar to the 2-dimcnsional Berczinskii-Kosterlitz- 
Thouless transition (Nelson and Kosterlitz 1977), i.e. by 
avoiding consideration of the Ginzburg-Landau free en- 
ergy functional altogether. It is possible, though it has 
not been proven, that the Ginzburg-Landau model and 
the vortex model belong to the same universality class 
and thus give the same critical exponents for the heat 
capacity and superfluid density. Both ground and space- 
based measurements of the critical exponents (Lipa et al. 
2003, 1996; Goldner and Ahlers 1992) have generally 
been consistent with renormalization-group estimations 
(Strosser and Dohm 2003; Campostrini et al. 2001; Klein- 
ert and den Bossche 2001). 

In an equilibrium phase transition the infinite network 
disappears below T c , but in a non-equilibrium phase tran- 
sition the tangle of microscopic vortices persists even 
into the ordered phase, due to critical slowing down. 
These vortices finally transform to conventional macro- 
scopic vortices, when the latter become well defined. In 
this language the Kibble-Zurek mechanism corresponds 
to a quench of the infinite vortex network across the 
non-equilibrium transition from the normal to the su- 
perfluid phase (Yates and Zurek 1998). What is impor- 
tant for us here, is the scaling law for the distribution 
of vortices. According to numerical simulations (An- 
tunes and Bettencourt 1998; Antunes et al. 1998) and 
the phenomenological vortex model (Williams 1993a, b, 
1999; Chattopadhyay et al. 1993) the scaling exponent 
5, which characterizes the distribution of vortex loops 
in Eq. 32), is close to the value 5 = 2/(d + 2) = 0.4 
obtained using the concept of Flory calculations for self- 
avoiding polymers. It is different from the value given 
by the Vachaspati-Vilenkin model for a Brownian-walk 
network. On the other hand, this model was applied to 
the Ginzburg region of critical fluctuations in 4 He, which 
might not be valid for 3 He. 



2.16.2. Scaling in equilibrium phase transitions 

Important differences exist between the phase transi- 
tions in the 4 He and 3 He liquids. These also are involved 
in a number of other phenomena. Let us therefore recall 
some of these differences. The temperature region of the 
critical fluctuations, where the simple Ginzburg-Landau 
theory does not work and one must introduce the thermal 
renormalization of the Ginzburg-Landau functional, can 
be derived from the following simplified considerations. 
Let us estimate the length Lt of the thermal vortex loop, 
ie. a loop whose energy is comparable to temperature: 



= k B T w k B T c 



(62) 



In the broken symmetry phase far below T c , this length 
is less than the coherence length £(T), which means 
that there are no real vortices with the energy of or- 
der k B T, and real vortices with higher energy are expo- 
nentially suppressed. When the temperature increases, 
one approaches the point at which the length Lt be- 
comes comparable to the coherence length £. This is the 
Ginzburg temperature Tq- u determined by the condition 
L T (T = T Gi )~Z(T = T Gi ): 



or 



1 



Ps (T Gi )£(T Gi ) ~ Mk 



T G i A ( ^bT c \ 2 fT r s 



(63) 



(64) 



Here as before Tp ~ ti 2 /ma 2 <~ 1 K is the degener- 
acy temperature of the quantum fluid, with a being the 
inter-atomic spacing. The region of critical fluctuations 
- the Ginzburg region T G i < T < T c - is broad for 
4 He, where T c ~ Tp, and extremely small for 3 He, where 
T c < 10~ 2 Tp. 

In the region of Ginzburg fluctuations the scaling expo- 
nents for the thermodynamic quantities, such as £(T) and 
p s (T), are different from those in the Ginzburg-Landau 
region, T c - T Gi < T c - T < T c : 

~ i 1 ~ y) ' Tc ~ Tci < Tc ~ T <<: Tc ' (65) 
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■! - - ) , T Gi < T < T c , (66) 



T c — T Gl < T c - T < T c , (67) 



1-f ) - T^i < ' T { . ((,8) 



There are two relations, the scaling hypotheses, which 
connect the critical exponents for £(T) and p s (T) in the 
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Ginzburg region, with the exponent of the heat capac- 
ity. In the Ginzburg region Tci < T < T c , the coherence 
length is determined by thermal fluctuations, or which is 
the same thing, by thermal vortices. This gives for the re- 
lation between the coherence length £(T) and superfiuid 
density p s (T) in the Ginzburg region: 

p s (T)Z(T) ~^~ P a (^j , T Gi < T < T c . (69) 

This equation gives the Joscphson scaling hypothesis: 
v = C 

Another relation between these exponents comes from 
a consideration of the free energy, which has the same 
scaling law as the kinetic energy of superflow: 



comparable to the time t which is left until the transition 
takes place. At this moment 
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(70) 



Here a is the the critical exponent for the heat capacity 
in the critical region, C(T) = -Td\F - (1 - T/T c )- a . 
Eqs. (69) and (70) give v = C = (2 - a)/3. 



2.16.3. Non-equilibrium phase transitions 

The formation of vortices during a rapid transition into 
the broken-symmetry phase is the subject of dynamic 
scaling and a poorly known area in the field of critical 
phenomena. Dynamic scaling is characterized by an ad- 
ditional set of critical exponents, which depend not only 
on the symmetry and topology of the order parameter, 
but also on the interaction of the order parameter with 
the different dynamic modes of the normal liquid. The 
question first posed by Zurek was the following: What is 
the initial density £ v of macroscopic vortices at the mo- 
ment when they become well defined? According to the 
general scaling hypothesis one has 



£v ~ £o(^To) A and w : 



1 



(71) 



where u> is the characteristic frequency of the dynamic 
process. 

In the time-dependent Ginzburg-Landau model 
(Eq. (47)) and also in its extension based on the renor- 
malization group approach, the exponent A is determined 
by the static exponents and by the exponent for the re- 
laxation time t = t (1 - T/T c )-^: 



A = 



1 + H 



(72) 



This follows from the following consideration. When we 
approach the critical temperature from the normal phase, 
at some moment tzurek the relaxation time r(t) becomes 



iZurck — t [ — ) 



(73) 



the vortex network is frozen out. After the transition 
it becomes "unfrozen" when r(t) again becomes smaller 
than t, the time after passing the transition. The initial 
distance between the vortices is determined as the co- 
herence length £ = £o(*/ t q) -1 ' at i = izurck, which gives 
Eq. (71) with A as in Eq. (72). For the conventional time- 
dependent Ginzburg-Landau model in Eq. (47) one has 
p, = 1, v = 1/2, and thus A = 1/4. In the Ginzburg 
regime, if one assumes that fi remains the same, while 
v fa 2/3, one obtains A w 1/3. 

In numerical simulation of a quench in the time- 
dependent Ginzburg-Landau model, the main problem 
becomes how to resolve between the microscopic vortices 
(ic. the zeroes in the order parameter) and real macro- 
scopic vortices. This requires some coarse-graining pro- 
cedure, which is not well established. Also speculations 
exist (Rivers 2000) that if the quench is limited to within 
the region of the Ginzburg fluctuations, then the network 
of microscopic vortices might effectively screen out the 
real macroscopic vortices, so that the density of the real 
vortices after the quench is essentially less than its Zurek 
estimate. This could explain the negative result of the 
recent pressure-quench experiments in liquid 4 He (Dodd 
et al. 1998), where the final state after the decompression 
was well within the region of critical fluctuations. 

To summarize, we note that vortices play an impor- 
tant role, not only in the ordered state, but also in the 
physics of the broken symmetry phase transitions. The 
proliferation of vortex loops with infinite size can be in- 
terpreted to destroy the superfiuid long range order above 
the phase transition. In a non-equilibrium phase transi- 
tion from the symmetric normal phase to the superfiuid 
state, the infinite vortex cluster from the normal state 
survives after the rapid quench and becomes the source 
for the remanent vorticity in the superfiuid state. This is 
an alternative picture in which one may understand the 
formation of vortices in a rapid non-equilibrium transi- 
tion. 



2.17. Implications of the quench-cooled experiments 

2.17.1. Topological-defect formation 

It is not obvious that a phenomenological model like 
the KZ mechanism, which is based on scaling arguments, 
should work at all: It describes a time-dependent phase 
transition in terms of quantities characterizing the equi- 
librium properties of the system. Numerical calculations 
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on simple quantum systems, where one studies the fluc- 
tuations in the amplitude of the system wave function 
while it is quench cooled below a second order phase tran- 
sition, have provided much evidence for the KZ model 
and appear to agree with its qualitative features. Most 
attempts of experimental verification suffer from short- 
comings. Measurements on superfluid 3 Hc arc the first to 
test the KZ model more quantitatively. In this case the 
experimental deviations from the ideal KZ setup include 
the presence of a strict boundary condition and a strong 
thermal gradient. At present time we can conclude that 
the experiments and the model are in reasonable har- 
mony, assuming that open questions from Sec. 2.14 can 
be answered satisfactorily, as seems likely. 

However, even good agreement leaves us with an in- 
teresting question: What is the microscopic basis for the 
applicability of the KZ mechanism to such experiments? 
A rapid quench through the superfluid transition is more 
amenable to microscopic analysis in the case of liquid 3 He 
than in most other systems, since the freeze out of order- 
parameter domains can be demonstrated with physically 
acceptable calculations. The consequences from this are 
exciting and the prospects for a better understanding of 
non-equilibrium phase transitions look promising. Al- 
though detailed agreement has not yet been reached be- 
tween experimental and theoretical work, nevertheless 
the effort by Aranson et al. (1999, 2001) illustrates that 
many aspects of the neutron absorption event in 3 Hc-B 
can be treated in realistic ways. 

Are there implications from such work to cosmological 
large-scale structure formation? The combined existing 
evidence from condensed matter experiments supports 
the KZ mechanism at high transition velocities. This 
says that the KZ mechanism is an important process in 
rapid non-equilibrium phase transitions and cannot be 
neglected. Whether it was effective in the Early Uni- 
verse and led to the formation of large-scale structure, 
rather than some other phenomenon like inflationary ex- 
pansion, is a separate question. It can only be answered 
by measurements and investigations on the cosmological 
scale which directly search for the evidence on topolog- 
ical defect formation or its exclusion. The fact that no 
direct traces from topological defects have been identi- 
fied, excepting the large-scale structure itself, is not yet 
sufficient proof that they can be discarded altogether. 
As measurements of the moving first-order AB interface 
in 3 He superfluids show, the interactions of the defects 
with the phase front and across it are complicated and 
no simple general rule exist on the final outcome. 

2.17.2. Phase transitions 

As discussed in Sees. 2.10.1 and 2.10.2, the KZ mech- 
anism provides an attractive explanation for the radia- 



tion induced first order transition from supercooled 3 He- 
A to 3 He-B. The nucleation of competing phases, with 
different symmetries and local minima of the energy func- 
tional, has been discussed both in superfluid 3 He (Volovik 
2003) and in the cosmological context (Linde 1990). 

Above the critical phase-transition temperature T c , liq- 
uid 3 He is in its symmetric phase: it has all the symme- 
tries allowed in non-relativistic condensed matter. The 
continuous symmetries, whose breaking are relevant for 
the topological classification of the defects in the non- 
symmetric phases of 3 He, form the symmetry group 

G = 50(3) L x 50(3) s x E/(1) N . (74) 

Here 5(7(3)l is the group of solid rotations of the coor- 
dinate space. The spin rotations of the group 5<3(3)s 
may be considered as a separate symmetry operation if 
one neglects the dipolar spin-orbit interaction. The mag- 
netic dipole interaction between the nuclear spins is tiny 
in comparison with the energies characterizing the super- 
fluid transition. The group C/(1)n is the global symmetry 
group of gauge transformations, which stems from the 
conservation of the particle number N for the 3 He atoms 
in their ground states. U(l) is an exact symmetry if one 
neglects extremely rare processes of excitations and ion- 
ization of the 3 He atoms, as well as the transformation 
of 3 He nuclei, in neutron radiation. 

Below T c the "unified GUT symmetry" 50(3) L x 
50(3) s x C/(l) N can be broken to the U(l) x U(l) sym- 
metry of the A phase or to the SO (3) symmetry of the 
B-phase, with a small energy difference between these 
two states, but separated by an extremely high energy 
barrier (~ lO 6 /^? 1 ) from each other. In this situation a 
thermally activated A^B transition becomes impossible. 
For the radiation-induced A^B transition two solutions 
have been suggested: 1) the baked Alaska configuration 
where the transition takes place inside a cool bubble iso- 
lated by a warm shell from the surrounding A-phase bath 
(Leggett 1984), and 2) the KZ mechanism, first proposed 
by Volovik (1996) and then put on quantitative ground 
by Bunkov and Timofeevskaya (1998a,b). 

This situation can be compared to that in the early 
Universe. It is believed that the 517(3) x 517(2) x 17(1) 
symmetries of the strong, weak, and electromagnetic in- 
teractions (respectively) were united at high energies (or 
at high temperatures). The underlying Grand Unifi- 
cation (GUT) symmetry (5f7(5), 5O(10), or a larger 
group) was broken at an early stage during the cool- 
ing of the Universe. Even the simplest GUT symmetry 
SU(5) can be broken in different ways: into the phase 
5/7(3) x SU(2) x [7(f), which is our world, and into 
5f7(4) x £7(1), which apparently corresponds to a higher 
energy state. In supersymmetric models both phases rep- 
resent local minima of almost equal depth, but are sepa- 
rated from each other by a high energy barrier. 

In the cosmological scenario, after the symmetry break 
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of the SU(5) GUT state, both new phases are created 
simultaneously with domain walls between them. The 
calculated probability for the creation of our world — the 
SU(3) x SU(2) x U(l) state — appears to be smaller than 
that of the false vacuum state of SU (4) x 17(1) symmetry 
at higher energy. Thus initially the state corresponding 
to our world appears to have occupied only a fraction 
of the total volume. Later bubbles of this energetically 
preferred state grew at the expense of the false vacuum 
state and finally completely expelled it. However, before 
that interfaces between the two states were created in 
those places where blobs of the two phases met. Such 
an interface is an additional topologically stable defect, 
which is formed in the transition process and interacts 
with other topological defects. 

3. VORTEX DYNAMICS AND QUANTUM FIELD 
THEORY ANALOGUES 

From the phenomcnological point of view, vortex dy- 
namics in 3 He-B is similar to 4 He-II, but the properties 
which govern the dynamics arise from the quite differ- 
ent microscopies of the p-wave Cooper-paired fermion 
superfluid. These properties turn out to have interesting 
analogies with various problems in quantum field theory 
(Volovik 2003). It is these connections which are in the 
focus of the discussion in this section. 

The fundamental starting point is the notion that the 
bosonic and fermionic excitations in superfluid 3 He are 
in many respects similar to the excitations of the ener- 
getic physical vacuum of elementary particle physics - the 
modern ether. This similarity allows us to model, with 
concepts borrowed from 3 He physics, the interactions of 
elementary particles with the evolving strings and do- 
main walls, which are formed in a rapid phase transition, 
for instance. Such processes become important after de- 
fect formation in the initial quench and give rise to the 
cosmological consequences which we are measuring to- 
day. 

The quantum physical vacuum - the former empty 
space - is in reality a richly structured and asymmet- 
ric medium. Because the new quantum ether is such 
a complicated material with many degrees of freedom, 
one can learn to analyze it by studying other materi- 
als, ie. condensed matter (Wilczek 1998). Fermi super- 
fluids, especially 3 He-A, are the best examples, which 
provide opportunities for such modelling. The most pro- 
nounced property of 3 He-A is that, in addition to the 
numerous bosonic fields (collective modes of the order 
parameter which play the part of gauge fields in elec- 
tromagnetic, weak, and strong interactions) it contains 
gaplcss fermionic quasiparticles, which are similar to the 
elementary excitations of the quantum physical vacuum 
(leptons and quarks). 



It is important that the quantum physical vacuum be- 
longs to the same class of fermionic condensed matter 
as 3 He-A: both contain topologically stable nodes in the 
energy spectrum of the fermionic excitations. As a re- 
sult both of these fermionic systems display, for example, 
the gravitational and gauge fields as collective bosonic 
modes. Other fermionic systems belong either to a class, 
which is characterized by Fermi surfaces (such as normal 
metals and the normal 3 He liquid), or to a class with a 
gap in the fermionic spectrum (such as conventional su- 
perconductors and 3 He-B). High-temperature supercon- 
ductors seem to belong to a marginal class, with topolog- 
ically unstable lines of gap nodes. Thus 3 He-A (together 
with 3 He-Ai) is the condensed matter system in which 
the properties of the physical vacuum can, in principle, 
be probed with laboratory experiments. 

3.1. Three topological forces acting on a vortex and their 
analogues 

Here we consider the experimentally observed forces 
which act on a moving vortex in superfluid 3 He. As listed 
in Sec. 2.14.3, there are 3 different topological contribu- 
tions to the total force. The more familiar Magnus force 
arises when the vortex moves with respect to the super- 
fluid vacuum. In the case of a relativistic cosmic string 
this force is absent, since the corresponding superfluid 
density of the quantum physical vacuum is zero. How- 
ever, the analog of this force appears if the cosmic string 
moves in a uniform background charge density (Davis 
and Shcllard 1989; Lee 1994). The other two forces of 
topological origin - the Iordanskii force and spectral flow 
force - also have analogues in the case of a cosmic string 
(Volovik and Vachaspati 1996; Volovik 1998; Sonin 1997; 
Wexler 1997; Shelankov 1998a,b). 

We start with the Iordanskii force (Iordanskii 1964, 
1965; Sonin 1975), which arises when the vortex moves 
with respect to the heat bath represented by the normal 
component of the liquid, or the quasiparticle excitations. 
The latter corresponds to the matter of particle physics. 
The interaction of quasiparticles with the velocity field of 
the vortex resembles the interaction of matter with the 
gravitational field induced by such a cosmic string, which 
has an angular momentum, - the so-called spinning cos- 
mic string (Mazur 1986). The spinning string induces 
a peculiar space-time metric, which leads to a difference 
in the time which a particle needs to orbit around the 
string at same speed, but in opposite directions (Harari 
and Polychronakos 1988). This gives rise to the quan- 
tum gravitational Aharonov-Bohm effect (Mazur 1986, 
1987, 1996). We discuss how the same effect leads to an 
asymmetry in the scattering of particles from the spin- 
ning string and to the Iordanskii lifting force which acts 
on the spinning string or on the moving vortex. 
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The spectral flow force, which also arises when the vor- 
tex moves with respect to the heat bath, is a direct con- 
sequence from the chiral anomaly effect. The latter vio- 
lates the conservation of fermionic charge. The anoma- 
lous generation by the moving vortex of fermionic charge 
or momentum (called "momentogenesis" ) leads to a net 
force acting on the vortex. The existence of this force 
was experimentally confirmed in the 3 He measurements 
of Bevan et al. (1997b). This phenomenon is based on 
the same physics as the anomalous generation of matter 
in particle physics and bears directly on the cosmologi- 
cal problem of baryonic asymmetry of our Universe: the 
question why the Universe contains so much more matter 
than antimatter ( "baryogenesis" ) . 

The experimental observation of the opposite effect 
to momentogenesis has been reported by Krusius et al. 
(1998): the conversion of quasiparticle momentum into a 
non-trivial order parameter configuration or "texture". 
The corresponding process in a cosmological setting 
would be the creation of a primordial magnetic field due 
to changes in the matter content. Processes, in which 
magnetic fields are generated, are very relevant to cos- 
mology since magnetic fields are ubiquitous now in the 
Universe. Our Milky Way and other galaxies, as well 
as clusters of galaxies, are observed to have a magnetic 
field whose generation is still not understood. One pos- 
sible mechanism is that a seed field was amplified by the 
complex motions associated with galaxies and clusters of 
galaxies. The seed field itself is usually assumed to be of 
cosmological origin. 

It has been noted that the two problems of cosmologi- 
cal genesis - baryo- and magnctogenesis - may be related 
to each other (Roberge 1989; Vachaspati and Field 1994, 
1995; Vachaspati 1994). More recently an even stronger 
reason for a possible connection was proposed (Joyce and 
Shaposhnikov 1997; Giovannini and Shaposhnikov 1997). 
In this same way their analogs are related in 3 He, where 
the order parameter texture is the analog of magnetic 
field, while the normal component of the superfluid rep- 
resents the matter: the moving vortex texture leads to an 
anomalous production of quasiparticles, while an excess 
of the quasiparticle momentum - the net quasiparticle 
current of the normal component - leads to the formation 
of textures. This mapping of cosmology to condensed 
matter is not simply a picture: the corresponding effects 
in the two systems are described by the same equations 
in the low-energy regime, by quantum field theory and 
the axial anomaly. 

3.2. lordanskii force 

3.2.1. Superfluid vortex vs spinning cosmic string 

To clarify the analogy between the lordanskii force and 
the Aharonov-Bohm effect, let us consider the simplest 



case of phonons propagating in the velocity field of the 
quantized vortex in the Bose superfluid 4 He-II. According 
to the Landau theory of superfluidity, the energy of a 
quasiparticle moving in the superfluid velocity field v s (r) 
is Dopplcr shifted: E(p) = e(p) + p-v s (r). In the case of 
the phonons with the spectrum e(p) = cp, where c is the 
sound velocity, the energy-momentum relation is thus 

(£-p-v s (r)) 2 = C y . (75) 

Eq. (75) can be written in the general Lorentzian form 
with p^ = {-E,p): 

g^PvPv = (76) 

where the metric is 

<7°° = 1, 9 0i = -vl, g lk = -c 2 P k + ■ (77) 

We use the convention to denote indices in the 0-3 range 
by Greek letters and indices in the 1-3 range by Latin 
letters. Thus the dynamics of phonons in the presence 
of the velocity field is the same as the dynamics of pho- 
tons in the gravity field (Unruh 1976): both are described 
by the light-cone equation ds = 0. The interval ds for 
phonons is given by the inverse metric which deter- 
mines the geometry of the effective space: 

ds 2 = g^dx^dx" , (78) 

where x = (t, r) are physical (Galilean) coordinates in 
the laboratory frame. 

A similar relativistic equation holds for the fermionic 
quasiparticles in superfluid 3 He-A in the linear approxi- 
mation close to the gap nodes. In general, i.e. far from 
the gap nodes, the spectrum of quasiparticle in 3 He-A is 
not relativistic: 

e 2 (p)=v 2 (p-pp) 2 + ^(lxp) 2 . (79) 
Pf 

Here vp{p — Pf) is the quasiparticle energy in the normal 
Fermi liquid state above the transition, with pp the Fermi 
momentum and vp — pp/m*; m* is the effective mass, 
which is of order of the mass of the 3 He atom; Aa is 
the so-called gap amplitude and the unit vector 1 points 
in the direction of the gap nodes. 

The energy in Eq. (79) is zero at two points p = eA 
with A = ppl and e = ±1. Close to the two zeroes 
of the energy spectrum one can expand the equation 
(E — p ■ v s (r)) 2 = e 2 (p) in p — eA and write it in a form 
similar to the propagation equation for a massless rela- 
tivistic particle in curved spacetime in the presence of an 
electromagnetic vector potential: 

g^{p ll -eA lx ){p v -eA v ) = Q. (80) 

Here A n = pp(l ■ v s ) and the metric is anisotropic with 
the anisotropy axis along the 1-vector: 
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(81) 



The quantities cm and cj_ correspond to the speeds of 
light propagating along or transverse to 1. 

For simplicity, let us turn back to the case of phonons 
and vortices in 4 He-II which is described by Eqs. (76), 

(77) . If the velocity field is generated by one vortex with 
v quanta of circulation, v s = vncj)/2 , nr 1 then the interval 

(78) in the effective space, where the phonon is propa- 
gating along geodesic curves, becomes: 



ds 2 = 1 



dt + 



2ir{c 2 -v 2 ) / 
dr 2 dz 1 



(82) 



The origin of the Iordanskii force lies in the scattering 
of quasiparticles for small angles, so large distances from 
the vortex core are important. Far from the vortex v 2 /c 2 
is small and can be neglected, and one has 



ds 2 



d4> 

dt + — 



-(dz 2 +dr 2 + r 2 dfi 2 ), u 



2ttc 2 



(83) 



The connection between time and the azimuthal angle 
<p in the interval suggests that there is a characteris- 
tic angular velocity u>. A similar metric with rotation 
was obtained for the so-called spinning cosmic string in 
3 + 1 space-time, which has the rotational angular mo- 
mentum J concentrated in the string core, and for the 
spinning particle in the 2+1 gravity (Mazur 1986, 1996; 
Staruszkievicz 1963; Deser et al. 1984): 



ds 2 = [dt + 



1 



- -{dz z + dr z +r 2 d^), u = 



AJG 

(84) 



where G is the gravitational constant. This gives the fol- 
lowing correspondence between the circulation vk around 
the vortex and the angular momentum J of the spinning 
string 



kv = 8nJG 



(85) 



Although we here consider the analogy between the spin- 
ning string and vortices in 4 Hc-II, there exists a gen- 
eral statement that vortices in any superfluid have the 
properties of spinning cosmic strings (Davis and Shel- 
lard 1989). In particular, the spinning string generates 
a density of angular momentum in the vacuum outside 
the string (Jensen and Kuvcera 1993). The density of 



angular momentum in the superfluid vacuum outside the 
vortex is also nonzero and equals at T = 



r x pv s = Hvubz, 



(86) 



where ijb is the density of elementary bosons in the su- 
perfluid vacuum: the density p/m^ of He atoms in su- 
perfluid 4 Hc-II or the density p/2m 3 of Cooper pairs in 
superfluid 3 He. 



3.2.2. Gravitational Aharonov-Bohm effect 

For the spinning string the gravitational Aharonov- 
Bohm effect is a peculiar topological phenomenon (Mazur 
1986) which can be modeled in condensed matter. On the 
classical level the propagation of particles is described 
by the relativistic equation ds 2 = 0. Outside the string 
the space metric, which enters the interval ds, is flat, 
Eq. (84). But there is a difference in the travel time for 
particles along closed paths around the spinning string in 
opposite directions. As can be seen from Eq. (84), this 
time difference is (Harari and Polychronakos 1988) 
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At large distances form the core the same equation is ap- 
proximately valid due to the equivalence of the metrics 
in Eqs. (83) and (84). The asymmetry between the par- 
ticles orbiting in different directions around the vortex 
implies that, in addition to the symmetric part of the 
cross section, 
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where a(6) is a scattering amplitude, there should be an 
asymmetric part of the scattering cross section, 
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The latter is the origin of the Iordanskii force acting on 
the vortex in the presence of a net momentum from the 
quasiparticles. Another consequence of Eqs. (83), (84) 
is displayed on the quantum level: the connection be- 
tween the time variable t and the angular variable <f> in 
Eqs. (83), (84) implies that the scattering cross sections 
of phonons (photons) from the vortex (string) should be 
periodic functions of energy, with the period equal to Hu. 

Calculations which allow us to find both symmetric 
and asymmetric contributions to the scattering of quasi- 
particles in the velocity field of the vortex have been per- 
formed by Sonin (1997) for phonons and rotons in 4 He-II 
and by Cleary (1968) for the Bogoliubov-Nambu quasi- 
particles in conventional superconductors. In the case of 
phonons the propagation is described by the Lorentzian 
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equation for the scalar field <&: g^d^dv $ = 0, with 
from Eq. (77). We are interested in large distances from 
the vortex core. Thus the quadratic terms v 2 /c 2 can be 
neglected and the equation can be rewritten as (Sonin 
1997) 
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This equation maps the problem under discussion to 
the Aharonov-Bohm (AB) problem for the magnetic flux 
tube (Aharonov and Bohm 1959) with the effective vec- 
tor potential A = v s , where the electric charge e is sub- 
stituted by the mass E/c 2 of the particle (Mazur 1987; 
Jensen and Kuvcera 1993; Gal'tsov and Letelier 1993). 
Actually v s plays the part of the vector potential of the 
so called gravimagnetic field (Volovik 2003). Because of 
the mapping between the electric charge and the mass of 
the particle, one obtains the AB expression (Aharonov 
and Bohm 1959) for the symmetric differential cross sec- 
tion 
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This expression was obtained for the scattering of parti- 
cles with energy E in the background of a spinning string 
with zero mass (Mazur 1987, 1996) and represented the 
gravitational AB effect. Note the singularity at — > 0. 

Returning back to vortices, one finds that the analogue 
with the spinning string is not exact. In more accurate 
calculations one should take into account that as distinct 
from the charged particles in the AB effect, the current in 
the case of phonons is not gauge invariant. As a result the 
scattering of the phonon with momentum p and energy 
E from the vortex is somewhat different (Sonin 1997): 
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The algebraic difference between the AB results in 
Eqs. (91) and (92) is {hc/2irE) sm 2 (irE/uj), which is in- 
dependent of the scattering angle e and thus is not impor- 
tant for the singularity at small scattering angles, which 
is present in Eq. (92) as well. For small E the result in 
Eq. (92) was obtained by Fetter (1964). The generaliza- 
tion of the Fetter result for quasiparticles with arbitrary 
spectrum e(p) (rotons in 4 Hc-II and the Bogoliubov- 
Nambu fermions in superconductors) was recently sug- 
gested by Demircan et al. (1995): In our notations it is 
(vk 2 p/&itvq) cot 2 (#/2), where vq = de/dp is the group 
velocity of a quasiparticle. 



3.2.3. Asymmetric cross section of scattering from a vortex 

The Lorentz force, which acts on the flux tube in the 
presence of an electric current, has its counterpart - the 



Iordanskii force, which acts on the vortex in the pres- 
ence of a mass current of the normal component. The 
Lorentz-type Iordanskii force comes from the asymmetric 
contribution to the cross section (Sonin 1997; Shelankov 
1998a, b), which has the same origin as the singularity at 
small angles in the symmetric cross section and leads to 
a non-zero transverse cross section. 

For the phonons in 4 He-II with the spectrum E(p) = 
cp the transverse cross section is (Sonin 1997) 
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At low E <C hw the result becomes classical: a± = 2t:c/uu 
does not contain the Planck constant h. This means 
that in the low energy limit the asymmetric cross sec- 
tion can be obtained from the classical theory of scatter- 
ing. In this case it can be generalized for an arbitrary 
spectrum E(p) of scattering particles (Sonin 1997). Let 
us consider a particle with the spectrum E(p) moving 
in the background of the velocity field v s (r) around the 
vortex. The velocity field modifies the particle energy 
due to the Dopplcr shift, E(p,r) = E(p) + p ■ v s (r). 
Far from the vortex, where the circulating velocity is 
small, the trajectory of the quasiparticle is almost a 
straight line parallel to, say, the axis y, with the dis- 
tance from the vortex line being the impact parame- 
ter x. It moves along this line with almost constant 
momentum p y w p and almost constant group velocity 
dy/dt — vq = de/dp. The change in transverse mo- 
mentum during this motion is determined by the Hamil- 
tonian equation dp x /dt = —dE/dx = —p y dv sy /dx, or 
dpx/dy = —(p/vc)dv sy /dx. The transverse cross sec- 
tion is obtained by integration of Ap x /p over the impact 
parameter x: 
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This result is purely classical: Plancks constant h drops 
out (it enters only via the quantized circulation vk which 
characterizes the vortex). 



3.2.4. Iordanskii force: quantized vortex and spinning string 

The asymmetric part of scattering, which describes the 
momentum transfer in the transverse direction, after in- 
tegration over the distribution of excitations gives rise 
to the transverse force acting on the vortex if the vortex 
moves with respect to the normal component. This is the 
Iordanskii force: 
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It is proportional to the density of mass current P n car- 
ried by excitations (matter). Of all the parameters de- 
scribing the vortex, it depends only on the circulation 
vn. This confirms the topological origin of this force. In 
the case of the equilibrium distribution of quasiparticles 
one has P n = p n v n , where p n and v n are the density and 
velocity of the normal component of the liquid. Thus 
one obtains Eq. (36) . To avoid the conventional Magnus 
force in this derivation, we assumed that the asymptotic 
velocity of the superfluid component of the liquid is zero 
in the frame of the vortex. 

Eq. (95) was obtained using the asymptotic behavior of 
the flow field v s , which induces the same effective met- 
ric (83) as the metric around the spinning string (84). 
We can thus apply this result directly to the spinning 
string. The asymmetric scattering cross-section of rel- 
ativistic particles from the spinning string is given by 
Eq. (93). This means that in the presence of momentum 
from matter the spinning cosmic string experiences a lift- 
ing force, which corresponds to the Iordanskii force in su- 
pcrfiuids. This force can be obtained from a relativistic 
generalization of Eq. (95). The momentum density P n 
of quasiparticles should be substituted by the component 
Tq of the energy-momentum tensor. As a result, for 2+1 
space-time and for low energy E, which corresponds to 
low temperatures T for matter, the Iordanskii force on a 
spinning string moving with respect to the matter is 
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Here u a is the 3-velocity of the string and Tif is the 
asymptotic value of the energy-momentum tensor of the 
matter at the site of the string. Using the Einstein equa- 
tions one can rewrite this as 
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where Rtf is the Riemannian curvature at the location of 
the string. This corresponds to the force which acts on 
a particle with spin J from the gravitational field owing 
to the interaction of the spin with the Riemann tensor 
(Thorne and Hartle 1985; Mino et al. 1997). 

Note that previously we have shown how ideas first de- 
veloped in the cosmological context (such as nucleation 
of bubbles of different broken symmetry phases during 
a rapid phase transition) could be applied to superfluid 
helium. Here we have the opposite case when ideas origi- 
nating in helium physics have found application in other 
systems described by quantum field theory. 

The Iordanskii force has been experimentally identi- 
fied in rotating 3 Hc-B (Bevan et al. 1997a). Eq. (42) 
describes the temperature dependence of the mutual fric- 
tion coefficient d±. Here the term mzC(T) arises from the 
spectral flow force which is discussed in Sec. 3.3.3 and 
vanishes at low temperatures. Thus at low temperatures 
d± w —p-a.jp (Kopnin et al. 1995), where p is the total 
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FIG. 33. Mutual friction parameters dii ± as a function of 
temperature in 3 He-B at 10 bar. The negative sign of d± — 1 
at T/T c ?s 0.4-0.5 constitutes experimental verification for 
the existence of the Iordanskii force. The curves are fits to 
measured data. (From Bevan et al. 1997a). 



density of the liquid. The negative value of d± arises 
entirely from the existence of the Iordanskii force. This 
is in accordance with the experimental data, which show 
that d± does have negative values at low T (Fig. 33). At 
higher T the spectral flow force dominates, which leads 
to the sign reversal of d± . This fact can be interpreted to 
represent experimental verification of the analog of the 
gravitational Aharonov-Bohm effect for a spinning cos- 
mic string. 



3.3. Spectral flow force and chiral anomaly 

3.3.1. Chiral anomaly 

In the standard model of electroweak interactions there 
are certain quantities, like the baryon number Qb, which 
are classically conserved but can be violated by quantum 
mechanical effects known generically as "chiral anoma- 
lies" . (Each of the quarks is assigned Qb = 1/3 while 
the leptons (neutrinos and the electron) have Qb = 0.) 
The process leading to particle creation is called "spec- 
tral flow", and can be pictured as a process in which 
fermions flow under an external perturbation from neg- 
ative energy levels towards positive energy levels. Some 
fermions therefore cross zero energy and move from the 
Dirac sea into the observable positive energy world. 

The origin for the axial anomaly can be seen from the 
behavior of a chiral particle in constant magnetic field, 
A=(l/2)Bxr. A chiral particle we call a particle with- 
out mass, but with spin a [=1/2]. It can be classified 
as a right or left particle, depending on whether its spin 
is parallel or antiparallcl to its momentum. The Hamil- 
tonians for the right particle with the electric charge eR 
and for the left particle with the electric charge eL are 
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FIG. 34. The energy spectrum of a chiral particle in constant 
magnetic fields along z (Landau levels). The plots on left and 
right show spectra for left and right particles, respectively. 
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As usual, the motion of the particles in the plane per- 
pendicular to B||z is quantized in Landau levels. Thus 
the free motion is effectively reduced to one-dimensional 
motion along B with momentum p z . Fig. 34 shows the 
energy spectrum where the thick lines represent the oc- 
cupied negative-energy states. The peculiar feature of 
the spectrum is that, because of the chirality of the par- 
ticles, the lowest (n = 0) Landau level is asymmetric. 
It crosses zero only in one direction: E = cp z for the 
right particle and E — —cp z for the left. If we now ap- 
ply an electric field E along z, particles are pushed from 
negative to positive energy levels according to the equa- 
tion of motion p z — ejiE z {p z — 6lE z ) and the whole 
Dirac sea moves up (down) creating particles and elec- 
tric charge from the vacuum. This motion of particles 
along the "anomalous" branch of the spectrum is called 
spectral flow. The rate of particle production is propor- 
tional to the density of states at the Landau level, which 
is 
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The production rate of particle number n = np> , + and 
of charge Q = hrCr + n\,e\, from vacuum is 
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This is an anomaly equation for the production of par- 
ticles from vacuum of the type found by Adler (1969) 
and by Bell and Jackiw (1969) in the context of neutral 
pion decay. We see that for particle or charge creation, 
without creation of corresponding antiparticles, it is nec- 
essary to have an asymmetric branch in the dispersion 



relations E(p), which crosses the axis from negative to 
positive energy. Additionally, the symmetry between the 
left and right particles has to be violated: eR ^ cl for 
charge creation and e\ ^ e 2 for particle creation. 

In the electroweak model there are two gauge fields 
whose "electric" and "magnetic" fields may become 
a source for baryoproduction: The hyperchargc field 
U{1) and the weak field SU{2). Anomalous zero mode 
branches exist in the core of a Z-string, where quarks, 
electrons and neutrinos are all chiral particles with known 
hypercharge and weak charges. If we consider a process 
in which one electron, two u-quarks and one (i-quark are 
created, then lepton and baryon numbers are changed by 
one unit while electric charge is conserved (Bevan et al. 
1997b). If we sum appropriate charges for all particles 
according to Eq. (100) the rate of this process is 

At? 

nbar = niept = ^2 (~^W ' *W + By • Ey) (101) 

where Np = 3 is the number of families (generations) of 
fermions, B^, and E a w are the colored SU(2) magnetic 
and electric fields, while By and Ey are the magnetic 
and electric fields of the U(l) hypercharge. While a color 
and hypercharge magnetic flux is always present in the 
Z-string core, a color and hypercharge electric field can 
also be present along the string if the string is moving 
across a background electromagnetic field (Witten 1985) 
or in certain other processes such as the decoupling of two 
linked loops (Vachaspati and Field 1994, 1995; Garriga 
and Vachaspati 1995). Thus parallel electric and mag- 
netic fields in the string change the baryonic charge and 
can lead to cosmological baryogenesis (Barriola 1995) and 
to the presence of antimatter in cosmic rays (Starkman 
1996). 

In superconductors and in superfluid 3 He an anoma- 
lous zero mode branch exists for fermions in the core of 
quantized vortices. For electrons in superconductors it 
was first found by Caroli et al. (1964), and for vortices in 
superfluid 3 He by Kopnin (1993). One of the physically 
important fermionic charges in 3 Hc-A, 3 Hc-B and super- 
conductors which, like baryonic charge in the standard 
model, is not conserved due to the anomaly, is linear mo- 
mentum. The spectral flow of momentum along the zero 
mode branch leads to an additional "lift" force which acts 
on a moving vortex. 

The analogy is clearest for the continuous vortex in 
3 He-A, which has two quanta of superfluid circulation, 
v = 2 (Blaauwgeers et al. 2000). This vortex is similar 
to the continuous Z-vortex in electroweak theory: it is 
characterized by a continuous distribution of the order 
parameter vector 1, which denotes the direction of the 
angular momentum of the Cooper pairs (Fig. 35). When 
multiplied by the Fermi wave number fcp = Pf/^, this 
vector acts on the quasiparticles like an effective "electro- 
magnetic" vector potential A = fcpl. The quasiparticles 
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FIG. 35. The order parameter 1-texture in the soft core of the 
continuous Anderson- Toulouse-Chechetkin vortex. 



in 3 He-A, which are close to the gap nodes, are chiral: 
they are either left or right handed (Volovik and Vachas- 
pati 1996). As follows from the BCS theory of 3 He-A the 
sign of the "electric" charge e, introduced in Sec. 3.2.1, 
simultaneously determines the chirality of the fermions. 
This is clearly seen from a simple isotropic example (with 
en = cj_ = c): 
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A particle with positive (negative) e at the north (south) 
pole is left-handed (right-handed). Here a is a Bogoli- 
ubov spin. 

For such gapless chiral fermions the Adler-Bell-Jackiw 
anomaly applies and the momentum ( "chiral charge" ) of 
quasiparticles is not conserved in the presence of "elec- 
tric" and "magnetic" fields, which are defined by 



E = fcpftl, B = fc F V x 1. 



(103) 



In 3 He-A each right-handed quasiparticle carries the mo- 
mentum pr, = pfI (we reverse the sign of momentum 
when it is used as fermionic charge), and a left-handed 
quasiparticle has pl = — Pr- According to Eq. (100) the 
production rates for right and left handed quasiparticles 
are (since e\ = = 1 in this case) 
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As a result there is a net creation of quasiparticle mo- 
mentum P in a time-dependent texture: 
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What we know for sure is that the total linear momen- 
tum is conserved. Then Eq. (105) means that, in the 
presence of the time-dependent texture, the momentum 
is transferred from the superfluid motion of vacuum to 
matter (it to the heat bath of quasiparticles which form 
the normal component). 
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FIG. 36. Mutual friction parameters du and (1 — d±) in 3 He-A 
at 29.3 bar. The curve is a theoretical fit to the data points. 
(From Bevan et al. 1997a) . 



3.3.2. Anomalous force acting on a continuous vortex and 
baryogenesis from textures 

Let us take as an example the simplest model, namely 
the vortex with a continuous soft core structure in su- 
perfluid 3 He-A. The core has the following distribution 
of the unit vector l(r), which points in the direction of 
the point-like gap nodes in the smooth core 



l(r) = z cos rj(r) + f sin rj(r) 
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where z,r,4> are the cylindrical coordinates. Here r](r) 
is a function which ensures that the 1-vector orientation 
changes within the smooth core from 1(0) = — z to l(oo) = 
z. The circulation of the superfluid velocity along a path 
far outside the soft core corresponds to v = 2: § dr ■ v s = 
2k. Such a continuous 1 texture thus represents a doubly 
quantized vortex without singularities. In practice the 
vortex measured in the NMR experiments is subject to an 
external magnetic field (B||z) and the 1-vector at infinity 
is fixed to the (x,y) plane (see Fig. 35 with l(oo) = x), 
but this does not change the topological structure of the 
vortex. 

When the continuous vortex moves in 3 He-A with 
a velocity vl it generates a time dependent 1 texture 
1 = I(r— vl*)- Hence both an "electric" and a "magnetic" 
field from Eq. (103) exist and this leads to "momentoge- 
nesis" . Integration of the anomalous momentum transfer 
in Eq. (105) over the cross-section of the soft core of the 
moving vortex gives an additional force which acts on the 
vortex due to spectral flow: 



ftP = ttHvCqz x (vl — v„ 
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Here z is the direction of the vortex, Co = k F /3ir 2 , and v n 
is the heat bath velocity. Thus we have obtained Eq. (37) 
with a temperature independent parameter C(T) = Cq. 
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Measurements of the mutual friction coefficients in 
3 He-A with the continuous vortex (Bevan et al. 1997a) 
provide experimental verification for the spectral flow 
force. According to Eq. (42) it should be that dj_ = 
(77J3C0 — p n )/p s - The value of m^Co is the total mass den- 
sity p in the normal phase. Its difference from p in the 
superfluid phase is thus determined by the effect of super- 
fluidity on the particle density which is extremely small: 
p - m 3 C ~ p(A /v F p F ) 2 = p(cj_/c||) 2 < p. Thus one 
must have d± w 1 for all practical temperatures. 3 Hc-A 
experiments at 29.3 bar and T > 0.82 T c are consistent 
with this conclusion within experimental uncertainty: it 
was found that |1 — d±\ < 0.005, as demonstrated in 
Fig. 36. 

3.3.3. Anomalous force acting on a singular vortex and 
baryogenesis with strings 

The discussion of spectral flow in the previous sections 
and particularly Eqs. (100) cannot be directly applied to 
the singular vortex structures which are found in 3 He-B 
and superconductors. The reason is that the deviation 
in the magnitude of the order parameter from its equilib- 
rium value in the cores of such vortices create a potential 
well for the core quasiparticlcs. In this well quasiparti- 
cles have discrete energy levels with some characteristic 
separation Hloq instead of continuous spectra (as a func- 
tion of p z ), as was considered above. Thus the theory of 
spectral flow becomes more complicated, but can still be 
constructed (Kopnin et al. 1995; Kopnin 2002). 

The basic idea is that discrete levels have some broad- 
ening H/t, resulting from the scattering of core excita- 
tions by the free excitations in the heat bath outside the 
core (or by impurities in superconductors). At low tem- 
peratures, when the width of levels is much less than 
their separation, i.e. loqt 3> 1, spectral flow is essentially 
suppressed, C(T) = 0. In the opposite case, lo^t <C 1, the 
levels overlap and we have a situation similar to spectral 
flow in a continuous spectrum: C(T) w Cq. We may con- 
struct an interpolation formula between these two cases: 

C(T) ~TT=iF- (108) 

In fact, both the d\\ and d± mutual friction coefficients 
are affected by this renormalization of the spectral flow 
force (Kopnin et al. 1995; Stone 1996): 

rf || -i(l-d ± ) = ^^^tanh^5 • (109) 
11 v ±J p s l + ioj T 2k B T v ' 

Let us derive this equation using the Landau-type 
phenomenological description for fermions in the vor- 
tex core, as developed by Stone (1996). For simplicity 
we consider the axisymmetric vortex core; the general 
case of the asymmetric core is discussed by Kopnin and 



Volovik (1998). The low-energy spectrum of Caroli-de- 
Gennes-Matricon quasiparticles around a vortex contains 
an anomalous branch of fermionic zero modes. In the case 
of the axisymmetric vortex, excitations on this branch are 
characterized by the angular momentum L z 

E(L z ,p z ) = -Lu ( Pz )L z . (110) 

For superconductors with a coherence length £ much 
larger than the inverse Fermi momentum, pp£ ^> 1, the 
electron wavelength is short compared with the core size, 
and the quasiclassical approximation is relevant. The 
quasiclassical angular momentum L z is a continuous vari- 
able; thus the anomalous branch crosses zero as a func- 
tion of L z at L z — and spectral flow can occur along 
this branch between the vacuum states with E < and 
the excited states with E > 0. Such spectral flow occurs 
during the motion of the vortex with respect to the nor- 
mal component, where it is caused by the interaction with 
impurities in superconductors or with the thermal scat- 
tering states in supcrfluids. In the quasiclassical approx- 
imation, the Doppler shifted spectrum of the fermions in 
the moving vortex has the form 

E(L Z , p) = -u> {p z )L z + (v s - v L ) • p ■ (HI) 

Here the momentum p is assumed to be at the Fermi sur- 
face: p = (pf sin 9 sin ip, p-p sin 9 cos ip, pf cos 9) . The az- 
imuthal angle ip is canonically conjugated to the angular 
momentum L z . This allows us to write the Boltzmann 
equation for the distribution function n(L z ,tp) at fixed 
Pz = Pf cos 9: 

d t n - oj d v n - d v ((v s - v L ) • p) d L n 
n(L z ,<p)-n cq (L z ,<p) 

T 

where the collision time t characterizes the interaction 
of the bound state fermions with impurities or with the 
thermal fermions in the normal component outside the 
vortex core. The equilibrium distribution function is: 

n eq (L z ,tp) = 

f{E - (v n - vl) • P) = f(-oJ L z + (v s - v„) • p) , (113) 

where f(E) = (1 + exp(E /T))^ 1 is the Fermi-function. 
Introducing the shifted variable 

l = L z -{w a - v n ) • p/lo , (114) 

one obtains the equation for n(l, p) 

d t n - uj d v n - <9 v ((v n - v L ) • k) din 
n(l,<p)- f(-upl)) 

T 

which does not contain v s . To find the force acting on the 
vortex from the heat bath environment, we are interested 
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in the evolution of the total momentum of quasiparticles 
in the vortex core: 
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It appears that the equation for P(p z ) can be written in 
closed form 



d t P(p z ) - uj z x P(p z ) 
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Here we take into account that J dldin is limited by the 
bound states below the gap A(T) of bulk liquid: above 
the gap A(T) the spectrum of fermions is continuous, ie. 
the inter-level distance ujq = 0. 

In steady state vortex motion one has dtP(p z ) = 0. 
Then, since /(A(T))-/(-A(T)) = tanh(A(T)/2T), one 
obtains the following contribution to the momentum from 
the heat bath to the core fermions due to the spectral flow 
of bound states below A(T) 




■ bsf 



PF 



PF 



dp z P(p z ) 
2tt t 



k , A(T) 
-tanh — x 
4 2T 



" F — ^ f Z 2 {(vL - VnVoT + 5 X (V L - V n )] . (118) 



PF 



27T 1 + W 2 T 2 



The spectral flow of unbound states above A(T) is not 
suppressed, the corresponding ujqt = 0, since the dis- 
tance between the levels in the continuous spectrum is 
u)o = 0. This gives 



FIG. 37. Experimental values of 1 — a' in 3 He-B at 10 bar 
compared with the theoretical result in Eq. (122). The full 
curve is for the theoretical value of A(T) and the broken curve 
is the fit with reduced A. (From Bevan et al. 1997a). 



Pz) = Pp/37r 2 . In the general case the result in Eqs. (120) 
and (121) is complicated since both co and r depend on 
p z . However, if one neglects this dependence and adds 
the Magnus and Iordanskii forces, one obtains Eq. (109) 
for all T. If the temperature is not too high, so that 
tanh A(T)/2T w 1, one obtains Eq. (108) for the rcnor- 
malizcd spectral flow parameter. 

Comparison of Eq. (109) with experiment should take 
into account that neither oj nor t arc known with good 
precision. However, the following combination of mutual 
friction coefficients does not depend on these parameters 
explicitly: 



F us f = 



l-a' = 



d± 



di + a- d^y 
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tanh 
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2fc B T 



(122) 



K 2 f PF dp z 



sin 2 9 1 



tanh^P^zx (v L -v n ) . (119) 



Thus the total nondissipative (transverse) and frictional 
(longitudinal) parts of the spectral-flow force are 
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In the limit lo^t -C 1, the friction force disappears, while 
the transverse spectral flow force is maximal and coin- 
cides with that obtained for the continuous vortex in 
Eq. (107) with the same value of C = (1/4) J m dp z (j>\- 



This equation (Kopnin et al. 1995) is compared to exper- 
imental data on mutual friction in 3 Hc-B (Bevan et al. 
1997a) in Fig. 37. The agreement is excellent in view of 
the approximations in the theory. 

This shows that the chiral anomaly is relevant for 
the interactions of condensed matter vortices (analogues 
of strings) with fermionic quasiparticles (analogues of 
quarks and leptons). For continuous vortices in 3 He- 
A the spectral flow of fermions between the supcrfluid 
ground state (vacuum) and the heat bath of positive en- 
ergy particles forming the normal component (matter) 
dominates at any relevant temperature. For singular vor- 
tices in 3 He-B it is important at T ~ T c and vanishes at 
T«T C . 

It is interesting to note that in a homogeneous ground 
state (vacuum) momentum (fermionic charge) is con- 
served separately in the ground state (vacuum) and in 
the heat bath of excitations (matter). Topological de- 
fects are thus the mediators for the transfer of momentum 
between these two subsystems. The motion of a vortex 
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across the flow changes both the topological charge of the 
vacuum (say, the winding number of superflow in a torus 
geometry) and the fermionic charge (angular momentum 
in the same geometry) . All these processes are similar to 
those involved in the cosmological production of baryons 
and can thus be investigated in detail. 

However, it is important to note that although spec- 
tral flow leads to anomalous creation of fermionic charge 
from vacuum, the total fermionic charge of vacuum plus 
matter remains conserved. There is also no intrinsic bias 
in this process for the direction of charge transfer: from 
vacuum to matter or in the opposite direction. The nec- 
essary conditions for this mechanism to start operating 
are that the system is in a non-equilibrium state and un- 
der the influence of symmetry breaking which biases the 
direction of charge transfer. How such conditions were 
realized in the early Universe is not clear at present - 
if we interpret the positive baryonic charge of matter in 
our Universe to arise from spectral flow. But it is com- 
monly believed that broken P and CP invariances play 
crucial roles. In contrast in the 3 He case the situation is 
completely under control: The relevant symmetry break- 
ing is achieved by applying rotation or a magnetic field, 
while the non-equilibrium state is generated by applying 
external superflow. 

The experimental verification of momentogenesis in su- 
pcrfiuid 3 He can thus be viewed to support these ideas of 
cosmological baryogenesis via spectral flow and points to 
a future where several cosmological problems are mod- 
elled and studied in light of the superfluid 3 He example 
(Volovik 2003). Now let us consider the opposite effect 
which leads to the analogue of magnetogenesis. 

3.4. Analog of magnetogenesis: vortex textures generated 
in normal-superfluid counterflow 

From Eq. (100) it follows that the magnetic field con- 
figuration can absorb fermionic charge. If this magnetic 
field has helicity, it acquires an excess of right-moving 
particles over left-moving particles: 

K-»l)a = ^A-(VxA) . (123) 

The right-hand side is the so called Chern-Simons (or 
topological) charge of the magnetic field. 

The transformation of particles into a magnetic field 
configuration opens the possibility to examine the cos- 
mological origin of galactic magnetic fields from a system 
of fermions. This is the essential step in the scenario de- 
scribed by Joyce and Shaposhnikov (1997). In this model 
an initial excess of right-handed electrons, eR, was as- 
sumed to have been generated in the early Universe. This 
excess would then have survived until the electroweak 
phase transition (at about 10~ 10 s after the big bang) 
at which point anomalous lepton (and baryon) number 



violating processes became efficient enough to erase the 
excess. However, it appears that well before the elec- 
troweak transition an instability developed, where the 
excess of the right electrons transformed into a hyper- 
magnetic field. Then, when the electroweak transition 
took place, it transformed part of the hypermagnetic field 
into electromagnetic field, so that the universe was be- 
stowed with a primordial (electro-) magnetic field. 

We can now discuss the corresponding process in 3 Hc- 
A - the transformation of fermionic charge to magnetic 
field. In our case it is the quasiparticle momentum which 
plays the part of the relevant fermionic charge. The net 
quasiparticle momentum is generated by the relative flow 
of the normal and superfluid components. This fermionic 
charge is transformed via the chiral anomaly to order pa- 
rameter texture which, as we have seen, plays the part 
of the magnetic field. The transformation occurs in the 
form of an instability which transfers the excess in mo- 
mentum to the formation of 1-textures. The process cor- 
responds to the counterflow instability observed in 3 He- 
A, which has been discussed both experimentally and 
theoretically (Ruutu et al. 1997b). Thus the 3 He-A anal- 
ogy closely follows the cosmological scenario described by 
Joyce and Shaposhnikov (1997). 

Let us discuss this instability. In the presence of coun- 
terflow, v = v n — v s , of the normal component of 3 Hc- 
A liquid with respect to the superfluid, the 1-vector is 
oriented along the flow, lo || v. We are interested in 
the stability condition for such homogeneous counterflow 
with respect to the generation of inhomogeneous pertur- 
bations 51, 

l = io + <ri(r,i) , (124) 

keeping in mind that the space and time dependences 
of 61 correspond to "hyperelectric field" E = kpd t Sl and 
"hypermagnetic field" B = fc F V x SI. 

It is important for our consideration that the 3 He-A 
liquid is anisotropic in the same manner as a nematic 
liquid crystal. For the relativistic fermions this means 
that their motion is determined by the geometry of some 
effective spacetime which in 3 He-A is described by the 
metric tensor in Eq. (81). As we have already discussed 
above, in the presence of counterflow the energy of quasi- 
particles is Doppler shifted by the amount p • v. Since 
the quasiparticles are concentrated near the gap nodes, 
this energy shift is constant and opposite for the two gap 
nodes: p • v m ±pp(l ■ v). The counterflow therefore 
produces what would be an effective chemical potential 
in particle physics, which has opposite sign for the right - 
and left-handed particles: 

Mr = -Ml = Pf(1o • v). (125) 
The kinetic energy of the counterflow is 

^kin = ^Vp„||V . (126) 
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Here the density of the normal component is a tensor in 
the anisotropic 3 He-A liquid, and only the longitudinal 
component p n \\ is involved. 

Let us consider the low-temperature limit T <C T c , 
where T c ~ Ao is the superfluid transition temperature. 
Then using the expression for the longitudinal density of 
the normal component of 3 He- A (Vollhardt and Wolflc 
1990) 



2 m 3 
Pn\\ = I" P 

m 3 



k B T 
Ao 



(127) 



and also the 3 He-A equivalent of the chemical potential 
(125) one obtains 



\rntkl^.v? = \^gT^l. 



(128) 



In the last equality an over-all constant appears to be the 
square root of the determinant of an effective metric in 
3 He-A: = l/ c || c ± — "J^f/Ao- ^ n relativistic theo- 

ries the rhs of Eq. (128) is exactly the energy density of 
the massless right-handed electrons in the presence of the 
chemical potential /ir. Thus the kinetic energy, stored 
in the counterflow, is exactly analogous to the energy 
stored in the right-handed electrons. The same analogy 
occurs between the net quasiparticle linear momentum, 
P = p n v, along 1 and the chiral charge of the right elec- 
trons, 



1 

Pf 



P lo • 



(129) 



The inhomogeneity which absorbs the fermionic charge, 
is represented by a magnetic field configuration in real 
physical vacuum and by a 51-texture in 3 He-A. However, 
Eq. (123) applies in both cases, if in 3 He-A we use the 
standard identification A = kpSl. 

Just as in the particle physics case, we now consider 
the instability towards the production of the "magnetic" 
texture due to the excess of chiral particles. This in- 
stability can be seen by considering the gradient energy 
of the inhomogeneous texture on the background of the 
superflow. In the geometry of the superflow, the tex- 
tural contribution to the free energy of the 51-vector is 
completely equivalent to the conventional energy of the 
hypermagnetic field (Volovik 1992) 



grad 




47re cff 



(1 x (V x 5l)f 



g ij g u F ik F jl = F a 



(130) 



Here Fjfc = V^Afe — VfcAj and we again have included 
the effective anisotropic metric in Eq. (81) appropriate 
for 3 He-A. 
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FIG. 38. NMR signal from an array of continuous ATC vor- 
tices, with the frequency shift from the Larmor value plotted 
on the horizontal axis. The large truncated peak on the right 
is generated by the bulk order-parameter texture, while the 
small satellite peak on the left is produced by the continuous 
soft-core vortex textures. The frequency shift of the satellite 
peak from the Larmor frequency indicates the type of vortex, 
while the intensity of the peak is proportional to the number 
of vortices of this particular type. 



It is interesting that the logarithmic factor in the 
gradient energy plays the part of the running coupling 
e~g = ( 1 /37r he) ln(Ao/T) in particle physics, where e c s 
is the effective hyperclcctric charge; while the gap am- 
plitude Ao, the ultraviolet cutoff, plays the part of the 
Planck energy scale. Now if one has the counterflow in 
3 He-A, or its equivalent - an excess of chiral charge pro- 
duced by the chemical potential p^ - the anomaly gives 
rise to an additional effective term in the magnetic en- 
ergy, corresponding to the interaction of the charge ab- 
sorbed by the magnetic field with the chemical potential. 
This effective energy term is: 

1 



Fes = {nR - n h )p R = ^MR A ■ (V x A) 

= !^-p(io-v)(ji-vx5i) . 

Ira 



(131) 



The rhs corresponds to the well known anomalous inter- 
action of the counterflow with the 1-texture in 3 He-A, 
where p is the mass density of 3 He (Volovik 1992) (the 
additional factor of 3/2 enters due to nonlinear effects). 

For us the most important property of this term is 
that it is linear in the derivatives of 51. Its sign thus can 
be negative, while its magnitude can exceed the positive 
quadratic term in Eq. (130). This leads to the helical 
instability where the inhomogeneous (51-field is formed. 
During this instability the kinetic energy of the quasi- 
particles in the counterflow (analog of the energy stored 
in the fermionic degrees of freedom) is converted into the 
energy of inhomogeneity V x 61, which is the analog of 
the magnetic energy of the hypercharge field. 
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FIG. 39. The height of the satellite peak, which is generated 
by continuous vortex lines, is here monitored during a slow 
increase of fi. Initially there are no vortex lines in the sample 
and the peak height is zero. When the velocity of the counter- 
flow v in the lo direction (corresponding to the chemical po- 
tential /iR, of the chiral electrons) exceeds some critical value, 
a textural instability occurs and the texture transforms into a 
configuration with a lower critical velocity. Simultaneously a 
larger number of vortex lines is formed. This means that the 
center of the sample cylinder becomes suddenly filled with the 
1-texture (or hypermagnetic field) which forms a central clus- 
ter of continuous vortex lines. (From Ruutu et al. (1996c)). 



anomaly is an important mechanism in the interaction of 
vortex textures (the analogue of the hypercharge mag- 
netic fields and cosmic strings) with fermionic excita- 
tions (analogue of quarks and leptons). These two ex- 
periments verified both processes which are induced by 
the anomaly: the nucleation of fermionic charge from 
vacuum, observed by Bevan et al. (1997b), and the in- 
verse process of the nucleation of an effective magnetic 
field from the fermion current as observed by Ruutu et al. 
(1997b). 

3.5. Vortex mass: chiral fermions in strong magnetic field 

Until now we have assumed that the mass My of a vor- 
tex can be neglected in experimental vortex dynamics. 
How well justified is this approach? The term Myd t VL 
in the force balance equation for the vortex contains the 
time derivative and thus at very low frequencies of vor- 
tex motion it can be neglected compared to other forces, 
which are directly proportional to vl . To be more quanti- 
tative, we estimate the vortex mass in the BCS superflu- 
ids and superconductors. There are several contributions 
to the vortex mass. It will also become clear that these 
are related to some most peculiar phenomena in quan- 
tum field theory. We start from the contribution, which 
is relevant for the Bose superfluid 4 He. 



When the helical instability develops in 3 He-A, the fi- 
nal result is the formation of a 1-texture, which represents 
a large jump from the vortex-free counterflow texture to- 
wards the global free energy minimum, the equilibrium 
vortex texture. The new texture contains a central part 
with a periodic 1-texture, where the elementary lattice 
cell represents a so-called Anderson- Toulouse-Chechetkin 
(ATC) continuous vortex (Fig. 35). ATC vortices give 
rise to a characteristic satellite peak in the NMR ab- 
sorption spectrum where their number is directly pro- 
portional to the height of the satellite peak, Fig. 38 
(Blaauwgeers et al. 2000). 

Experimentally this is observed when the rotation ve- 
locity is slowly increased, starting from a state with 
vortex-free counterflow (ie. with fermionic charge, but 
no hypermagnetic field). The helical instability is then 
observed as a discontinuity, when the vortex satellite is 
suddenly formed (Fig. 39). Its peak height jumps from 
zero to a magnitude which approaches that in the equi- 
librium vortex state, which means that the counterflow is 
essentially reduced and the number of vortex lines is close 
to that in equilibrium. Most of the counterflow (fermionic 
charge) thus becomes converted into the vortex texture 
(magnetic field). 

Together with the results of Bevan et al. (1997b), 
the textural superflow instability shows that the chiral 



3.5.1. "Relativistic" mass of the vortex 

In the hydrodynamic theory the mass of the vortex is 
nonzero owing to the compressibility of the liquid. Since 
sound propagation in fluids is similar to light propaga- 
tion in vacuum, the hydrodynamic energy of the vortex, 
soliton or other extended object moving in the liquid is 
connected with its hydrodynamic mass per unit length by 
the "relativistic" equation E — ms 2 , where the speed of 
sound s substitutes the speed of light (Davis 1992; Duan 
1995; Kao and Lee 1995; Iengo and Jug 1995; Wcxlcr 
and Thouless 1996). Thus the hydrodynamic mass of a 
vortex loop of length L at T = is according to Eq. (18): 

M compr - — = In j • (132) 

For Fermi supcrfluids s is on the order of the Fermi ve- 
locity vp ~ Pf/ 171 ( m is the mass of the electron in met- 
als or of the 3 He atom), and the estimate for the hy- 
drodynamic mass of a small vortex loop of length L is 
-Mcompr ~ pa 2 LhiL/^, where a is the interatomic dis- 
tance. For superfluid 4 He, where the core size £ <~ a, we 
can roughly speaking associate the hydrodynamic vortex 
mass <~ pa 2 L with the "mass of the liquid in the vortex 
core volume" . However for 3 He superfluids and for super- 
conductors, where £ 3> a, the hydrodynamic vortex mass 
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is much less than pt; 2 L and other contributions become 
more important. 



3.5.2. Contribution from bound states to the mass of a singular 
vortex 

It appears that the most important contribution to 
vortex mass originates from the quasiparticlcs occupying 
the bound states in the vortex core and thus forming the 
normal component concentrated in the core. For the vor- 
tices in conventional superconductors and in 3 He-B this 
contribution to the vortex mass depends on ujqt and in 
the clean-limit case it is proportional to the mass of the 
liquid in the vortex core, as was first found by Kopnin 
(1978) for superconductors and by Kopnin and Salomaa 
(1991) for superfluid 3 He-B: M Kop nin ~ pi 2 L. This core 
mass is essentially larger than the logarithmically diver- 
gent contribution, which comes from the compressibility. 
In spite of the logarithmic divergence, the latter contains 
the speed of sound in the denominator and thus is smaller 
by the factor (a/£) 2 <C 1, where a is the interatomic dis- 
tance. The compressibility mass of the vortex dominates 
in Bosc supcrfluids, where the core size is small, £ ~ a. 

According to Kopnin's theory the core mass comes 
from the fermions trapped in the vortex core (Kopnin 
1978; Kopnin and Salomaa 1991; van Otterlo et al. 1995; 
Volovik 1997; Kopnin and Vinokur 1998). This Kopnin 
mass of the vortex can be derived using a phenomeno- 
logical approach. Let us consider the limit of low T in 
the superclean regime loqt 3> 1 (Volovik 1997). If the 
vortex moves with velocity vl with respect to the super- 
fluid component, the fermionic energy spectrum in the 
vortex frame is Doppler shifted and has in Eq. (Ill) the 
form: E = E (q) — p vl, where q represents the fermionic 
degrees of freedom in the stationary vortex. The summa- 
tion over the fermionic degrees of freedom leads to the 
extra linear momentum of the vortex cx vl: 

P = J2p°(- E ) = &( k ' v l)<TO - M K op„i„v L , 



M Ko pnin = i^pi<5(£ ) • 



(133) 



For the axisymmctric vortex, where Eq = —L z uJo(p z ) and 
J2 q = f dL z dp z dz/2n, one has 
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(134) 



Eq. (134) can also be obtained from the time dependent 
kinetic equation (112). It is the coefficient in the contri- 
bution to the longitudinal force, which is linear in exter- 
nal frequency u>: Fjj f — — iwMKopninVL- 

Note that this vortex mass is determined in essentially 
the same way as the normal component density in the 



bulk system. The Kopnin vortex mass is nonzero if the 
density of fermionic states is finite in the vortex core, 
which is determined by the inter-level spacing loq in the 
core: iV(0) oc 1/u)q. That gives for the Kopnin vortex 
mass per unit length an estimation: Mxopnin ~ Pp/^o ~ 
p£ 2 . A method to measure this contribution for the 3 He- 
B vortex with spontaneously broken axisymmetry of the 
vortex core, by driving the core into rapid rotation with 
an NMR excitation field, has been suggested by Kopnin 
and Volovik (1998). The much stronger connection to 
the normal component fraction in the core of a vortex 
with continuous structure, like the dominant vortex line 
in 3 Hc-A, will be examined in the next subsection. 



3.5.3. Kopnin vortex mass in the continuous-core model: 
connection to chiral fermions in magnetic field 

The continuous-core vortex in 3 He-A is the best ex- 
ample - a model case - which helps to understand the 
vortex core mass. The continuous-core model can also 
be applied to other Fermi superfluids and superconduc- 
tors: the singular core can thereby be smoothed, so that 
the 1/r-singularity of the superfluid velocity is removed, 
by introducing point nodes in the superfluid energy gap 
in the core region (Fig. 40). As a result the super- 
fluid/supcrconducting state in the vortex core acquires 
the properties of 3 He-A with its continuous vorticity and 
point gap nodes (Volovik and Mineev 1982; Salomaa and 
Volovik 1987). After that one can easily separate dif- 
ferent contributions to the vortex mass. Actually this is 
not only a model: Spontaneous smoothing of the velocity 
singularity occurs in the cores of both types of 3 Hc-B vor- 
tices (Salomaa and Volovik 1987), while in heavy fermion 
and high-T c superconductors such smoothing can occur 
due to the admixture of different pairing states in the 
vortex core. 

For the smoothed singly quantized vortices of 3 He-B 
and superconductors one has two 1- vectors: li and I2, 
each for one of the two spin projections. The simplest 
distribution of both fields is given by Eq. (106) with such 
r?(r) that li(0) = i 2 (0) = -z and li(oo) = -12(00) = ? 
(Salomaa and Volovik 1987). The region of radius R sc , 
where the texture of 1- vectors is concentrated, represents 
the smoothed nonsingular core of the vortex. 

In the continuous vortex, the normal component in the 
core can be considered to be a local quantity, determined 
at each point in the vortex core. This consideration is 
valid for a smooth core with radius R sc » £, where the 
local classical description of the fermionic spectrum can 
be applied. The main contribution comes from the point- 
like gap nodes, where the classical spectrum has the form 

Eq = \J v^(p — pf) 2 + A§(p x i) 2 and A is the gap am- 
plitude. In the presence of a gradient in the 1- field, which 
acts on the quasiparticles as an effective magnetic field, 
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FIG. 40. Smoothing operation of the singular vortex core, (a) 
In the singular vortex the gap is continuously reduced and 
becomes zero exactly on the vortex axis (at r = 0). (b) For 
some vortices it is energetically favorable to avoid the van- 
ishing order parameter at r — 0. Instead, within the smooth 
core, r < R sc , point-like gap nodes appear in the spectrum of 
fermions (Volovik and Mineev 1982). The unit vectors li and 
I2 show the directions to the nodes at different r. Close to 
the gap nodes the spectrum of the fermions is similar to that 
in 3 He-A (example at r = 0). 



this gapless spectrum leads to the nonzero local DOS, 
discussed in Sec. 3.3.1 for the relativistic chiral fermions. 
To apply the DOS in Eq. (99) to the case of anisotropic 
3 He-A, one must make the covariant generalization of the 
DOS, by introducing the general metric tensor and then 
substituting it by the effective 3 He-A metrics which de- 
scribes the anisotropy of 3 He-A. The general form of the 
DOS of the chiral fermions in the curved space is 



N(0) = 



(135) 



where F ik is defined on p. 57. One now can apply this to 
3 He-A where |e| — 1, and the metric tensor is given by 
Eq. (81). Neglecting the dependence of the small velocity 
field v s in the smooth core, one obtains the following local 
DOS for the fermions in the 1-texture at T = 0: 



N(0,v) = 



P% 



2tt 2 A, 



-1(1 x (Vxl)| 



This DOS can be inserted in the expression for the local 
density of the normal component at T = (see Eq. (5.24) 
in the review by Volovik (1990)): 



(Pn)y(r) = W>jpV(0,r) . 



(137) 



For the axisymmetric continuous vortex (Eq. (106)) one 
has 



N(0,v) = 



Pi 



2tt 2 A, 



■ sin?7|9 r ?7|. 



(138) 



The integral of this normal density tensor over the 
cross section of the soft core gives the Kopnin mass of 
the vortex in the local density representation 

M Kop nin - J d 2 r(v L ■ l) 2 p 2 F N(0, r) 

= ^ K - A Jd 2 r^ V (r)\d rV \, (139) 



where vl is the unit vector in the direction of the vortex 
velocity vl . The same equation for the mass can be 
obtained from Eq. (134), using the exact expression for 
the inter- level distance Ldo(p z ) (Volovik 1997). 

Since vf/Aq ~ £ one obtains that the Kopnin mass 
of the continuous vortex <~ p£R sc , i-e. it is linear in the 
dimension R sc of the core (Kopnin 1995). Thus it follows 
that the area law for the vortex mass is valid only for 
vortices with a core size of order £ (i.e. R sc ~ £). 

Note that the vortex mass discussed above comes from 
the normal component trapped in the vortex and which 
thus moves with the vortex velocity v n = vl- In this 
consideration it is assumed that oj t » 1. In this limit 
the normal component in the core and the normal com- 
ponent in the heat bath do not interact with each other 
and thus can move with different velocities. The local 
hydrodynamic energy of the normal component trapped 
by the vortex is 

F = ^(Pn)«(r)(vt - v s ) 4 (v L - v a )j . (140) 

This can be rewritten in a form, which is valid also for 
chiral fermions: 



F = 



Mr + A 
8ir 2 



1\ V'/^/^/v , 



(141) 



where, as before in Eq. (125), the chemical potential of 
the left and right fermions in 3 Hc-A are expressed in 
terms of the counterflow: /lir = — /xl =pf(1-v). Eq. (141) 
represents the magnetic energy of the chiral particles 
with finite chemical potential in a strong magnetic field 



(136) 3.5.4. Associated hydrodynamic mass 



Recently the problem of another vortex mass of hy- 
drodynamic origin was raised by Sonin et al. (1998). It 
is the so-called backflow mass discussed by Baym and 
Chandler (1983), which also can be proportional to the 
core area. Here we compare these two contributions in 
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the superclean regime and at low T -C T c using the model 
of a continuous core. The associated (or induced) mass 
appears when, say, an external body moves in the super- 
fluid. This mass depends on the geometry of the body. 
For the moving cylinder of radius R it is the mass of the 
liquid displaced by the cylinder, 

M associatod = ttR 2 P , (142) 

which is to be added to the actual mass of the cylinder to 
obtain the total inertial mass of the body. In supcrfiuids 
this part of the superfluid component moves with the 
external body and thus can be associated with the normal 
component. A similar mass is responsible for the normal 
component in porous materials, in aerogel for instance, 
where some part of the superfluid is hydrodynamically 
trapped by the pores. It is removed from the overall 
superfluid motion and thus becomes part of the normal 
component. 

In the case when a vortex is trapped on a wire of ra- 
dius R wc 3> £, such that the wire replaces the vortex 
core, Eq. (142) gives the vortex mass due to the backfiow 
around the moving core. This is the simplest realization 
of the backfiow mass of the vortex discussed by Baym 
and Chandler (1983). For this vortex with a wire-core 
the Baym-Chandler mass is the dominant mass of the 
vortex. The Kopnin mass, which can result from nor- 
mal excitations trapped near the surface of the wire, is 
essentially less. 

Let us now consider the Baym-Chandler mass for 
the free vortex at T = 0, using again the continuous- 
core model. In the wire-core vortex this mass arises 
from the backfiow caused by the inhomogcncity in p s : 
p s (r > Rwc) = p and p s (r < R wc ) = 0. Similar but 
less severe inhomogeneity of p s = p — p n occurs in the 
continuous-core vortex due to the nonzero local normal 
density in Eq. (137). Owing to the profile of the local su- 
perfluid density, the external flow is disturbed near the 
core according to the continuity equation 

V-(p s v s ) = . (143) 

If the smooth core is large, R sc 3> £, the deviation of 
the superfluid component inside the smooth core is small 
from its asymptotic value outside the core: 5p s = p — p s ~ 
(t;/R sc )p <C p and can be considered as a perturbation. 
Thus if the asymptotic value of the velocity of the super- 
fluid component with respect of the core is v sU = — vl, 
the disturbance <5v s = of the superflow in the smooth 
core is given by: 

pVH = v i s0 V j (p n ) ij . (144) 

The kinetic energy of the backfiow gives the Baym- 
Chandler mass of the vortex 

M BC = t- / d2r ( V *) 2 • ( 145 ) 

v s0 J 



In the simple approximation that the normal component 
in Eq. (137) is considered isotropic, one obtains 

Mbc = T p J d2rpl[r) ~ ?e ■ (146) 

The Baym-Chandler mass does not depend on the core 
radius R sc , since the large area R 2 C of integration in 
Eq. (146) is compensated by the small value of the normal 
component in the rarefied core, p n ~ p(£/i? sc ). That is 
why this mass is parametrically smaller than the Kopnin 
mass in Eq. (139), if the smooth core is large: R sc 3> £. 

In conclusion, both contributions to the mass of the 
vortex result from the mass of the normal component 
trapped by the vortex. The difference between Kopnin 
mass and Baym-Chandler backfiow mass is only in the 
origin of the normal component trapped by the vortex. 
The relative importance of the two masses depends on 
the vortex core structure: (1) For the free continuous 
vortex with a large core size i? sc ^> £, the Kopnin mass 
dominates: M Kop nin ~ pR sc £, > M B c ~ pS 2 ■ (2) For 
the circulation trapped around a wire of radius i? wc 3> £, 
the Baym-Chandler mass is proportional to the core area, 
M-qc ~ pR wc , and is parametrically larger than the Kop- 
nin mass. (3) For the free vortex core with a core radius 
R ~ £ the situation is not clear since the continuous core 
approximation does not work any more. But extrapola- 
tion of the result in Eq. (146) to R <~ £ suggests that 
the Baym-Chandler mass can be comparable with the 
Kopnin mass. 

3.5.5. Topology of the energy spectrum: gap nodes and their 
ramifications 

Zeroes in the fermionic spectrum, such as Fermi sur- 
faces (surfaces of zeroes) and Fermi points (point zeroes) 
play an extremely important role in condensed matter 
and in analogue models of the low-energy physics of the 
quantum vacuum (Volovik 2003). In condensed mat- 
ter, the gapless fermions interacting with the Bose fields 
of the order parameter lead to the anomalous behavior 
of superfluids and superconductors at low temperatures, 
T < T c , such as spectral flow in vortex dynamics, non- 
analytic behavior of the current and gradient energy, non- 
linear and nonlocal Meissner effect, etc. The counterpart 
of this behavior in high energy physics manifests itself 
in the axial anomaly, baryogenesis, zero charge effect, 
running coupling constants, photon mass, etc. It is the 
zeroes in the fermionic spectrum, through which the con- 
version of the vacuum degrees of freedom into that of 
matter takes place. 

Similar zeroes, but in real space, exist in the cores 
of topological defects, especially in quantized vortices 
(or cosmic strings in high energy physics). Actually the 
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real-space zeroes and the momentum-space zeroes are de- 
scribed by the same topology extended to the 8=4+4 di- 
mensional space. For example from the topological point 
of view, the Fermi-surface represents the vortex singu- 
larity of the Green's function in the iv, p space, where 
uj is the Matsubara frequency. The Green's function 
G(lo,p) = l/(iuj+VF(p—pF)) displays a vortex in the ui,p 
plane with the winding number v = 1. This makes the 
Fermi surface topologically stable and robust under per- 
turbations of the Fermi system. Even though the pole in 
the Green's function can disappear under some pertur- 
bations, the Fermi surface will survive in the marginal 
and Luttinger superfluids. The latter thus belong to the 
same class of Fermi systems as the Landau Fermi-liquid. 

In the same manner, superfluids and superconductors 
with a non-vanishing gap behave in the vicinity of the 
vortex core like superfluid 3 Hc-A in bulk - point-like 
gap nodes appear to be the common feature. Owing to 
the common topological origin of the point nodes, the 
fcrmions near the gap nodes in gapless Fermi superflu- 
ids and superconductors and the low energy fermions lo- 
calized in the cores of vortices in conventional gapped 
superconductors produce similar anomalous effects. 

4. FINAL REMARKS 

We have examined here a few cases in which superfluid 
3 He has been used as an analogue system, to model new 
concepts in quantum field theory The focus in Chap. 2 
is on the question whether the observed vortex forma- 
tion in a rapid 2nd order phase transition is described 
correctly by the Kibble-Zurek scaling theory in 3 He-B. 
In 3 He-A vortex formation (Ruutu et al. 1997b) and vor- 
tex dynamics (Eltsov et al. 2002) are very different. And 
so are the measuring techniques which are required to 
study these features. In fact, the rotating measurements 
on the KZ mechanism in Chap. 2 cannot be repeated in 
3 Hc-A, because vortex-free flow can be reached only up to 
< 1 mm/s while, as seen in Fig. 13, the required thresh- 
old velocity v cn is at least twice higher. No such funda- 
mental limitations restrict calorimetric measurements on 
neutron absorption events in 3 He-A, but present measur- 
ing techniques are not sufficiently developed to tackle this 
question. The consequences from spectral flow and chiral 
anomaly in the vortex dynamics of especially 3 He-A have 
been described in Chap. 3. A number of other analogue 
models involving both 3 He-A and B are described in the 
monograph by Volovik (2003). 

What is the role of analogue studies in examining the 
validity of a new concept in physics? If such a test gives 
a positive result, then good grounds may exist to claim 
that the concept has fundamental value as a physical 
model. However, even if a particular model has been 
shown to work in superfluid 3 He this does not imply uni- 



versal applicability. Analogue studies are not a replace- 
ment for cosmological observations or a means to spare 
investments in high-energy accelerators! For instance, 
from a condensed matter experiment no direct conclu- 
sions can be drawn on the presence or absence of the KZ 
mechanism and its consequences in the Early Universe. 
Although superfluid 3 He experiments in Chap. 2 support 
the validity of the KZ mechanism, at present this sce- 
nario is not a popular explanation for the anisotropy of 
the cosmic background radiation, which was precipitated 
by the inhomogeneity of the Early Universe when it had 
cooled to the point where it became transparent to radi- 
ation. Theories based on the presence of an epoch with 
exponentially accelerated inflationary expansion appear 
to provide better fits. Nevertheless, defects are formed in 
abundance in non-equilibrium phase transitions in con- 
densed matter physics and one is left wondering why 
cosmological transitions should be any different. In this 
sense it remains an open mystery why defects would not 
have formed, or if indeed they were formed, what was 
their role in shaping the Early Universe. 

Irrespective of its cosmological origin, in condensed 
matter physics the KZ mechanism has today become an 
important new concept in rapid non-equilibrium phase 
transitions. It is one of the few intrinsic processes of 
the bulk material which accounts for defect formation 
in a wide spectrum of different systems in the absence 
of boundaries, impurities, or other non-ideal constraints. 
The strongest experimental support to date for the KZ 
mechanism comes from the neutron absorption measure- 
ments in 3 He-B where different types of defects arc ob- 
served to form (usual mass-flow vortices, spin-mass vor- 
tices, AB interfaces). In these measurements it is also di- 
rectly seen that later phase transition fronts, which sweep 
through the system, interact with the existing defects and 
that the final outcome depends on the relative velocity of 
the phase front with respect to the defects. The amount 
of work performed on the KZ mechanism in the form of 
simulations, measurements, and analysis by far exceeds 
other examples of analogue studies. It serves as a good 
example of what can be achieved with analogue studies 
in the unification of physical principles. 
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